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INDUSTRIAL  ROBOTICS:  STATUS  AND  SOME  TRENDS  OF  DEVELDFMENT 

917F0203A  Leningrad  IZVESTIYA  VYSSHIKH  UCHEBNYKH  ZAVEDENIY : 
PRIBOROSTROYENIYE  in  Russian  No  10,  1990  (manuscript  received  9  Jul  90) 
pp  3-8 

[Article  by  A.  M.  Petrina  and  A.  S.  Yushchenko,  Moscow  State  Technical 
University  imeni  N.  E.  Bauman] 

[Text]  A  survey  of  the  current  status  of  industrial 
robotics  in  the  USSR  and  abroad  and  also  analysis  of 
some  trends  of  its  development  is  presented. 

The  use  of  robotics  today  is  one  of  the  main  methods  of  increasing  the 
volume  of  the  products  of  industrial  production,  of  increasing  their 
quality  and,  which  is  especially  important,  of  a  flexible  response  to 
the  variable  requirements  of  the  market.  Industrial  robots  play  a  key 
role  in  the  designs  of  future  plants— automated  plants— being  developed 
today.  Robotics  has  been  transformed  into  a  major  sector  of  industry  in 
the  developed  capitalist  countries.  Such  well  known  corporations  as 
Fanuk,  Fujitsu  (Japan),  General  Electric  and  IBM  (United  States), 
Olivetti  (Italy),  Volkswagon  (West  Germany),  ASEA  (Sweden)  and  also  many 
medium  and  small  enterprises  are  involved  in  research,  development, 
production  and  marketing  of  robotics.  More  than  700  companies, 
including  approximately  400  companies  in  Japan,  50  companies  in  the 
United  States,  and  250  companies  in  Western  Europe,  are  engaged  in 
manufacture  of  industrial  robots  (PR)  in  the  overall  complexity  of 
capitalist  countries.  A  significant  part  of  research  in  robotics  is 
performed  at  the  leading  technical  universities.  Departments  on 
robotics  have  been  created,  for  example,  at  the  Massachusettes  Institute 
of  Technology,  and  Carnegie-Mellon  and  Stanford  Universities.  The 
research  is  financed  mainly  by  industrial  enterprises,  interested  in 
introduction  of  new  technology  into  manufacturing.  However,  significant 
funds  for  development  of  robotics  are  also  allocated  from  the  state 
budget  of  the  leading  capitalist  governments.  All  this  leads  to  steady 
growth  of  both  the  number  of  industrial  robots  used  in  industry  and  to 
an  increase  of  the  quality  characteristics  of  robotic  systems. 

According  to  data  of  the  French  Federation  on  Robotics ,  there  were 

44 , 000  industrial  robots  in  Western  Europe ,  38 , 000  in  the  United  States , 

and  210,000  in  Jap>an  by  the  end  of  1989  [1].  One  should  bear  in  mind 
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the  conditional  nature  of  these  figures,  related  to  different 
definitions  of  industrial  robots  in  different  countries.  However,  the 
trend  of  development  of  industrial  robots  can  be  established.  According 
to  forecasts  of  the  early  1980's,  the  manufacture  and  introduction  of 
robotic  devices  should  be  performed  on  a  considerably  greater  scale  and 
should  have  a  significant  influence  on  their  rate  of  employment  in 
industrial  production.  Thus,  according  to  data  of  the  Institute  of 
Robotics  of  Carnegie-Mellon  University,  industrial  robots  should  free  7 
million  workers  from  production  in  the  United  States.  According  to  the 
forecast  of  the  second  half  the  1980 's,  compiled  by  the  American 
economists  V.  Leontyev  and  F.  Dachin,  no  fewer  than  2  million  workers 
will  be  freed  due  to  robotics  by  the  year  2000.  Other  forecasts  are 
also  subject  to  correction  toward  a  decrease.  This  variation  of  the 
prospects  is  mainly  related  to  the  difficulties  in  introduction  of  the 
manufactured  product.  The  new  forecasts  proceed  from  the  fact  that  the 
rates  of  growth  of  the  worldwide  capitalist  market  of  industrial  robots 
have  stabilized  and  comprise  15-20  percent  annually.  A  period  of 
euphoria  has  been  replaced  by  a  more  sober  attitude  toward  the  aspects 
of  the  application  of  robotics. 

The  forecasts  concerning  the  structure  of  robotic  systems  and  their 
areas  of  application  have  been  largely  confirmed.  As  one  would  assume, 
second-generation  robots,  which  have  external  information  sensors  that 
permit  them  to  correct  their  motions  to  some  extent,  comprised  a 
significant  part  of  the  industrial  robot  stock  by  the  mid-1980's.  These 
are  mainly  robotic  vision  systems  and  force-moment  sensitization.  The 
main  sphere  of  application  of  industrial  robots  still  remains  machine 
building,  primarily,  the  automotive  industry,  and  electronics  and 
electrical  engineering  enterprises.  However,  the  demand  for  robotic 
systems,  designed  for  automation  of  operations  in  nonmachine-building 
sectors,  has  now  increased  sharply.  These  are  construction,  transport, 
the  coal -mining  industry,  light  industry  and  agriculture.  According  to 
data  of  the  Japanese  Association  on  Industrial  Robots,  which  conducted 
research  on  the  needs  of  industry  for  second-generation  industrial 
robots,  the  main  part  of  these  industrial  robots  in  machine  building 
should  be  engaged  in  assembly  operations  (20.9  percent ) .  At  the  same 
time,  32.2  percent  of  industrial  robots  are  required  for  automation  of 
operations  in  extreme  environments,  8.3  percent  are  required  for 
automation  of  transport  operations,  12.4  percent  are  required  for  repair 
and  preventive  maintenance,  and  6.1  percent  are  required  for  use  in 
emergency  situations  [2]. 

The  output  of  industrial  robots  in  the  USSR  during  the  12th  Five-Year 
Plan  should  comprise  100,000  units  [3],  which  is  double  the  data  for  the 
10th  Five-Year  Plan.  Thus,  there  is  the  tendency  to  maintain  the  high 
rates  of  growth  of  the  industrial  robot  stock.  Considerable  experience 
in  robotization  of  the  automotive  industry,  machine  tool  building  and 
instrument  building  has  been  accumulated.  A  total  of  28  percent  of 
industrial  robots  is  now  involved  in  welding  operations,  30  percent 
service  stamping  machines  and  presses,  and  22  percent  service  machine 
tools.  The  operations  are  mainly  automated  due  to  first-generation 
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robots,  i.e.,  rigidly  programmed  robots,  69  percent  of  them  have  cyclic 
control  and  28  percent  have  position  control.  A  total  of  3  percent 
comprise  robots  with  contour  control  that  use  microprocessors  [4].  At 
the  same  time,  this  figure  now  comprises  31  percent  for  the  industrial 
robots  used  in  the  West  and  in  Japan.  Thus,  there  is  some  lag  in  the 
structure  of  the  industrial  robot  stock;  second-generation  robots  have 
not  yet  found  broad  introduction  in  industry.  One  of  the  causes  of  this 
is  the  low  reliability  of  industrial  robots,  manufactured  at  some 
enterprises.  According  to  the  technical  assignment,  it  usually 
comprises  approximately  2,000  hr  of  mean  times  between  failure,  whereas 
this  figure  is  equal  to  10,000  for  Western  industrial  robots.  In 
practice,  the  mean  times  between  failure  is  rarely  less.  Another 
problem  is  the  high  cost  of  robotics.  The  cost  of  loading-transport 
industrial  robots  now  comprises  11,000-15,000  rubles,  while  that  of 
production  industrial  robots  comprises  35,000-  95,000  rubles.  The  level 
of  prices  for  robotics  increases  by  30-40  percent  annually.  Because  of 
this,  the  level  of  the  economic  impact  of  the  use  of  industrial  robots 
in  the  USSR  remains  very  low.  In  the  opinion  of  industrial  workers, 
transition  to  new  management  conditions  may  result  in  a  significant 
decrease  of  plans  for  introduction  of  robotics  at  enterprises.  To  avoid 
this,  a  serious  change  of  the  structure  of  robotic  systems  manufactured 
by  industry  toward  an  increase  of  second-generation  industrial  robots, 
improvement  of  the  reliability  of  industrial  robots,  and  of  their 
efficiency  are  required.  Estimates  show  that  the  necessary  efficiency 
will  be  achieved  with  an  increase  of  the  number  of  realistically  freed 
workers  of  three-fivefold  per  industrial  robot. 

The  worldwide  trend,  related  to  expansion  of  the  traditional  sphere  of 
application  of  robotics  in  the  nonmachine-building  sector  of  the 
national  economy,  has  also  been  manifested  in  Soviet  robotics  f5].  This 
direction  is  of  important  social  significance,  since  the  percentage  of 
unskilled  heavy  manual  labor  is  very  high  in  the  USSR,  compared  to  the 
other  industrially  developed  countries.  This  is  primarily  true  of  the 
coal  and  mining  industry,  ferrous  and  nonferrous  metallurgy, 
construction,  agriculture,  and  the  light  and  food  industry.  The 
importance  of  using  robotics  in  these  sectors  is  determined  by  the  fact 
that  a  large  part  of  the  country’s  labor  resources  is  engaged  in  them, 
while  the  degree  of  automation  is  lower  than  that  in  machine  building. 
Studies  are  being  conducted  intensively  in  this  sector  and  the  first 
models  of  similar  robots  are  beginning  to  be  introduced  in  the  national 
economy. 

Let  us  present  several  examples  in  order  for  the  reader  to  gain  an  idea 
of  the  technical  capabilities  of  modem  industrial  robots  [6,  7],  The 
M5030  Hitachi  welding  robot  (Japan)  with  five  degrees  of  freedom  and 
with  capacity  of  3  kg  performs  welding  at  a  rate  of  2.4  m/min  and  with 
accuracy  of  ±0.2  mm.  The  IRB  2000  ASEA  industrial  arc-welding  robot 
(Sweden)  with  capacity  of  10  kg  has  a  maximum  speed  of  3  m/s  and  a 
positioning  error  of  ±0.1  mm,  and  also  a  control  system  for  12  variables 
simultaneously  (6  are  related  to  the  robot  and  6  are  related  to  external 
devices).  The  KUKA  robot  (West  Germany)  with  capacity  of  150  kg,  which 
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has  seven  degrees  of  freedom,  is  used  for  spot  welding  in  the  automotive 
industry.  An  example  of  a  modern  industrial  robot  for  loading  machine 
tools  and  pallets  is  the  IBM  3000  industrial  robot  with  capacity  of 
30  kg,  with  positioning  error  of  ±0.4  mm,  with  a  speed  of  1.5  m/s,  and 
with  working  zone  radius  of  1.95  m.  Among  second-generation  robots 
designed  for  welding,  one  can  note  the  SCARA  AR-i  270  Stenzel  (West 
Germany)  industrial  robot,  distinguished  by  special  accuracy,  with 
capacity  of  2  kg,  with  speed  of  4.1  m/s,  and  with  positioning  error  of 
±0.05  mm.  The  Sony  Corporation  (Japan)  has  delivered  to  the  market 
two-armed  assembly  robots  with  length  of  the  arms  of  250  and  350  mm, 
capacity  of  1  and  5  kg,  respectively,  and  with  speed  of  5.3  m/s. 

Intel-8086,  Intel-80186,  and  Motorola  68000  microprocessors  are  used  as 
the  control  systems  of  these  robots.  Programming  is  in  high-level  VAL, 
VAL  II,  AL,  AML,  HELP,  JARS,  and  MSL  languages,  and  also  in  universal  C, 
Fortran,  and  Basic  programming  languages. 

The  MP-9S  and  MP-11  (Minavtoprom  [USSR  Ministry  of  the  Automotive 
Industry]),  the  RF-202,  RF-204,  Granat-2.5  and  Granat-10  (USSR  Ministry 
of  the  Radio  Industry]),  Elektronika-NTsIM-O . 1  and  also  the  M20P40.01 
industrial  robots  produced  by  MSPO  [Moscow  Union  of  Consumers  Society] 
Krasnyy  proletariy  among  the  industrial  robots,  manufactured  in  the  USSR 
and  having  characteristics  corresponding  to  the  worldwide  level  can  be 
named;  a  number  of  robots  is  manufactured  in  the  USSR  by  licenses 
purchased  from  foreign  companies,  including  KUKA  and  Fanuc. 

The  control  systems  of  Soviet  industrial  robots  are  implemented  on 
16-bit  K1810  and  K586  microprocessors  and  one-chip  KM1816  computers. 

The  main  computer  is  the  NTs80-01DN  microcomputer,  constructed  on 
KR1801NM2  microprocessors;  K580  microprocessors  are  used  to  control  the 
drives.  A  new  fourth-generation  control  system  IS. 02,  which  is  a 
combination  of  standardized  modules  and  which  is  designed  for  group  and 
autonomous  position-contour  control  of  manipulators  in  a  flexible 
manufacturing  system,  has  been  developed.  This  system  considerably 
exceeds  Soviet  analogues  (MPSU  and  Sf era-36)  in  speed,  capacity  of  the 
control  program  (526  kbytes),  the  number  of  digital  inputs-outputs  and 
reliability.  Its  mean  times  between  failure  is  10,000  hr. 

The  tendency  toward  an  increase  of  the  role  of  the  electric  drive  has 
been  noted  in  the  field  of  industrial  robot  drives,  which  has  begun  to 
be  used  not  only  for  small-  and  medium-capacity  robots  (10-25  kg),  but 
for  large-capacity  industrial  robots  (160-200  kg),  i.e.,  for  industrial 
robots  which  were  earlier  constructed  on  the  basis  of  electrohydraulic 
drives  [8].  Electric  drives  have  been  developed  in  the  USSR,  which  are 
not  inferior  to  foreign  analogues  in  their  parameters.  These  are  the 
EPB-2  drive  with  transistor  converter  and  DVU  and  2DVU  synchronous 
electric  motors,  and  also  the  Razmer  2M  electric  drive  with  transistor 
converter  and  with  asynchronous  motors  of  the  AI  series.  The  most 
promising  for  industrial  robots  are  drives  based  on  contactless  electric 
motors  (DC,  three-phase  asynchronous  motors),  which  have  a  longer 
service  life  than  the  DC  motors  previously  used  here.  Reducerless 
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( "direct  action" )  drives  based  on  high-torque  contactless  motors  must  be 
developed  to  realize  high-performance  industrial  robots,  operating  at 
speeds  up  to  10  m/s.  These  studies  are  being  conducted  both  in  the  USSR 
and  abroad.  The  specifics  of  a  reducerless  drive  consists  in  the  need 
to  construct  an  industrial  robot  control  system  with  regard  to  the 
dynamics  of  the  manipulator  mechanism,  which  in  turn  requires 
construction  of  high-performance  computers  and  real-time  operating 
systems.  Difficulties  arise  in  design  of  a  drive  for  assembly  robots, 
since  in  this  case  the  industrial  robot  should  reproduce  the  given  laws 
of  variation  of  forces,  along  with  reproduction  of  the  given  laws  of 
motion . 

Robotic  vision  systems  and  force-moment  sensitization  are  mainly  used  in 
industrial  robot  sensitization  systems. 

The  robotic  vision  systems  (STZ),  used  under  production  conditions, 
ordinarily  operate  with  images  of  256  x  256  or  512  x  512  elements,  the 
brightness  of  each  of  which  is  encoded  by  6  or  8  bits.  An  essential 
parameter  of  robotic  vision  systems  is  their  speed.  A  large  volume  of 
data  on  the  image  (up  to  256  kbytes),  rigid  time  restrictions  on  the 
permissible  time  of  processing  one  image  impose  their  own  imprint  on  the 
architecture  of  the  hardware  and  on  the  structure  of  the  software  of 
robotic  vision  systems.  The  apparatus  of  a  robotic  vision  system 
usually  consists  of  several  modules,  operating  synchronously  or  in 
parallel ,  performing  the  main  volume  of  image  processing  with  equivalent 
speed  of  10  million  image  elements  per  second.  Examples  of  these 
systems  may  be  the  SERIA  150/151  of  Imaging  Technology,  Inc.  (United 
States)  and  the  DT  2861,  DT  2858,  and  DT  7020  modules  of  Data 
Translation,  Inc.  (United  States) . 

The  indicated  modules  perform  all  the  preliminary  data  processing, 
related  to  viewing  large  files  of  points,  using  point  image  algorithms, 
realized  at  the  hardware  level,  filtration,  calculation  of  brightness 
gradients  and  histograms,  contour  description,  computation  of 
convolution  and  so  on.  The  final  image  description  operation  in  some 
database  and  decision-making  base  is  performed  on  a  personal  computer  or 
on  a  specialized  controller. 

Among  the  robotic  vision  systems  manufactured  in  the  USSR,  one  should 
mention  the  Videosens,  MB- 196  and  Vedoskan  systems,  which,  although 
somewhat  inferior  to  the  above  mentioned  microprocessor  robotic  vision 
systems  in  their  parameters,  process  an  image  within  1-2  seconds  at 
considerably  lower  cost. 

The  force— moment  sensitization  systems  of  industrial  robots , 
manufactured  in  the  USSR  and  abroad  are  completely  functional  high-class 
instruments ,  which  measure  six  projections  of  the  force  vector  and  force 
moment,  conversion  of  analog  data  to  digital  form,  processing  and 
transmission  of  it  to  the  robot  controller  or  to  other  equipment  through 
a  serial  interface  channel. 
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Examples  of  force-moment  systems  may  be:  the  Dinasens-200  NTs 
Robototekhnika  (USSR)  and  the  FS6-120A  of  the  Astek,  Inc.  Company 
(United  States). 

It  should  be  noted  that  the  use  of  real-time  sensory  systems  (for 
example,  during  assembly,  requires  development  of  special  control 
algorithms  and  an  increase  of  the  requirements  on  the  characteristics  of 
the  controllers. 

Among  other  problems  that  require  scientific  studies,  one  should  mention 
those  of  increasing  the  speeds  of  industrial  robots  to  10-15  m/s, 
increasing  the  capacity  to  natural  mass  ratio  of  industrial  robots,  and 
reducing  energy  consumption  and  materials  consumption.  Systematic  study 
of  the  technological  and  economic  feasibility  of  using  robotics  in 
modern  industrial  production  is  necessary  to  improve  the  effectiveness 
of  introducing  industrial  robots. 

The  authors  are  grateful  to  corresponding  member  of  the  USSR  Academy  of 
Sciences  Ye.  P.  Popov,  and  also  to  colleagues  of  the  Scientific  Training 
Center  Robototekhnika  G.  V.  Pismennyy  and  V.  A.  Polskiy  for  assistance 
in  compilation  of  this  survey. 
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MORPHOLOGICAL  ANALYSIS  OF  FLEXIBLE  MANUFACTURING  SYSTEMS 
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[Article  by  0.  M.  Kalin,  Kiev  Polytechnical  Institute] 

[Text]  Analysis  of  the  structural -configuration 
layouts  of  flexible  manufacturing  systems,  based  on 
the  use  of  morphological  Zvicky  analysis,  is 
considered.  Classification  features  of  structural- 
configuration  layouts  of  flexible  manufacturing 
systems  are  determined. 

The  subjective,  empirical  approach  now  predominates  in  design  of 
structural -configuration  layouts  of  flexible  manufacturing  systems 
(GPS),  which  reduces  their  effectiveness  due  to  constriction  of  the  area 
of  searching  for  possible  solutions.  Therefore,  the  problem  of  studying 
processes  of  objective,  formal  structurization  of  manufacturing  systems 
based  on  morphological  methods  is  timely. 

Analysis  of  the  structural -configuration  layouts  of  flexible 
manufacturing  systems,  based  on  the  use  of  morphological  Zvicky  analysis 
[2]  and  further  development  of  it  in  [3,  5].  Two  qualitatively 
different  steps — morphological  analysis  and  morphological  synthesis — are 
distinguished  in  the  morphological  method.  The  step  of  morphological 
analysis  consists  in  many-sided  consideration  of  the  technical  system 
and  description  of  it  in  some  space  of  features.  The  main  problem  of 
this  step  is  careful  study  of  the  design  principles  of  the  system,  of 
the  structural  components  and  different  forms  of  their  interaction  in 
the  system.  The  purpose  of  morphological  analysis  is  development  of 
morphological  classifiers,  which  encompasses  all  the  possible  structural 
solutions  of  the  technical  system.  The  step  of  morphological  synthesis 
consists  in  analysis  and  selection  of  a  set  of  possible  solutions  of 
those  which  correspond  to  the  advanced  requirements.  The  purpose  of 
morphological  synthesis  is  working  out  optimization  algorithms  that 
permit  one  to  orient  oneself  in  the  space  of  possible  solutions. 
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The  most  acceptable  is  the  method  of  the  morphological  Zwicky  box  among 
the  known  morphological  methods  for  problems  of  formulating  the 
structural-configuration  layouts  of  flexible  manufacturing  systems. 
According  to  this  method,  the  technical  system  is  studied  by 
distinguishing  a  number  of  typical  features,  to  which  formal  variables 
xi  correspond  and  on  which  the  structural  solution  is  dependent. 
Gradations  are  found  for  each  feature,  i.e.,  mutually  exclusive  values 
of  the  same  variable.  The  qualitative  state  of  the  structural  component 
is  reflected  by  an  ordinal  number  of  the  i-th  variable,  the  number  of 
which  varies  from  1  to  ki .  Each  of  the  possible  j-th  realizations, 
corresponding  to  the  i-th  value  of  the  variable,  is  selected  for 
formation  of  the  structural  solution. 

Graphically,  the  process  of  formation  of  structural  solutions  can  be 
represented  in  the  following  manner.  Data  on  the  features  of  variables 
and  of  their  gradations  are  represented  in  the  form  of  a  morphological 

table: 


Structural  solutions  are  formulated  by  interrogating  the  gradations  of 
all  features  and  of  including  them  alternately  in  discrete  sets  that 
characterize  one  or  another  variants  of  the  structural  designs  of  the 
system .  For  example , 


< Xi\  X23,  X32 


Xzn  >  ■ 


It  is  obvious  that  a  complete  structural  solution  will  be  formulated 
only  if  each  feature  is  represented  by  one  of  its  gradations.  The 
dimensionality  of  a  specific  set  of  possible  structural  solutions  is 
generally  determined  by 


n=  hku 

i=  i 
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where  n  is  the  number  of  features  and  ki  is  the  number  of  gradations  in 
the  i-th  feature. 


Construction  of  the  morphological  classifier  is  a  subjective  procedure 
and  is  entirely  dependent  on  the  qualifications  of  the  compiler.  This 
process  can  be  based  on  analysis  of  the  structural-configuration  layouts 
of  existing  flexible  manufacturing  systems  and  those  operated  in 
industry.  Conducting  this  analysis  means  objective  determination  of  the 
main  structural  design  principles  of  the  flexible  manufacturing  system 
by  the  expert  method  and  determination  of  their  quantitative  analysis  in 
the  form  of  the  applicability  factor,  which  characterizes  the  fraction 
of  the  flexible  manufacturing  system  with  given  structural  feature  in 
the  total  set  of  those  studied.  This  approach  guarantees  determination 
of  all  the  structural  features  of  the  flexible  manufacturing  system  and 
provides  a  panorama  of  their  use  in  practice. 


The  production  design  characteristics  of  manufacturing  objects  affect  to 
a  considerable  degree  the  selection  of  one  or  another  principle 
structural -configuration  layout  of  a  flexible  manufacturing  system. 
Specifically,  the  design,  overall  dimensions  and  shape  of  objects  affect 
the  selection  of  the  automated  transport -warehousing  system  of  the 
flexible  manufacturing  system.  Taking  this  characteristic  into  account, 
the  entire  combination  of  studied  flexible  manufacturing  systems  must  be 
divided  into  two  groups  upon  morphological  analysis:  flexible 
manufacturing  systems,  designed  to  manufacture  objects  of  prismatic 
shape  (housing  parts) — GPS(Pr) — and  flexible  manufacturing  systems, 
designed  to  manufacture  parts  of  the  body  of  revolution  type — GPS (Tv). 
The  applicability  of  the  structural  solutions  fjjp  and  f^ ,  respectively, 

is  determined  within  the  limits  of  each  of  the  groups. 


The  following  classification  features  are  determined  as  a  result  of  the 
structural -configuration  solutions  of  more  than  200  flexible 
manufacturing  systems,  data  on  which  are  presented  in  [1,  7,  8]: 
xi — the  type  of  flow  scheme  of  blanks,  X2 — the  principle  of  organization 
of  accumulation  of  production  objects,  X3 — the  principle  of  arrangement 
of  the  main  production  equipment  with  respect  to  the  flow  of  production 
objects,  X4 — the  principle  of  organization  of  the  automatic  checking 
system,  X5 — the  principle  of  organization  of  the  automated  tool  support 
system,  and  xg — the  principle  of  organization  of  the  automated  waste 
removal  system. 

The  given  classification  features  fully  characterize  the  technological 
system  of  a  flexible  manufacturing  system,  which  includes,  according  to 
GOST  [State  Standard]  26228-85,  an  automated  transport-warehousing 
system,  automatic  checking  system,  tool  support  system  and  waste  removal 
system.  The  problem  subsequently  reduces  to  study  of  the  set  of 
structural-configuration  layouts  of  the  flexible  manufacturing  system 
according  to  the  adopted  classification  features,  determination  of  the 
gradations  of  these  features,  and  also  according  to  the  quantitative 
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analysis  of  the  applicability  of  structural  solutions  by  coefficient  f^ 
and 


The  results  of  the  analysis  are  presented  in  Tables  1-6.  It  must  be 
noted  that  analysis  of  the  principles  of  organization  of  automated  tool 
support  system  (Table  5)  and  of  a  waste  removal  system  (Table  6)  showed 
that  the  class  of  manufacturing  objects  does  not  affect  the  organization 
of  these  systems;  therefore,  the  analysis  was  performed  for  a  general 
combination  of  flexible  manufacturing  systems  without  separation  into 
GPS(Pr)  and  GPS (Tv) . 

The  analysis  permitted  us  to  formulate  a  morphological  table  of  the 
principal  structural -configuration  layouts  of  a  flexible  manufacturing 
system,  which  includes  the  entire  set  of  the  possible  structural- 
Qonf iguration  layouts  of  flexible  manufacturing  systems.  Performing 
alternate  interrogation  of  the  gradations  of  all  the  structural 
features,  one  can  find  the  set  of  variants  of  the  structural  design  of 
the  system 


W=  {<0i} 


(where  i  =  1,  N) ,  each  of  which  is  a  sample  of  the  gradations  of  the 
features  of  variables,  for  example, 


<Br  =  {*13,  *22.  *33,  A2,  *53.  *62}- 


This  structural  solution  is  characterized  by:  the  use  of  a  rectilinear 
broken  type  of  blank  flow,  by  centralized  separate  arrangement  of  blanks 
and  of  machined  parts,  by  location  of  the  main  production  equipment 
under  the  main  transport  path  (overhead  transport  is  used),  by 
implementation  of  automatic  checking,  by  organization  of  developmental 
support  based  on  the  main  transport  system  of  the  flexible  manufacturing 
system  with  the  flow  of  tools  from  a  central  tool  warehouse  (TsIS) 
directly  to  the  machine  tool  without  intermediate  storage  on  a  separate 
section  ( IU) ,  and  removal  of  wastes  on  the  basis  of  the  main  automated 
transport— warehouse  system  of  the  flexible  manufacturing  system. 
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Table  1 .  Types  of  Blank  Flow  Layouts 
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KEY: 


1.  Type  of  layout 

2 .  Representation 

3 .  Linear 

4.  Rectilinear 

5 .  Reversible 

6.  Rectilinear  broken 

7 .  Closed 

8.  Circular 


9  Circular  branched 

10.  Circular  matrix 

11.  Multirow 

12.  Matrix  rectilinear 

13.  Matrix  reversible 

14.  Network 

1 5 .  Radial 


The  problem  of  morphological  analysis  terminates  in  design  of  a 
morphological  classifier.  The  problem  of  morphological  synthesis  can 
generally  be  formally  described  as  [5]. 

The  specific  function  f(tfi)  is  determined  on  set  W,  i.e.,  an  algorithm 
for  calculating  this  function  f(aq)  for  any  uq  6  W  is  constructed. 

Generally  =  {M©i),  k=\~m)  is  a  set  of  criteria,  part  of  which  is 

minimized,  and  part  of  which  is  maximized.  The  problem  of  morphological 
synthesis  consists  in  removal  of 


/o{coi}-»-opt/cDi€=tt7; 
fj((»i)>0,  fp((DiX0, 


where  fo  is  the  optimal  value  of  the  criterion  of  quality,  fj(tfi)  and 
fp(^i)  are  restrictions,  and  j  =  1,  k,  p=k+  1,  m. 

The  following  methods  can  be  used  to  solve  the  problem  of  morphological 
synthesis  with  the  vector  criterion  of  quality:  sequential  analysis  and 
elimination  of  variants  of  design  solutions  [4],  competing  points  [6], 
and  filtration  for  tree-like  design  [5], 
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Table  2.  Principles  of  Accumulation  of  Blanks  and  Parts 
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KEY: 


1.  Principle  of  accumulation 
and  location 

2 .  Representation 

3.  Centralized 

4 .  Common 

5 .  Separate 

6 .  Decentralized 

7 .  Combination 


8.  Accumulation  on  transport 
system 

9 .  Note 

10.  Automated  work  position 

1 1 .  Automated  storage 

12.  Storage  device  common  for 
blanks  and  parts 

13.  Storage  device  separate  for 
blanks  and  parts 


Table  3.  Principles  of  Arrangement  of  Production  Equipment 
With  Respect  to  Blank  Flow 
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KEY: 

1.  Principle  of  arrangement 

2 .  Representation 

3.  Work  positions  arranged  in  direction  from  main  transport  path 

4.  With  common  branching  for  input  and  output  of  blanks  and  parts 

5.  With  separate  branching  for  input  of  blanks  and  output  of  parts 

6.  Work  positions  are  located 

7.  Under  main  transport  path 

8.  On  main  transport  path 
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Table  4.  Principles  of  Organization  of  Automatic  Checking 
System  of  Flexible  Manufacturing  System 
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i. 
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Table  5.  Principle  of  Organization  of  Automated 
Tool  Support  System  of  Flexible  Manufacturing  System 
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(Key  on  following  page) 


16 


KEY: 


1.  Principle  of  organization 

2.  Notation  of  variable 

3.  Value  of  applicability  factors 

4.  Special  automated  tool  support  systems 

5.  Central  tool  storage-machine  tool 

6.  Central  tool  storage-tool  section-machine  tool 

7 .  Organization  of  tool  support  based  on  automated  transport 
warehouse  system  of  flexible  manufacturing  system 

8.  Delivery  of  tool  manually  on  carts 

Table  6.  Principle  of  Organization  of  Automated 
Waste  Removal  System  of  Flexible  Manufacturing  System 
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KEY: 

1.  Principle  of  organization 

2.  Notation  of  variable 

3.  Value  of  applicability  factors 

4.  Special  automated  waste  removal  systems 

5 .  Removal  of  wastes  through  automated  transport  warehousing 
system  of  flexible  manufacturing  system 

6.  Removal  of  tanks  containing  wastes  using  boom  cranes  and  hand 

carts 
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[Text]  Methods  are  proposed  for  synthesis  of 
correction  algorithms  for  control  of  assembly 
processes  and  parametric  correction  of  simulation 
models  of  flexible  manufacturing  systems,  based  on 
hierarchical  Petri  nets.  Problems  related  to 
development  of  effective  methods  of  control  of 
assembly  processes  and  to  analysis  of  them  by 
simulation  are  considered. 

The  assembly  process  as  an  object  of  automation  is  characterized  by 
cyclicity,  hierarchy  of  representation  of  production  processes, 
asynchronism  and  parallelism  in  operation  of  assembly  equipment.  This 
multiplan  layout  and  complexity  of  assembly  manufacturing  determines  the 
need  to  use  methods  of  simulation  in  problems  of  analysis,  design  and 
control.  The  most  effective  means  of  solving  problems  of  simulation  of 
objects  of  a  given  class  is  the  apparatus  of  Petri  nets  (SP).  Practical 
application  of  Petri  nets  requires  specific  modification  of  their 
classical  apparatus.  In  the  given  case,  this  modification  is  directed 
toward  solution  of  the  following  problems:  adequate  description  of  the 
simulated  object,  the  possibility  of  designing  automated  procedures  of 
synthesis  of  control  models  of  assembly  processes,  and  implementation  of 
automated  control  of  simulation  experiments. 

Problem  1.  It  is  suggested  that  a  modification  of  the  Petri  net — a 
hierarchical  painted  net  (IR-net),  which  has  the  properties  of  animation 
and  safety — be  used  as  a  model  of  the  discrete  assembly  process.  The 
positions  of  the  net  in  this  model  are  compared  to  the  operations  of  the 
process,  transfers  are  compared  to  events  of  a  change  of  operations,  and 
the  colors  of  marks  are  related  to  elements  of  material  and  information 
flows.  The  characteristic  feature  of  the  proposed  modification  is 
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construction  of  a  known  correct  (animated  and  safe)  net  with  as  large  a 
number  of  elements  as  desired  from  some  simple  nets,  which  show  the 
properties  of  the  object  of  simulation,  the  correctness  of  which  has 
been  established.  Let  us  call  the  IR-net  a  structure  consisting  of  a 
finite  set  of  correctly  formulated  net  structures  (blocks)  N  =  (Nj ) , 
where 


{Pu  Ti,  Fu  Hu 
Xi,  i|)»,  noO  > 

T  i=T  i*\)T  iu\}T  iS, 

Tiu  =  {tu | VpeP; :  Fi (p,  f«)=0},  Tis—  {/s| Vp<=Pt :  Hi{ts,  p)=  0}, 

on  which  the  operation  of  substitution  of  the  net  block  is  determined. 
This  set  is  partially  ordered  and  graphically  represented  by  a  tree,  to 
the  root  of  which  corresponds  an  elementary  net  or  correctly  formulated 
block.  The  elementary  net  is  a  strongly  connected  automatic  net  of  two 
positions  and  two  transitions  with  initial  distribution  of  all  colors  of 
the  marks  in  one  of  the  positions.  The  substitution  operation  is  an 
operation  of  superposition  of  the  net  structure  of  the  block  on  a 
fragment  of  the  net  with  a  doubler  position,  substitution  of  its  input 
and  output  transitions-sources  tu  and  discharges  ts  of  the  block  and 
removal  of  the  doubler  position  from  the  received  net.  The  block  is 
considered  correctly  formulated  if  substitution  of  it  instead  of  the 
position  of  the  elementary  net  generates  an  animated  and  safe  Petri  net. 
it  is  shown  that  any  Petri  net,  belonging  to  a  class  of  correct  IR-nets, 
is  animated  and  safe. 


a) 
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The  form  of  representation  of  correct  IR-nets  is  determined.  It  is 
suggested  that  six  typical  blocks  that  give  the  main  ratios  of  the 
relationship  of  operations  of  the  manufacturing  process  be  used,  an 
automatic  block  (Figure  1,  a);  an  alternative  block  (Figure  1,  b);  a 
cyclical  block  (Figure  1,  c ) ;  a  parallel  block  with  the  absence  of 
interactions  (Figure  1,  d) ;  a  parallel  block  with  synchronization  by 
event  (Figure  1,  e);  and  a  parallel  block  with  competition  (Figure 
1,  f).  The  type  of  net  structure  of  the  root  of  the  vertex  of  the  tree 
of  the  IR-net  predetermines  the  class  of  the  simulated  processes.  The 
elementary  net  of  the  root  vertex  corresponds  to  cyclical  processes, 
while  a  net  structure,  represented  by  an  automatic,  alternative  or 
cyclic  block  with  introduced  inhibitor  arcs  that  guarantee  its  safety 
corresponds  to  a  conveyer  net. 

Problem  2.  Management  of  the  manufacturing  system  in  terms  of  a  Petri 
net  reduces  to  the  need  to  manage  transitions  of  the  net  as  a  function 
of  the  markings  of  its  positions.  The  problem  of  synthesizing  the 
control  model  consists  in  assigning  the  representation  (PR),  which 
determines  the  order  of  response  of  competing  transitions,  on  a  set  of 
transitions  T  of  a  model  of  the  process  N.  The  control  algorithm  is 
given  in  the  form  of  a  control  net  Ny  =  (N,  X,  f,  Y,  f,  PR),  where  X  is 
the  set  of  states  of  the  control  object,  Y  is  the  set  of  control 
y 

actions,  p:T  -»  2  is  a  function  of  clear  agreement  of  the  set  of 
transitions  and  subset  of  states,  f:/(  -*  Y  is  a  function  of  the  clear 
conformity  of  marking  of  the  net  and  the  set  of  controls,  and 
PR:(P  x  Q)  x  T  -*  {0,  1,  ...}  is  the  priority  function  of  response  of 
transitions  of  the  net . 

Problem  3.  Automated  simulation  (AIM)  in  control  and  design  problems  of 
assembly  plants  is  related  to  implementation  of  direction  runs  of  the 
simulation  model  (IM)  to  find  a  combination  in  the  given  range  of 
variation  of  input  parameters  that  provides  the  best  value  of  some  index 
of  quality  (performance,  cost  and  so  on).  The  conditional  structure  of 
the  AIM  is  represented  by  the  simulation  model  itself  and  by  the 
experiment  control  device — the  search  model  (FM)  (Figure  2). 

Hierarchical  representation  of  the  simulation  model  based  on  the 
proposed  net  structures  guarantees  direct  interaction  with  algorithms  of 
the  search  model  and  the  possibility  of  parametric  correction  of  the 
model. 


Input 

data 


Figure  2 
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The  search  model  algorithm  realizes  some  optimization  method.  Its 
specifics  consists  in  calculating  the  index  of  quality  on  the  basis  of 
the  results  of  the  run  of  the  simulation  model.  The  relationship  of  the 
parameters  represented  in  the  simulation  model,  the  length  of  the  run  of 
the  model  and  other  characteristics  require  correct  and  effective 
organization  of  management  of  the  series  of  experiments.  The 
correctness  of  interaction  of  the  search  and  simulation  models  is 
maintained  by  specially  developed  control  procedures.  Based  on 
decomposition  and  analysis  of  the  structure  of  the  discrete  optimization 
algorithms  realized  in  the  search  model,  methods  of  improving  the 
efficiency  of  automated  decision-making  during  control  of  experiments 
were  worked  out. 

The  proposed  modifications  of  net  methods  and  approaches  to  use  of  them 
permit  an  increase  of  the  effectiveness  of  solving  complex  object 
control  problems,  and  also  design  and  automated  simulation  of  assembly 
plants . 
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[Text]  Implementation  of  contour  control  in  multilink 
articulated  manipulators  based  on  solution  of  the 
inverse  kinematic  problem  is  considered.  Algorithms 
are  proposed  for  solving  kinematic  controls  based  on 
combination  use  of  a  modified  "number  by  number 
method.  The  characteristics  of  hardware  and 
microprogram  support  of  contour  control  based  on  a 
functionally  oriented  processor  in  the  SChPU  4S-02  are 
described. 

When  developing  NP  systems  for  contour  control  of  industrial  robots 
(PR) ,  main  attention  is  devoted  to  algorithms  for  solving  trajectory 
problems,  which  determine  to  a  significant  degree  the  accuracy  of 
reproduction  of  the  trajectory  and  the  speed  of  the  industrial  robot  [1, 
2] .  A  distinguishing  feature  of  NP  systems  of  multilink  manipulators  is 
the  need  to  solve  two  main  problems  during  control  of  motion:  the 
problem  of  interpolation  that  calculates  the  position  of  the  industrial 
robot  grab  in  Cartesian  coordinates  and  the  inverse  kinematic  problem 
(OKZ),  which  calculates  the  angular  coordinates  of  the  industrial  robot 
according  to  the  given  Cartesian  coordinates  of  the  grab.  Whereas  the 
interpolation  algorithms  are  now  well  known  and  are  widely  used  in 
machine  tool  NP  systems  [1,  2],  algorithms  for  solving  the  inverse 
kinematic  problem  in  real  time  have  been  very  weakly  developed. 

Numerous  gradient  and  other  iterative  algorithms  for  solving  the  OKZ  [3] 
contain  a  large  number  of  operations  of  multiplication,  division  and 
computing  the  trigonometric  functions  and,  therefore,  are  used  only  for 
simulation  of  the  motions  of  the  industrial  robot.  In  this  regard, 
preference  is  given  to  analytical  solution  of  the  OKZ  [4]  with 
subsequent  tabular  implementation  of  the  main  trigonometric  functions  m 
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NP  systems  of  5-6-degree  manipulators.  A  disadvantage  of  this  approach 
is  the  large  volume  of  tables  and  at  the  same  time  the  long  cycle  time 
of  solving  the  OKZ  (up  to  75-125  ms). 

Thus,  working  out  analytical  or  tabular-analytical  methods  of  solving 
the  kinematic  equations  of  multilink  manipulators,  which  reduce  the 
volume  of  tables  and  the  time  of  solving  the  OKZ,  and  also  those  adapted 
for  effective  implementation  on  functionally  oriented  processors  (POP), 
performed  within  the  NP  system,  is  promising. 

Mathematical  postulation  of  the  problem  is  dependent  on  the  design  of 
the  industrial  robot.  Let  us  consider  as  an  example  a  system  of 
equations  that  describe  the  kinematics  of  one  of  the  relatively  simple 
robots  of  the  TUR-10  type  (Figure  1): 


A  [L2  cos(Q2)  L3  cos  (Q2  +  Q3)  +  L4  cos  (Q2  +  Q3  +  Q4)  ]  cos  (Q1 ) ; 

Y  =  [ L2  cos  (Q2)  — L3  cos  (Q2  +  Q3)  4-  L4  cos  (Q2  +  Q3+ Q4)  J  sin  (Q 1 ) ; 

Z=L1  +  L2  sin  (Q2)—  L3  sin  (Q2  +  Q3)  +L4  sin  (Q2+Q3+Q4) ;  ( 1 ) 

.4  =  Q2-  Q3-f  Q4;  • 

0  =  Q1; 

T=Q5, 


where  X,  Y,  and  Z  are  the  Cartesian  coordinates  of  the  grab,  A,  0,  and  T 
are  the  orientation  angles  (Euler  angles),  Lj  (j  =  1,  4)  are  the  lengths 
of  the  links,  and  Qi  (i  =  1,  5)  are  the  angles  between  the  links. 


z  I 


Figure  1.  Kinematic  Layout  of  TUR-10  Robot 
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The  expressions  look  so  simple  for  the  orientation  angles  due  to  the 
simplicity  of  the  design  of  the  robot. 

Analysis  of  relations  (1)  shows  that  the  fraction  of  trigonometric 
operations  and  multiplication  operations  is  high  both  when  solving  the 
direct  kinematic  problem  (PKZ)  (the  Cartesian  coordinates  are  determined 
by  the  angles  of  rotation  of  the  links)  and  when  solving  the  inverse 
kinematic  problem.  Thus,  the  number  of  "long"  operations  when  solving 
the  OKZ,  depending  on  the  method  of  solution  and  on  the  design  of  the 
industrial  robot,  fluctuates  in  the  range  of  12-16  transcendental 
functions,  and  30-40  multiplication  and  division  operations  and 
extracting  the  square  root  for  different  types  of  six-link  manipulators 
[3,  4]. 


Ordinarily,  trigonometric  transformations  are  determined  as  standard 
transformations,  and  different  methods  are  used  to  implement  them  [5]. 
At  the  same  time,  as  follows  from  analysis  of  relations  (1),  more 
complex  operations  of  transformation  and  rotation  of  the  coordinate 
system  can  be  distinguished  as  typical  ( POLDEK— transformation  from  a 
polar  to  a  Cartesian  coordinate  system,  DEKPOL— transformation  from  a 
Cartesian  to  a  polar  coordinate  system,  and  ROTATE — rotation  of  a 
Cartesian  coordinate  system) .  The  algorithm  for  solving  the  inverse 
kinematic  problem  can  be  represented  with  regard  to  the  introduced 
operations  in  the  form  presented  in  Figure  2. 


Figure  2.  Algorithm  for  Solving  Inverse 
Kinematic  Problem 
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The  use  of  POLDEK  and  DEKPOL  commands  reduces  the  number  of  "long" 
operations,  which  results  in  a  reduction  of  the  time  of  fulfilling  the 
algorithm.  The  use  of  the  "number  by  number"  method  [2,  5], 
specifically,  Voider  and  Medjit  methods,  is  promising  for  fulfillment  of 
these  operations.  Although  the  given  methods  have  a  disadvantage  in 
speed  to  some  other  methods  of  calculation  of  trigonometric  functions,  a 
significant  advantage  is  achieved  in  the  given  case  due  to  multiple  use 
of  all  the  functional  capabilities  of  the  operations.  For  example,  the 
POLDEK  operation  is  "equivalent"  to  simultaneous  fulfillment:  sin,  cos 
and  two  multiplications. 

The  NP  systems  of  industrial  robots  axe  now  constructed  on  the  basis  of 
universal  16-bit  microprocessors,  the  architecture  of  which  is 
insufficiently  effective  for  solving  a  class  of  trajectory  problems  [1], 
This  is  caused  by  the  small  word  length  and  rigid  fixed  list  of 
commands.  An  effective  method  of  solving  the  performance  of 
general-purpose  microcomputers  is  connection  of  functionally  oriented 
processors.  The  use  of  functionally  oriented  processors,  oriented 
toward  solution  of  only  trigonometric  problems,  is  undesirable  in 
universal  terminal  control  systems  (TSU)  of  modules  of  flexible 
manufacturing  systems,  which  include  not  only  robots,  but  machine  tools 
and  other  auxiliary  equipment. 

The  4S-02  terminal  control  system  has  been  developed  at  the  Leningrad 
Polytechnical  Institute  for  control  of  the  modules  of  a  flexible 
manufacturing  system  that  includes  the  SChPU  Elektronika  NMS- 12402 
control  module  (based  on  the  K1801VM1  microprocessor)  and  a  functionally 
oriented  processor-active  memory  module  (MAP)  [6] .  The  distinguishing 
features  of  active  memory  modules  are  the  increased  word  length  of  the 
operations  block  (48),  the  memory  capacity — 4K  48-bit  words,  hardware 
support  of  realization  of  complex  logic  transformations,  and  dynamic 
microprogramming  that  introduces  specialized  commands. 

The  indicated  commands — POLDEK,  DEKPOL  and  ROTATE — are  realized  at  the 
microprogram  level  for  effective  solution  of  trajectory  problems  in  the 
active  memory  module.  The  iteration  procedures  that  realize  these 
commands  are  reflected  in  the  table.  The  given  relations  are  variants 
of  known  Voider  and  Medjit  algorithms  [5],  modified  for  more  complete 
consideration  of  the  architectural  features  of  the  active  memory  module: 
48-bit  words,  2-bit  shifts  and  other  features  of  the  operations  block. 

The  most  significant  characteristics  of  the  algorithms  are  enumerated 
below.  A  specially  developed  microprogramming  automation  system  SYMON 
was  used  for  microprogram  realization  of  typical  commands  [7]. 

Inclusion  of  the  given  operations  into  the  command  system  of  the  active 
memory  module  enhances  the  level  of  its  internal  language,  while 
microprogram  realization  reduces  the  time  of  fulfillment  of 
trigonometric  operations  by  a  factor  of  2-2.5,  compared  to  program 
realization. 
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Operation 

Recursion  relation 

Initial 

condition 

Result 

DE'KPOL(A.  B) 

Zi+x=2(Zi+UiXt) 

Xt+l-Xt-UiZ{- 2-« 

Ui= — sgn(Zi) 

0,_,=2[0,- 

arclg(2->)] 

(A)-Xt-X 
(Bi =z(,=y 

0c =0 

• 

(A)  =  A'„  =KVAl2 3-t1'2 

(B) =0„  =  arctg(y/Aj 

POLDEK  (A,  B) 

0,+i=0f— Ui  arctg(2~') 
t/i  =  sgn(<t>i) 

Z;-.1=2(Z, -!-£/.  A';) 
X,+}=Xi—UiZr2-i' 

yn=z„-  2~« 

i.\)  =  Xu=R 
Z[,=0 

(B)  =0o=0 

,•  (A)—Xn=KR  cos(0) 
lB)-y,-/C/?sin(0) 

ROTATE  (A,  B) 

0, +l=0,— Ui  arctg(2-') 
Ui=  sgn(0i) 

Z  ,■ + 1 = 2  (Z  i + U  %  X  ,• ) 
Xw-Xr-UtZr  2-« 

yn=z»-2-n 

1 

(A)  =A’o—A 

RC  =  Y0=Y  \ 

.  (B)«Oo=^ 

i 

CA)  =A’„  = 

.\jA  cos(0)— 

— Ksir(0)] 

(B)  =  yn  = 

=K[y  sin(0)  + 
+Asin(0)] 

Note.  K  is  a  Voider  constant  (K  «  1.65);  A,  B— addresses 
of  operands  participating  in  commands,  third  independent 
variable  that  participates  in  ROTATE  command,  is  transmitted 
through  one  of  the  RONs  [not  further  identified]  of  the 
active  memory  module. 

The  characteristic  features  realizing  the  above  operations  in  the  active 
memory  module  reduce  to  the  following: 

1.  The  algorithm  mainly  corresponds  to  a  Medjit  algorithm,  which, 
compared  to  a  Voider  algorithm,  guarantees  greater  accuracy  and  speed  of 
computations,  and  also  has  hardware  support  in  the  active  memory  module 
in  the  form  of  shift  microoperations  by  two  bits  simultaneously. 

2.  Unlike  the  Medjit  algorithm,  one  cycle  was  organzied  for  the  POLDEK 
and  ROTATE  operations:  the  rotations  are  performed,  beginning  with  the 
largest  angles  (as  for  the  Voider  algorithm).  Upon  completion  of  the 
cycle ,  the  result  xs  shifted  by  N  bits  to  the  right.  There  are  no 
losses  of  accuracy,  since  additional  bits  can  be  stored  due  to  the  long 
word  length  of  the  active  memory  module  in  calculation  with  accuracy  up 
to  24  bits. 

3.  The  angles  of  rotation  can  be  given  in  degrees,  radians  or  other 
units  (sensor  codes)  by  switching  by  program  to  one  or  another  table  of 
constants.  This  permits  one  in  many  cases  to  do  away  with  additional 
scaling  of  operands  when  shaping  drive  control  signals  and  "saves" 
five-six  multiplication  operations.  There  are  no  restrictions  on  the 
value  of  the  angles  of  rotation.  The  constants  correspond  to  the  Voider 
algorithm.  A  longer  word  length  permits  one  to  select  the  number  of 
additional  bits,  necessary  for  the  required  accuracy. 
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4.  Program  adjustment  for  the  accuracy  of  the  calculations  is  possible 
and  a  proportional  increase  of  the  time  of  fulfillment  due  to  accuracy 
is  provided. 

5.  There  is  a  K-bit  shift  to  the  value  K  =  P,  where  P  is  ordinarily 
selected  equal  to  N/2,  but  it  can  be  changed  by  program  as  a  function  of 
the  required  accuracy. 

We  note  that  a  long  word  length  of  the  active  memory  module  permits  one 
to  store  the  values  of  coordinates,  squares  and  fourth  powers  of 
coordinates  in  one  cell,  necessary  for  realization  of  operations  in 
blocks  8  and  10  of  the  algorithm  for  the  inverse  kinematic  problem 
without  a  loss  of  accuracy  (see  Figure  2).  A  total  of  three-four  memory 
cells  would  be  required  to  store  the  intermediate  values  when  realizing 
the  algorithm  on  a  general-purpose  microcomputer,  which  would 
considerably  reduce  the  speed  of  performing  the  operations. 

Problems,  similar  to  robot  control  problems,  are  also  encountered  for 

multicoordinate  machine  tools  having  rotary  coordinates.  Whereas 

accuracy  for  performing  16-bit  operations  is  ordinarily  sufficient  for 

robots,  20-24  bits  are  required  for  a  machine  tool.  There  is  a 

nonproportional  increase  of  the  completion  time  with  an  increase  of 

accuracy  when  the  given  algorithms  are  realized  on  a  microcomputer. 

Thus,  the  completion  time  is  T  (16)  =  2.07  ms  with  accuracy  of  16  bits 

3 

and  T  (24)  =  4.1  ms  with  accuracy  of  24  bits  for  the  Elektronika-60 
3 

computer  with  realization  according  to  Medjit  [2,  5].  The  corresponding 
times  are  equal  to  T  (16)  =  0.5  ms  and  T  (24)  =  0.8  ms  for  an  active 

M  M 

memory  module.  The  total  time  of  solving  trajectory  problems  for  a 
manipulator  of  the  type  given  in  Figure  1  comprises:  2.6  ms  for  the 
direct  kinematic  problem  and  4  ms  for  the  inverse  kinematic  problem 
(with  accuracy  of  16  bits).  The  characteristics  for  the  Elektronika-NMS 
12402  microcomputer  are  similar:  10  ms  for  the  direct  kinematic  problem 
and  16  ms  for  the  inverse  kinematic  problem. 

As  indicated  earlier,  speed  is  increased  due  to  the  use  of  an  algorithm 
(see  Figure  2)  containing  a  fewer  number  of  "long"  operations  than 
ordinarily  used  algorithms,  as  a  result  of  using  complex  POLDEK  and 
DEKPOL  commands,  as  a  result  of  using  modifications  of  the  Voider  and 
Medjit  algorithms  that  permit  a  reduction  of  the  time  of  fulfilling  the 
algorithm  due  to  the  long  word  length  of  the  active  memory  module  and 
also  due  to  transition  of  the  algorithms  to  the  microprogram  level  when 
solving  the  direct  and  inverse  kinematic  problems. 

Thus,  realization  of  contour  control  is  feasible  on  the  basis  of  complex 
operations  of  transformation  of  coordinates,  realized  at  the 
microprogram  level  in  a  functionally  oriented  processor,  and  which 
considerably  increase  the  performance  of  a  general-purpose 
microcomputer . 
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PROGRAM-ADAPTIVE  ALGORITHMS  FOR  POSITIONAL  AND  CONTOUR  CONTROL  OF  ROBOT 
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PRIBOROSTROYENIYE  in  Russian  No  10,  1990  (manuscript  received  28  Feb  90) 
pp  24-31 

[Article  by  A.  V.  Timofeyev  and  A.  B.  Godun,  Leningrad  Institute  of 
Aviation  Instrument  Building] 

[Text]  Algorithms  for  adaptive  control  of  a  robot, 
based  on  correction  of  programmed  motions  at  fixed 
parameters  of  the  control  law,  are  considered.  The 
calculating  relations  and  experimental  results  of 
simulation  of  the  developed  algorithms  are  presented. 

The  efficiency  and  reliability  of  traditional  robot  program  control 
systems  are  reduced  significantly  under  indeterminant  production 
conditions  of  operating  them.  This  is  manifested  in  a  decrease  of 
accuracy,  occurrence  of  vibrations,  soft  failures  and  hard  failures. 

One  of  the  effective  methods  of  increasing  the  reliability  and  quality 
of  functioning  of  a  robot  is  the  use  of  adaptation  devices  in  their 
control  system,  which  provide  automatic  adaptation  to  various  types  of 
disturbances  and  noise.  Known  methods  of  adaptation  are  based  on  the 
sue  of  algorithms  for  identification  of  unknown  factors  or  on 
self-adjustment  of  the  parameters  of  the  control  laws.  However,  they 
are  complicated,  which  makes  it  difficult  to  implement  them  in  a  robot 
control  system. 

New  adaptive  control  algorithms,  based  on  correction  of  the  programmed 
motions  (PD)  of  the  robot  as  a  function  of  the  nature  and  level  of 
perturbations  and  noise  occurring  under  real  operating  conditions  of  the 
robot,  are  proposed  in  this  article. 

Postulation  of  problem  of  program-adaptive  control  of  robot.  The 
dynamics  of  a  wide  range  of  robots  is  described  by  a  system  of 
differential  equations  of  type 


A(q,  l)q(r)+b(q,  q,...,  q<r~l\  g)=u+n,  (1) 
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where  q  —  q(t)  is  an  m— dimensional  vector  function  of  the  generalized 
coordinates,  q( ) }  is  its  j-th  derivative  by  time  t,  j  =  2...r,  (  =  ((t) 
is  the  p-dimensional  vector  function  of  the  robot  parameters,  u  =  u(t) 
is  an  m-dimensional  vector  function  of  control  actions,  A(q,  ( ) , 
b(q. . .q* r' 1  * ,  0  are  given  m  x  m  and  m  functions,  respectively,  t  -  Jr(t) 
is  an  m-dimensional  vector  function  of  constant  external  perturbations, 
about  which  it  is  known  only  that  ||*(t)||  <  car,  where  ct  is  some  positive 
value.  The  order  of  r  vector  differential  equations  (1)  is  dependent  on 
the  type  of  motors  used  in  the  robot.  In  the  case  of  torque  motors, 
r  =  2;  when  using  DC  motors,  r  =  3. 

Let  it  be  required  to  switch  the  robot  during  time  T  =  tp  -  tQ  from  a 

given  initial  state  to  the  desired  final  state,  determined  by  the 
following  boundary  conditions 


?(M=9o---?(r_1)(M=?o(r-I);  q{tT)=qT  ...<7(r-,)(*r)=<7r(r  >> 
with  given  accuracy  € .  In  this  case,  the  purpose  of  control  is  given  by 

ll£(Mli<e,  <3> 

here  . E(t)  =  {e(f)(f)}i=o  ;  e(t)=q(t)—q?(t).  In  practice,  only  part  of  the 

combination  of  boundary  conditions  rather  than  all  of  them,  usually 
boundary  and  specific  conditions,  are  frequently  given 

q(to)  =<7o.  q(lT)  ~qr,  (4) 


or 


q(to)=q0,  q(t0)=q0;  q(tT)=qT,  q(tT)=qT, 


(6) 


MMIKe,  lk(/T)ll<e. 


(7) 


Solution  of  this  problem  is  naturally  divided  into  two  steps.  The 

programmed  motion  Qi>  (0  >  [fo.  /r],  which  satisfies  boundary  conditions 

of  type  (2),  (4),  and  (6),  is  synthesized  in  analytical  form.  In  the 
second  step,  the  control  law  that  guarantees  fulfillment  of  specific 
conditions  of  type  (3),  (5),  and  (7)  is  synthesized.  The  control  laws 
that  solve  the  formulated  problems  will  be  called  terminal  laws.  They 
are  primarily  oriented  toward  use  in  robot  position  control  systems. 

The  specific  condition  in  contour  control  systems  generally  has  the  form 


l|£(/)li<e;  h=J70,  *r].  (8) 


If  truncated  boundary  conditions  (4)  and  (6)  are  given,  the  specific 
condition  has  the  form 


lk(OII<e;  £r] 


(9) 


or 


MOIKe,  De(#)IKe;  t<=[t0,  tT], 


(10) 


This  means  that  the  contour  control  laws  guarantee  tracking  of  a  given 
programmed  motion  qp(t)  and  possibly  of  its  derivatives  in  the  sense  of 
(8),  (9),  and  (10)  with  given  accuracy  over  the  entire  range  of  motion 

[t0,  trp]  . 

Solution  of  terminal  (position)  and  contour  control  problems  is 
complicated  by  the  fact  that  in  practice  there  may  be  initial, 
parametric  or  constantly  acting  perturbations.  They  have  an  unfavorable 
influence  on  the  quality  of  control,  and  result  in  a  decrease  of 
accuracy,  in  autooscillations  or  even  in  a  loss  of  stability. 

Adaptation  devices  are  necessary  to  compensate  for  the  indicated 
perturbations.  Effective  methods  of  adaptation,  based  on  identification 
of  unknown  perturbing  factors  or  on  self-adjustment  of  the  parameters  of 
the  control  laws,  are  suggested  in  [1-3]. 

A  new  method  of  adaptation  in  positioning  and  contour  control  systems, 
the  characteristic  feature  of  which  includes  the  fact  that  the 
programmed  motion  is  corrected  instead  of  the  parameters  of  the  control 
law,  is  studied  in  the  article.  This  permits  one  to  use  ordinary 
positioning  and  contour  control  systems  with  standardized  regulators  of 
actuating  drives,  supplementing  them  only  by  programmed  motion 
correction  programs,  which  implement  the  program-adaptive  control 
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algorithms  suggested  below.  For  definiteness,  robots  with  DC  motors 
whose  dynamics  is  described  by  differential  equation  (1)  at  r  =  3  and 
specific  conditions  of  type  (5),  (7)  and  (9),  (10)  are  considered. 

Program-adaptive  position  control  algorithms.  Let  us  take  as  an  example 
the  problem  of  position  (terminal)  control  with  a  specific  condition  of 
type  (5)  or  (7).  The  parametric  optimization  method  of  programmed 
motions,  proposed  in  [2,  3],  can  be  used  for  analytical  synthesis  of 
programmed  motions  with  regard  to  given  boundary  conditions  (4)  or  (6). 
However,  it  requires  solution  of  the  variational  problem,  which 
complicates  analytical  synthesis  of  programmed  motion.  Therefore,  let 
us  consider  simpler  algorithms  for  analytical  synthesis  of  programmed 
motion.  The  desired  programmed  motion  for  boundary  conditions  (4)  is 
described  by  a  linear  function  of  type 

qP  (t)  —  Qo+  IQt—Qo)  (t — to) IT  (11) 


at  t  €  [tQ,  tT]  and  by  a  third-power  polynomial  for  boundary  conditions 

(6) 


Qp  (£)  —  <7o  +  Xl  {t — to)  +X2  ( t — £o)2/2  +  5<3  (t  to)3 /Q. 


(12) 


G  2  •  12 

Here  t^[to,  tT],  Ki  =  <7o.  *2=  —  (7r~  <7o) - f  (7r  +  2?o).  *3  =  yT  X 

X  (<7r  —  7o)  +"T ?<>)• 


Let  us  consider  the  stabilizing  control  law 


u=A(q,  l)  [qP  +  Ti{q—qp)  +T2(q—qp)  +T3(q—qp)]+b  (q,  <7,  <7-  D¬ 


(13) 


in  which  Ti  ,  ^ ,  and  T3  are  m  *  m  matrix  amplification  factors  in  the 
feedback  channels,  selected  from  the  condition  of  providing  the  desired 
nature  of  the  transition  process  (PP).  It  is  sufficient  to  take 
matrices  T  j ,  T2 ,  and  ^  as  diagonal  to  compensate  for  cross-connections 
in  a  closed  system.  By  substituting  (13)  into  (1)  at  r  =  3,  we  find  the 
equation  of  the  transition  process  in  a  closed  system 


E—TE+A  1  ( q ,  g)n, 
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o/o 


where  £=  {e<'>(/)}2t=0,  e=q—qp,  r=(°  0  /  ) 

'  r3  r2  r  J 


The  following  estimate  in  component  form  is  valid: 


ll« (Oil  j  . 
k(OII  \<c 
.Ik’ (Oil) 


&  (£q) 

e(t0) 

e  Uo)  i 


exp  (—7  (t  -  t0 )) 


C]- 


Here  t  >  to ,  7  is  the  module  of  the  most  essential  part  of  eigennumbers 
of  matrix  T,  c  is  a  positive  constant,  dependent  only  on  T,  and  /i  is  a 
positive  number,  ft  <  min  \i(A),  where  1 i (A)  are  the  eigennumbers  of 

i 

matrix  A(q,  () . 


If  the  initial  perturbations  are  equal  to  zero,  i.e.,  if  s(:0)  =e(t0)  =e{t0), 
selection  of  the  amplification  factors  Tj,  ^ ,  and  rg  from  the  condition 


Y^C[i,-ICK6-1 


(14) 


ensures  fulfillment  of  specific  condition  (5)  for  programmed  motion  (11) 
and  of  specific  condition  (7)  for  programmed  motion  (12). 

Let  us  consider  the  case  when  the  initial  perturbations  are  distinct 
from  zero  and  that  their  level  ||E ( to )  | j  <  f o ,  where  eo  is  some  positive 
integer.  Let  us  assume  that  the  constantly  effective  perturbations  are 
ff(t)  =  0.  Let  us  then  select  amplification  factors  Ti,  1^2 ,  and  I^g,  at 
which  the  following  inequality  will  be  fulfilled,  for  the  given  times  of 
motion  T,  accuracy  e  and  for  the  level  of  initial  perturbations  e  o 


.  1  r  n  «0 

r  In  £ 


(15) 


Control  law  (13)  ensures  fulfillment  of  specific  condition  (5)  for 
programmed  motion  (11)  and  of  specific  condition  (7)  for  programmed 
motion  ( 12 ) . 

Let  7r(t)  £  0.  Then  for  given  T,  e,  and  eo,  selection  of  the 
amplification  factors  from  the  condition  (found  on  the  assumption  that 
exp(yr)  >  1  +  ?T) 

Y>  (ceo+cn-'cxT—e  +  y  (e— ce0— cu~I'c-xTj2~+  I'lTz  * 16  ^ 
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guarantees  fulfillment  of  (5)  for  programmed  motion  (11)  and  of  (7)  for 
programmed  motion  (12).  It  may  happen  in  practice  that  the  real  level 
of  initial  perturbations  (for  example,  due  to  the  inaccurate  initial 
setting  of  the  robot)  will  be  greater  than  given  £q.  Specific 
conditions  (5)  and  (7)  are  then  violated  and  adaptation  of  the  control 
system  is  necessary.  In  this  case,  one  could  use  the  self-adjusting 
method  of  the  control  law  by  correction  of  parameters  T|,  T2 ,  and  T3  or 
of  vector  £.  But  this  method  of  adaptation  with  a  large  level  of 
initial  and  external  constantly  acting  perturbations  may  result  in 
violation  of  the  restrictions  on  control.  Because  of  this,  the  control 
law  may  be  unrealizable. 

A  simpler  and  more  effective  adaptation  algorithm  in  these  cases  may  be 
one  based  on  correction  of  the  programmed  motion  at  fixed  structure  and 
parameters  of  control  law  (13).  The  essence  of  this  algorithm  includes 
the  fact  that  the  programmed  motion  is  adjusted  by  formulas  (11)  or 
(12),  into  which  their  true  values  determined  by  the  position  and  speed 

sensor  are  substituted  instead  of  qo  and  qo  >  at  initial  moment  of  time 
to.  This  single  correction  of  the  programmed  motion  is  sufficient  to 
guarantee  fulfillment  of  specific  condition  (5)  or  (6)  at  the  level  of 
constantly  effective  perturbations 


Let  us  denote  by  tfc  some  moment  of  violation  of  specific  condition  (5) 
or  (7)  with  respect  to  programmed  motion  (11)  or  (12),  respectively, 
constructed  as  a  result  of  the  first  correction  of  the  programmed 
motion.  This  violation  may  occur,  for  example,  due  to  too  high  a  level 
of  constantly  acting  perturbations,  which  results  in  violation  of 
inequality  (17).  There  can  be  an  adaptation  to  these  unmonitored 
perturbations  due  to  correction  of  the  programmed  motion.  To  do  this, 
the  programmed  motion  must  be  constructed  according  to  formulas  (11)  or 

(12),  having  substituted  the  values  q(tk)  and  q(tk),  recorded  by  the 

4 

position  and  speed  sensor  at  moment  of  time  tk ,  instead  of  qo  and  qo , 
and  dimensionality  tT  -  tfe  instead  of  T.  One  can  show  that  specific 

conditions  (5)  and  (7)  will  be  fulfilled  after  the  final  number  of 
corrections  of  programmed  motion. 

Program-adaptive  contour  control  algorithms.  Let  us  now  consider  the 
problem  of  contour  control  of  a  robot  with  specific  condition  (9)  or 
(10).  The  necessary  condition  for  fulfillment  of  (9)  and  (10)  is  the 
low  level  of  initial  perturbations  compared  to  the  given  machining 
accuracy  of  the  programmed  motion.  It  is  necessary  that  |Je ( to )  ||  <  fo  < 
<  £  in  the  case  of  specific  condition  (9).  Let  us  select  the  program 
motion  in  the  form  of  (11)  and  let  us  use  control  law  (13).  Specific 
condition  (9)  will  be  fulfilled  if  the  levels  of  constantly  active  and 
initial  perturbations  satisfy  the  requirement 
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ceo  +  cy-'n-'cn^e. 


Otherwise,  it  may  be  violated  at  some  moment  of  time  t^ .  In  this  case, 

one  must  construct  an  auxiliary  (corrected)  programmed  motion  by  formula 

(11),  having  substituted  in  it  the  values  of  q(tfc)  and  t  -  t,  , 

1  K 

respectively,  instead  of  qo  and  T.  The  programmed  motion  corrected  in 
this  manner  is  substituted  into  control  law  (13)  and  we  check 
fulfillment  of  (9)  with  respect  to  the  initial  programmed  motion  (11). 

If  it  is  violated,  the  auxiliary  programmed  motion,  subsequently  used  in 
control  law  (13),  is  again  corrected.  One  can  show  that  the  number  of 
these  corrections  will  be  finite.  The  experimental  results  of 
simulation,  considered  subsequently,  also  indicate  this. 

The  program-adaptive  contour  control  algorithm  is  constructed  in  similar 
fashion,  but  programmed  motion  of  type  (12)  is  used  and  corrected  in 
control  law  (13). 

Results  of  simulation  of  program-adaptive  robot  control  algorithms.  The 
proposed  program-adaptive  control  algorithms,  based  on  correction  of 
programmed  motion,  were  simulated  in  the  problem  of  control  of  an 
industrial  robot  of  type  (TUR-2.5)  with  actuating  drives  based  on  DC 
motors  to  perform  manufacturing  operations  of  laser  machining  of 
materials  and  the  operation  of  grinding  the  edge  of  a  blank  during 
manufacture  of  footwear  [4].  The  adaptation  was  performed  both  with 
respect  to  errors  of  the  initial  setting  of  the  manipulator  and  with 
respect  to  constant  perturbations.  Not  only  control  law  (13),  but  the 
programmed  control  law  of  the  following  type  were  simulated 


up-A(qp,  | )qP  +  b(qp,  qPt  qpy  g),  (18) 


and  the  programmed  control  law  with  linear  feedback  was  simulated  by 
position  and  speed  of  the  following  type 


u=uP+T2(q—qp)+r3(q—qp).  <19> 


Control  of  the  TOR-2. 5  robot  includes  control  of  the  manipulator 
directly  according  to  degrees  of  mobility  qi  and  q2  and  control  of  the 
working  member  (linear  displacement  and  rotation  of  the  grab).  The 
first  case  is  of  the  greatest  applied  interest.  Let  us  present  an 

explicit  form  of  matrix  function  A(q,  ( )  and  vector  function  b(q,  q,  q, 
0 ,  used  in  control  laws  (13),  (18)  and  (19),  to  make  it  more  specific. 
With  regard  to  the  dynamics  of  DC  motors,  we  find 
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A(q,  l)=LaK3-l(Jp  +  p-1i \~lAm(q,  § m)). 

b  ( q ,  q,  q,  1)  =  {UKnrxP~\^Am{q,  ?m)  + 

+  R„Kz-l{Jp+p-'i\-lAm(q,  lm)))q  + 

+ UKr' (KTpq+p-lT\-lbm(q,  q>  I «))  + 

+R,,KjCl(KTpq+p-l'n~1bm{q .  <7,  lm))+Kipq> 


where 


r  an 

~  LQ21 

1al2  sin  ( 
. 

— a,2Si 


Am(q<  %m) 

bm  (q 


ai2cos(t71— q2) 
cos  (?!— 92)  a22 


1  (91 — q2)q22+b\  sin^  +  Kritfi 
>sin(«7i — 92)^i2  +  ^2sin  q2  +  Knq2 
ai\=m\T\2+  (W2  +  ffio)/i2;  <212=^21  =  tn2lir2-{-mol\l2;  fl22~m2r22  +  mo^2~? 
&i  =  (/n2r2  +  m1/i  +  mo/i)  g;  b2  =  (m2r2+ m0l2) g\ 
lm={mu  m2,  r,,  r2,  Zj,  /2,  g,  «o,  Kti,  Ktt}- 


Here  mi ,  li ,  m2,  and  I2  are  the  mass  and  length  of  the  first  and  second 

links  of  the  manipulator,  respectively,  ri  and  r2  are  the  distances  from 
the  centers  of  gravity  to  the  rotational  axes  of  the  first  and  second 
links,  mo  is  the  mass  of  the  weight,  g  is  the  acceleration  of  gravity, 
and  J,  K^,,  p,  7,  L^,  R^,  and  are  the  constant  diagonal  2x2 

matrices  of  the  robot  parameters.  After  the  values  of  7  that  satisfy 
conditions  (14),  (15)  or  (16)  as  a  function  of  the  postulated  problem 
have  been  selected,  the  amplification  factors  Tj  =  diag( 7i  i ( J  ’  )  can  be 
calculated  by  the  given  eigennumbers  7T ( T )  <  -7,  r  =  1-6  and  by  the 

formulas 


Yii(1)—  £  As,'— ;:-M  (r)> 

y..(2) = — A3i-2  (r)  x3;-i  (r)  -hi- 2  (r)  Ui  (H  (H  hi  (H , 

Yii(3>  =  /-3,'— 2  (T)  A3f— 1  (T)hi  (f)  . 


at  i  =  1,  2. 

Let  us  consider  the  results  of  typical  experiments  in  the  problem  of 
position  control  with  specific  condition  (5).  Programmed  motion  was 
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synthesized  in  the  form  of  (11),  where  q(t0)  =  (1,  1;  1,  1),  q(tT)  = 

=  (1,0;  1,0),  T=  1  s,  and  the  required  positioning  accuracy  is  e  = 

=  0.01  rad.  The  values  of  Tj  =  -151,  T2  =  751,  T3  =  -1251  (I  is  a  2  x  2 
unit  matrix),  corresponding  to  l^T)  =  -5,  r  =  1...6,  were  taken  as  the 

amplification  factors  in  (13). 


Figure  1 


Graphs  of  the  variation  of  generalized  coordinates  for  three  types  of 
control  law:  stabilizing,  programmed  and  programmed  with  linear 
feedback  (trajectories  1,  2,  3,  respective) — are  presented  in  Figure  1. 
The  level  of  the  initial  perturbations  e(to)  =  (0.22;  0.22)  considerably 
exceeded  e  due  to  the  errors  of  the  initial  setting  of  the  manipulator, 
the  real  motion  q(t)  lies  entirely  in  the  e-neighborhood  of  the 
corrected  programmed  motion  after  one  correction  and,  accordingly,  (5) 
is  fulfilled.  Similar  results  were  found  in  simulation  of  program- 
adaptive  contour  control  of  a  robot  with  specific  condition  (9)  and 
programmed  motion  (11). 

The  graphs  of  variation  of  generalized  coordinates  for  a  stabilized 
control  law  without  corrections  of  programmed  motions  and  with 
correction  of  programmed  motion  (trajectories  1  and  2,  respectively)  in 
the  presence  of  constant  perturbations  I,  the  level  C7T  of  which  exceeded 
the  permissible  level  (17),  are  shown  in  Figure  2.  The  programmed 
motion  (11)  was  corrected  at  moments  of  perturbation  (9),  and  after  this 
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(four  corrections),  the  specific  condition  (9)  was  fulfilled  all  the 
time.  The  use  of  control  law  (13)  without  corrections  of  programmed 
motion  resulted  in  violation  of  (9)  with  €  =  0.01  rad,  which  indicates 
the  need  to  correct  programmed  motions  at  a  high  level  of  constant 
perturbations . 
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[Article  by  V.  F.  Filaretov  and  A.  F.  Polyanskiy,  Far  Eastern 
Polytechnical  Institute  imeni  V.  V.  Kuybyshev] 

[Text]  The  problems  of  designing  a  two-way  copying 
manipulator  system  with  active  reflection  of  force, 
self-adjusting  upon  variation  of  load  parameters,  are 
considered.  The  regulator  is  designed  according  to 
the  quadratic  criterion  of  quality. 

A  human  operator  performs  the  decision-making  and  planning  functions  of 
the  motion  trajectories  of  the  actuating  member  (10)  in  copying 
manipulators  (KM) ,  and  motions  are  distributed  by  separate  degrees  of 
freedom  directly  during  operation.  Two-way  systems  (SDD) ,  located  in 
each  degree  of  freedom  of  the  copying  manipulator  and  which  are  systems 
of  the  actuating  level  of  control,  should  provide  precise  tracking  of 
the  assigned  trajectories  of  motion,  and  also  reproduction  of  forces  on 
the  given  member,  developed  in  the  degrees  of  freedom  of  the  actuating 
member.  The  dynamic  effects  of  the  mutual  influence  between  the  degrees 
of  freedom  of  the  manipulator  have  a  strong  influence  on  the  quality  of 
operation  of  the  two-way  system  during  motion. 

A  two-way  system  was  designed  in  [1]  with  regard  to  variation  of  the 
load  parameters.  It  was  assumed  that  the  variable  moments  of  inertia 
and  also  the  moments  of  gravitational  forces  have  the  main  influence  on 
the  dynamics  of  the  two-way  system.  However,  besides  inertial  and 
gravitational  forces,  Coriolis  and  centrifugal  forces,  which  can  not  be 
disregarded  at  high  speeds  (more  than  0.5  m/s)  of  the  motion  of  the  grab 
of  the  copying  manipulator  [2],  also  affects  the  links  of  the 
manipulator . 

The  problem  of  designing  self-adjusting  (SN)  two-way  systems  is  posed 
and  solved  in  this  paper  with  regard  to  the  total  dynamics  of  the 
actuating  members  of  the  copying  manipulator.  It  is  assumed  that  the 
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given  member  is  light  and  its  influence  on  the  dynamics  on  any  two-way 
system  is  determined  only  by  the  constant  moment  of  inertia. 

The  moment  effect  of  any  degree  of  freedom  of  the  copying  manipulator  on 
a  two-way  system  is  generally  determined  by  the  expression 


Q  =  H'  (<?i . . .  qN)an{t)  +  h(q\ ...  qx\  qi ...  q n\  <7i  •  •  •  q%-\,  qw  ■  ■  •  <7x)> 


(1) 


where  H' (q^...c^j)  is  the  moment  of  inertia  of  the  kinematic  circuit  of 

the  actuating  member  with  load  with  respect  to  the  considered  degrees  of 

freedom,  h(qr  .  .qN;  is  a  term  determined  by 

the  effect  of  gravitational,  Coriolis  and  centrifugal  forces,  and  also 
by  the  components  of  forces  of  inertia,  which  were  not  included  in  the 
first  term,  q^.-.q^  are  generalized  coordinates  of  the  kinematic  circuit 

of  the  actuating  member,  «H(t)  =  qi ,  i  is  the  number  of  the  considered 

degree  of  freedom,  and  N  is  the  number  of  degrees  of  freedom  of  the 
manipulator . 

Let  us  consider  a  two-way  system,  the  actuating  elements  of  which  are  DC 
motors  with  independent  excitation.  The  dynamics  of  these  motors  in  the 
direction  of  the  operator  and  load,  respectively,  is  described  by  the 
following  differential  equations: 


U\  ( t )  =  (i?l+l*lp)*al  (0  Kmlt'lOon(0» 

Ma Bl  (t)  =  JW.1  (0  —hh^-oa  (0  +Mt1  (aon)-Mon  (0 Ihi  (  2  > 

u2(t)  =  (R2+L2p)is<2{t)  +Ka2haH(t), 

Mto2(t)=K«2i*z(t)=hhaH(t)  +MT2(aH)  +Qlh, 

where  R)  and  R2  are  the  effective  resistance  of  the  armature  circuits  of 
the  motors,  Lj  and  L2  are  the  inductance  of  the  armature  circuits,  KM1 

and  are  the  torque  factors,  and  are  the  coefficients  of  the 

counter ,  emf ,  J  j  and  J  2  are  the  moments  of  inertia  of  the  rotors  of  the 
motors,  of  the  rotating  parts  of  the  kinematic  gears  and  of  the 
corresponding  links,  ii  and  i2  are  the  gear  ratios  of  the  reduction 

gears,  Mxl(aon),  Mt 2(aH)  are  the  dry  friction  moments,  U|(t)  and  u2(t) 

are  the  control  voltages,  i ^ ^ ( t ) ,  and  ^ ^ ^  are  the  currents  of  the 

motor  armatures,  M^^t)  and  M^ft)  are  the  moments  developed  by  the 


41 


motors,  ffon(t)  and  are  the  coordinates  of  the  operator  and  load 

shafts,  MQn ( t )  is  a  moment  developed  by  the  operator  and  p  is  a 
differentiation  symbol. 


According  to  the  methods  suggested  in  [ 1 ] ,  let  us  assume  that  a  ,  a  , 

on’  on 

ffon  and  a  Qn  are  independent  input  variables  in  the  mode  of  induced 

motion  of  the  operator  shaft  of  the  two-way  system.  Taking  this  into 
account  and  also  taking  into  account  relations  (l)-(3),  one  can  describe 
the  dynamics  of  the  two-way  system  in  the  mentioned  mode  with  respect  to 

phase  coordinates  (FK)  -Vi  =  aH(r),  .*2  =  ai r(/),  Xs^aH(t),  x4=  MQn ( t )  by  the 

following  system  of  differential  equations: 


—  x2, 


■T 


x2  =  xz, 


Ky\iK^i  R„  (J2  -r  H)  L2H  ,  ATM2 

■*3  —  -  Lf  (J2  +H)X 2  traPHT)  '  itLz  (J«  —  H)  ‘h 

-  Mr. - f2/Z~2  —  -  fh 

U  {J2  +  H)  '  + 

X^  —  X4  I]  -j-  Ui  l\  ^  aon  j  Li  J i  ^  *on 

-t -  t\J\0.Qn  —  it  (R^Lx  —  p)  Mn, 


(4) 


where  H  =  H'/i22» 

Since  the  load  shaft  of  the  two-way  system  should  accurately  reproduce 
all  motions  of  the  operator  shaft,  the  desired  values  for  phase 
coordinates  xi,  X2,  and  X3  are  given  in  the  following  manner:  xi*  = 

=  ctOn(0>  x2 *  =  aon(t),  x3*  =  aon(0-  When  giving  the  desired  value  for  x4 ,  one 

should  proceed  from  the  fact  that  the  operator  should  sense  the  moment 
Qo ,  which  is  determined  only  by  the  object  of  manipulation.  Therefore, 

let  us  assume  that  x4*=nQ0=n[J0a,i(t) +h0],  where  Qo  =  Q  -  0^  Ql^  is  a 

component  of  moment  Q,  determined  by  the  parameters  of  the  links  of  the 

-v  /v 

actuating  member,  Jo  and  ho  are  the  components  of  elements  H'  and  h, 
respectively,  determined  by  the  parameters  of  the  manipulation  object, 
and  n  is  a  scaling  factor. 
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Since  q.  =  a„,  q.  =  then  h  =  h(tf„,  «„) .  Assuming  that  the  two-way 
Hi  H  H  rl 

system  operates  sufficiently  accurately,  i.e.,  variables  and 

hardly  differ  from  «on  and  aon,  respectively,  function  h(«H>  »H)  can  be 
expanded  into  a  Taylor  series  [3] 


h(a.H,  aH)=h(a0n  aon)  {aH— aon)  +^(aon.  a0ii)  (««— aon) +  C?  = 

-h(-2oa,  aon)e2  +  M<X0n,  a0n)ei+G, 


(5) 


where 


,  ,  •  v  _  dh 

(aon?  aon)  ^ 

H 


a  =ct  ,  a  —a 
H  on  H  on 


.  ,  *  X _  dh  I 

*  L  (aon»  aonJ  —  I  i 

H  ia  *=a  ,  a  = 

H  On  H  on 


G  =  ft(aHj  Oh) — h(ct0Tii  cc0n)^2  ^(Oom  Cton)fl* 


the  term  hn(ff„,  it,,)  can  also  be  expanded  into  components  similar  to  (5) 
OH  H 


^o(oh,  Oh)  =ho(a  on>  Cton)  62  “h^O  (cton>  C£on)£i“rGi 


0> 


(6) 


where 


^0  (®on»  ®on) 


dh* 


d* 


dh0 


H  la,,— a  ,  a  =  a 

H  on  H  on 


»  ^‘0  (aon»  aon) 


a m — ®  *  ®u — c 

ri  on  **  on 


Go  =  ho(o,H,  O.H ) — ho(aon,  a0n)82  A.(aon,  «on)ei. 


System  (4),  written  with  respect  to  errors  of  the  phase  coordinates, 
assumes  the  following  form  with  regard  to  the  type  of  desired  values  of 
the  phase  coordinates  and  expressions  (5)  and  (6): 


£1  —  01282, 

82  =  £2363* 

83  =  ^3181  +0.3282-}“  #3363  +  ^32^2  +/§  +  O3,  ,  j 

£4  =  o4,  e,  +  04262  +  04383  +  04484  +  641  o,  +  bi2u2 +f4  4-  co4>  i 
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where 


n  _  <7—1  „  _  n  ^M2^w2“t-  |  R2h  *r  (h  ~r  ?.)  ]  /ft 

U-12  —  u-23  —  1  *  aZi  71 - :  ,  ,  r  »  “  32  =  — - : — 


*33 


®-t{J*+H)'  MV"") 

I  R2  H~  hiq  \  f  R  \ 

_  —  |  )  J  fl« - *  [  4  a3i  +  77“  +  'Ao )  > 


/Z  Jq  CL-yt  ~j~  — 7 -  hn  ‘ 


dt! - 


R, 


b  —  '•„  I,  =  —  n 


K* 


J 0  |  a3»  i  / 


L  I  J_  /  _i_  A 
i  •'o  •  A£o 


r'i  a  ' N  m2  u 

4,44  “  ~  ’  &32  ~  ULt  Ul  ~  H)  ’  541  “ 


'l^Ml  I,  _  „  I  J, 

^  »  “42  ~~  ”*^0“32’ 


/a  §32^12*  /*4  S 42^t2»  §32  /,  +  /yj  i  <?42  —  ^0^32» 


CO3  -  ^3]  CXon  "f"  ^32^on  “I”  ^33 Cion  "h  /?Z31^  +  ^32^?, 


d/:>4  - 


^Ma^o>2 

L*x  (72  4-  //) 


’  ^32  ( is  y2  -i-  h  )  *  ^33—  1. 


m 
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W-* 


m„  = 


’  -33  .-«)  ’ 

C04  =  ^4 1  ®on  +  di2an„  +  dna,m -f  duMTl  -\-m4lG  +  in42G  4- m^G0 + muG0, 

fl  _  ,2  V»l  .  ,  ,  .2  /?,  ,  ,  ,'  ,  /?.  \  .  ] 

41  1  i.]  dJgd-ji,  =  il  J0\^d3i  T-  ~j~^)  "T  J(\  . 

di3=ii2Jlt  d4i  =  il  (Ri/Li+p) ,  nuA  =  —nJ0m3u  m42  =  nj0in32, 
mi3  =  —nR\iLu  m44  =  —n. 


It  is  obvious  that  system  (7)  is  nonstationary,  since  parameters  H,  Jo, 
h,  ho,  J,  Jo,  H,  Jo,  h,  ho,  i,  -lo  are  essentially  variables. 

Let  us  assume  that  the  rates  of  variation  of  these  parameters  are  much 
less  than  the  rates  of  the  transient  processes  in  the  system,  i.e.,  let 
us  assume  that  the  parameters  of  the  two-way  system  are  quasi-steady, 
and  the  coefficients  of  the  equations  of  system  (7)  then  remain  constant 
at  small  time  intervals.  In  this  case,  system  (7)  will  have  the 
following  form  in  vector  form  at  each  of  the  indicated  intervals  from 
the  entire  time  range  of  operation  of  the  two-way  system: 


E(t)=AE(i)+BU(t)+f(E)  +  W(q,  q,  q),  E(t0)  =  X (t0) —X* (t0) ,  (8) 

here  E(t)=X(t) _ X *(t),  £(()eS4  is  a  vector  of  the  errors  of  the  phase 

coordinates,  X(t)  £  R4  is  a  vector  of  the  phase  coordinates  of  the 
system,  X*(t)  £  R4  is  the  vector  of  the  desired  values,  U(t)  £  R2  is  a 
vector  of  the  control  actions,  B  €  R4*2  is  a  matrix  of  the  amplification 
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factors  during  control  actions,  A  €  R4*4  is  a  matrix  of  the  dynamic 
properties  of  the  two-way  system,  f(E)  6  R4  is  a  nonlinear  vector 

function,  and  W(q,  q,  q)  6  R4  is  a  vector  of  external  effects. 

Let  us  assume  that  matrices  A  and  B  vary  discretely  upon  transition  from 
one  time  interval  to  another. 

One  can  use  the  method  proposed  in  [4]  to  design  a  regulator  of  a 
two-way  system  on  any  of  the  indicated  time  intervals  in  view  of  the 
steady  nature  of  the  parameters  of  equation  ( 8 ) .  According  to  this 
method,  quasi-optimal  control  for  stationary  system  (8),  which  minimizes 
the  quadratic  criterion  of  quality,  is 


/= 0,5  J  [£T  (0  (DE  (t)  +  t/T  {t)  W(<)  ]  dt. 


t0  and  ti  are  the  initial  and  final  integration  time,  $  is  a  positive 
semidefinite  matrix,  and  $  is  a  positive  definite  matrix 


U (t)  = — D~  E (i) — DHf(E) — DUW, 

where  q-  =  4f-IBTK,  qh _ £>b _  — (AT — KB'F_1Bt) K. 
Matrix  K  is  determined  from  a  Riccati  equation 


KA  +  ATK— KB'F"‘BTK+ O  =  0. 


In  scalar  form,  equation  (9)  can  be  rewritten  in  the  form 


W;  =  Kll  6a  +  KnO-H  +  Ki3O.au  +  K14O.H  +  Klaxon  +  KieSv  +  (  10  ) 

+  Knh  +  Kigh  +  Ki9ho+ Knoho+ Kin  sgn  oni  +  Km  sgn  aon  4-  A’mccon, 


where 


8a - El  —  (Xon — OlH',  8N 


84  =  M0n— llQo, 


Kn  =  Dn  4-  Dii +  m:vK)  4-  Du  (w4)X  4-  mvj.  4-  tfz43X()  4- 
Kn  —  —  Dn  +  D%  \mMh  4-  (h  4-  /•)]  4-  D%  j muh  4-  rni2  ( fi  —■  1.)  4- 

~r  mi3^o  “f-  ^44  (^"0  4”  >,)].  Kn  —  Ki 2  —  Dnd3x  —  Dndu, 

Ku  =  -  Dn  4-  Dntn^h  +  D^m^h  4-  mA4fi0),  Kn  =  -  Kn  -  D%d3!-  Oft*42, 

K is  —  D14,  Kp  —  Di3ffi3i  Duffle,  K =  —  Dnm33  — 
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K[ g  —  —  Dtjm.l3,  K  no  —  —  Duffle,  Km  —  (  Dizg^i  *  Kl-a, 

Ku2  =  —  DuduMrit  K ii3  =  —  Dud33  —  D^d^, 

M'tj  and  are  the  steady  frictional  moments  of  motion,  and  1=1,  2. 

When  realizing  equation  (10),  the  terms  containing  differentiation  in 
gg 2 ,  g4  2>  and  d4  4  can  be  disregarded,  since  they  have  no  influence  on 
the  dynamic  properties  of  the  two-way  system  after  the  beginning  of 

motion  of  the  system  [4].  Moreover,  terms  proportional  to  a  ^  must  be 

disregarded,  since  the  value  of  a  Qn  can  not  be  measured  for  use  in 
control . 

Having  designed  a  stationary  regulator  for  each  of  the  considered  time 
intervals  of  operation  of  the  two-way  system,  knowing  the  corresponding 
intervals  of  the  values  of  the  variable  parameters,  one  can  calculate 

the  dependence  of  coefficient  Kim  (1  =  1,  2,  m  =  1,  13)  of  the  regulator 
of  a  two-way  system  on  these  variable  parameters.  The  functions 
described  in  this  manner  are  written  to  the  computer  memory  and  are  then 
used  to  organize  adjustment  of  coefficient  Kim  of  the  regulator  of  the 
self-adjusting  two-way  system  as  a  function  of  the  current  values  of  the 
indicated  variable  parameters  of  the  manipulator. 

However,  it  is  rather  complicated  to  design  and  then  to  realize  a 
self-adjusting  two-way  system  with  a  large  number  of  variable 
parameters,  by  which  Kim  should  be  adjusted.  All  the  enumerated 
parameters  are  functions  of  the  generalized  coordinates  of  the  actuating 

member  qj ,  of  their  derivatives  qj ,  qj  ( j  is  the  number  of  the  degree  of 

freedom,  j  =  1 ,  N,  j  ^  i)  and  of  the  mass  of  the  object  of  manipulation 
mo.  Therefore,  the  dependence  of  coefficients  Kim  of  the  synthesized 
regulator  of  a  two-way  system  on  the  earlier  indicated  12  variables  can 

be  transformed  as  a  function  of  coordinates  qj ,  qj ,  qj  and  mo .  If  only 
three  transferable  degrees  of  freedom,  which  mainly  determine  variation 
of  the  load  parameters  of  the  two-way  system  [5],  are  taken  into 
account,  the  values  of  Kira  will  be  dependent  on  a  maximum  of  7  variable 

parameters  rather  than  on  12,  since  j  =  1,  3,  j  ^  i.  Since  the  mass  mo 
during  operation  of  a  copying  manipulator  ordinarily  remains  fixed, 
after  calculation  of  it,  for  example,  by  the  relations  found  in  [6]  on 
the  initial  step  of  motion  of  the  manipulator,  continuous  adjustment  of 

the  amplification  factors  Kim  to  functions  of  only  six  variables  qj ,  qj , 

and  qj  is  required. 
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Studies  of  the  effectiveness  of  the  proposed  approach  were  conducted  on 
the  example  of  a  two-way  system  of  the  first  degree  of  freedom  of  a 
copying  manipulator,  the  kinematic  layout  of  the  actuating  member  of 
which  is  presented  in  Figure  1 .  The  considered  actuating  member  has  the 
following  parameters:  lengths  of  links  l!  =  12  =  0.5  m,  mass  of  links 
mi  =  10  kg  and  m2  =  5  kg,  moments  of  inertia  of  the  links  with  respect 
to  the  axes  of  the  joints  J|  =0.85  kg*m2,  J2  =  0.42  kg*m2,  and  mo  = 

=  1  kg.  The  numerical  values  of  the  parameters  of  the  fixed  part  of  the 
two-way  system  are  presented  in  the  table. 
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The  mode  of  stimulated  motion  of  the  operator  shaft  of  the  first  degree 

of  freedom  according  to  the  law  aon=  1  (0/(0,167p+ l)3  at  different  modes 

of  motion  of  the  second  degree  of  freedom  of  the  copying  manipulator  was 
considered  during  studies  of  a  two-way  system.  Curves,  corresponding  to 
variation  of  the  angle  of  rotation  of  the  second  degree  of  freedom  by 

the  law  a/f2==  1 ,5  [  1  +  sin  (it/ + n/2)  ]  ,  are  presented  in  Figure  2,  while 

curves  corresponding  to  variation  of  the  angle  of  rotation  by  the  law 

ctH2=  1.5H  +  sin(jii — n;/2)] .  are  presented  in  Figure  3.  The  number  1  denotes 

curves  for  the  case  when  a  self-adjusting  two-way  system  is  accomplished 
only  by  the  variable  inertial  parameters  similar  to  [ 1] .  The 
centrifugal  and  Coriolis  forces  during  shaping  of  the  control  signal 
(control  law)  were  disregarded.  The  curves  denoted  by  the  number  2 
correspond  to  the  case  when  the  parameters  of  the  regulator  of  the 
two-way  system  are  adjusted  as  a  function  of  variation  of  all  the  load 
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parameters  (with  regard  to  the  moments  of  high-speed  forces)  according 
to  the  approach  considered  in  this  paper. 


Figure  2 


Analysis  of  the  graphs  presented  in  Figures  2  and  3  permits  one  to 
conclude  that  the  method  of  synthesis  and  design  of  a  self-adjusting 
two-way  system,  proposed  in  this  paper  and  that  takes  into  account  the 
total  dynamics  of  the  actuating  member  of  the  copying  manipulator, 
provides  a  significant  improvement  of  the  quality  of  operation  of  the 
system  with  significant  and  rapid  variation  of  all  the  load  paramters. 
However,  the  synthesized  regulators  are  rather  complex  to  realize  and 
practice.  Therefore,  one  can  attempt  to  simplify  self-adjusting 
regulators  without  significant  deterioration  of  the  quality  of  the  main 
operating  modes  of  two-way  systems. 


The  results  of  the  studies  showed  that  parameters  h  and  ho  can  be 
disregarded  during  formation  of  control  law  (10)  at  speeds  of  motion  of 
the  working  member  of  the  manipulator  up  to  5  m*s.  Moreover, 

«v 

coefficients  Kim  can  be  adjusted  only  by  the  values  of  H  and  Jo,  but  h 

and  ho  must  be  included  in  control  law  (10).  The  errors  of  the 
synthesized  two-way  system  are  denoted  by  the  number  3  upon  introduction 
of  the  indicated  simplifications  in  Figures  2  and  3.  It  is  obvious  from 
these  figures  that  curves  3  hardly  differ  from  curves  2  for  the 
considered  operating  modes,  which  correspond  to  a  two-way  system  with 
accurate  and  complete  self -adjustment  for  all  variable  load  parameters. 
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Figure  3 

Of  course,  introduction  of  the  indicated  simplifications  does  not 
guarantee  high-quality  operation  of  two-way  systems  in  all  operating 
modes  of  copying  manipulators;  therefore,  the  efficiency  must  be  checked 
at  all  possible  motions  of  the  actuating  member  when  using  the 
simplified  regulator.  However,  studies  and  analysis  of  curves  2  and  3 
permit  one  to  hope  that  the  coefficients  of  the  regulator  can  sometimes 
be  adjusted  by  only  two  parameters  H  and  Jo .  This  considerably 
simplifies  practical  realization  of  self-adjusting  two-way  systems 
without  significant  deterioration  of  the  indicators  of  its  operating 
quality. 
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[Text]  An  adaptive  robotic  complex  with  robotic 
vision  and  force— moment  sensitization  systems  has  been 
developed  for  theoretical  and  experimental  studies  of 
the  problems  of  using  sensory  information  in  control 
of  industrial  robots. 

Modern  experience  of  practical  use  of  industrial  robots  (PR)  showed  the 
need  to  develop  robotic  complexes  (RTK) ,  oriented  toward  fulfillment  of 
specific  manufacturing  operations  and  built  into  the  hierarchical 
structure  of  the  management  system  of  a  flexible  automated  plant. 

Special  sttention  is  devoted  to  giving  these  complexes  adaptive 
properties  to  reduce  expenditures  for  organization  of  the  working 
environment  and  to  readjust  the  manufacturing  equipment. 

A  number  of  problems  such  as  interaction  of  the  robotic  complex  with  the 
upper  level  control  system,  synchronization  of  simultaneous  operation  of 
several  robots,  processing  and  utilization  of  various  sensory 
information  for  control  of  robots,  depiction  of  the  operating 
environment  in  the  form  of  a  model  of  the  external  world,  planning  the 
motions  of  robots  with  regard  to  variations  of  this  model,  teaching 
robots  and  so  on  must  be  solved  in  development  of  adaptive  robotic 
complexes.  An  experimental  adaptive  robotic  complex  was  developed  to 
conduct  theoretical  and  experimental  studies  in  this  field.  The  complex 
is  oriented  toward  completion  of  mechanical  assembly  operations , 
machining,  painting  and  soldering.  The  work  is  being  conducted  in 
several  directions.  The  basic  direction  is  related  to  development  of 
the  algorithms  and  software  of  robots,  sensory  systems  and  the  control 
computer  of  the  robotic  complex.  Much  attention  is  devoted  to  design  of 
multicomponent  force-moment  sensors  and  of  passive  centering  devices  for 
assembly  of  mechanical  connections. 
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The  developed  robotic  complex  includes:  a  control  computer,  two 
industrial  robots  of  type  KM-01  (PUMA  560),  television  and  photoarray 
robotic  vision  systems,  two  force-moment  sensitization  systems,  a  number 
of  passive  centering  devices  and  a  lighting  system.  When  studying 
individual  manufacturing  operations,  the  robotic  complex  is  delivered 
with  additional  fittings  and  tools. 

The  central  part  of  the  developed  robotic  complex  is  an  Elektronika  15 
VUMS-28  computer.  Connection  between  separate  subsystems  of  the  complex 
is  provided  by  using  a  serial  interface  of  the  Styk  S2  type  and  a 
parallel  interface  of  the  IEPR  type  (see  figure).  The  software  of  the 
robotic  complex  supports  the  following  main  functions:  coordination  of 
the  joint  operation  of  two  manipulators,  synchronization  of  the 
operation  of  individual  subsystems,  compilation  of  industrial  robot 
control  programs  in  ARPS  language,  and  loading  of  them  into  the 
industrial  robot  controller  (UU)  Sfera-36,  and  transfer  of  data  files 
and  control  commands  of  subsystems. 


Block  Diagram  of  Experimental  Adaptive  Robotic  Complex: 

UU — Sfera-36  controller  of  RM-01  industrial  robot; 

TV — television  camera;  K  STZ — robotic  vision  system 
controller;  FM — photomaster  camera;  K  EM — photomaster 
controller;  SMD1 — force-moment  sensor  in  base  of  assembly 
table;  SMD2 — force-moment  sensor  in  wrist  of  industrial 
robot;  KSMD — force-moment  sensor  controller;  MK  SSO — 
force-moment  sensitization  system  microprocessor  controller; 
MTsU — mechanical  centering  device;  S — industrial  robot  grab; 
SO — lighting  system;  TO — manufacturing  equipment;  EVM  VU — 
upper-level  control  computer 

(Key  on  following  page) 


52 


KEY: 


1 .  Styk  S2  interface 

2.  Styk  S2  interface 

3.  Robotic  vision  system 
computer 

4.  Robotic  complex  computer 

5 .  Upper-level  control 
computer 

6 .  Television  camera 

7.  Robotic  vision  system 
controller 

8 .  IRPR  interface 

9 .  Force-moment  sensor 
controller 

10.  Force-moment  sensor 


1 1 .  Lighting  system 

12.  IRPR  interface 

13.  Force-moment  sensitization 
system  microprocessor 
controller 

14.  Sfera-36  controller 

15.  Manufacturing  equipment 

16 .  Controller 

17.  Photop)attern  controller 

18.  Industrial  robot  grab 

19.  PUMA  560  industrial  robot 

20 .  Photopattern 

21.  Mechanical  centering  device 


The  complex  is  controlled  through  special  commands,  written  in  a 
separate  file  or  introduced  by  the  operator  from  the  keyboard  of  the 
robotic  complex  computer  in  the  form  of  symbolic  lines.  The  commands 
are  executed  in  the  interpretation  mode.  Studies  are  now  under  way  to 
develop  a  robotic  complex  control  language. 

Robotic  vision  systems  (STZ).  Two  robotic  vision  systems  are  used  in 
the  developed  robotic  complex:  a  television  and  photopattern  system. 

The  television  robitic  vision  system  is  designed  to  identify  parts,  to 
determine  their  position  and  orientation,  and  also  for  visual  inspection 
of  products.  It  includes  a  fixed  television  camera  with  standard 
television  signal,  an  Elektronika  NMS-1201.02  microcomputer,  and  an 
interface  module  for  the  computer  to  communicate  with  the  camera.  The 
interface  module  enters  a  frame  of  the  image  measuring  256  *  256 
elements  and  writes  it  to  the  computer  RAM  in  a  format  of  1  bit/element. 
The  image  is  converted  to  the  binary  system  according  to  a  program- 
assigned  threshold  having  64  levels. 

The  photopattern  camera  is  mounted  on  the  industrial  robot  grab  and  is 
used  to  organize  the  external  feedback  according  to  the  position  of  the 
grab  with  respect  to  the  object.  The  pointing  function  is  realized  for 
a  given  moment.  Besides  a  camera,  the  photopattern  robotic  vision 
system  includes  a  photopattern  controller  with  integration  module  with 
the  Q-bus  of  the  Sfera-36  controller.  An  image  frame  measuring  32  x  32 
elements  with  64  levels  of  brightness  is  written  to  the  RAM  of  the 
controller . 

The  robotic  vision  system  is  programmed  to  perform  a  specific  operation 
using  tool  software  (PO) .  The  latter  consists  of  universal  and 
specialized  devices  and  functions  under  control  of  the  RAFOS  operating 
system . 

The  specialized  part,  which  is  the  core  of  the  image  processing  software 
in  the  robotic  vision  system  of  an  assembly  robotic  complex,  includes 
image  analysis  and  synthesis  systems,  functioning  modes  of  the  robotic 
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vision  system  and  communication  modes.  The  basis  of  each  system  is  a 
set  of  libraries  that  guarantee  development  and  study  of  visual  analysis 
procedures,  and  also  creation  of  applications  software.  The  analysis 
system  includes  image  shaping  libraries  (formed  by  image  input  and 
output  modules),  description  libraries  (modules  for  construction  of 
models,  initialization  of  data  structures,  analysis  of  connectedness, 
description  of  components,  formation  of  a  list  and  image  of  the  model) 
and  image  interpretation.  The  image  synthesis  library  includes  input, 
primitive,  filling  and  image  modules;  the  communication  libraries 
including  loading,  dialogue  and  systems  subroutine  modules;  and  the 
teaching  libraries  include  mode,  statistics  and  image  modules. 

The  program  modules  of  the  libraries  process  the  images,  shaped  by  both 
television  and  photopattern  cameras.  The  operating  mode  of  the  modules 
is  determined  by  the  values  of  the  parametes  assigned  during  generation 
of  the  robotic  vision  system  software. 

To  implement  the  functions  of  the  robotic  vision  system,  the  image  to  be 
analyzed  is  represented  in  the  form  of  an  oriented  root  tree  of 
components ,  i . e . ,  of  image  domains ,  corresponding  to  individual  parts  or 
fragments  of  them.  Each  vertex  of  the  tree  representing  a  component  is 
characterized  by  a  set  of  parameters  that  determine  the  absolute  and 
relative  position  and  orientation  of  the  component  and  that  guarantees 
identification  of  it.  The  "encompassing-embedded  component"  ratios  are 
given  between  the  vertices  of  the  tree.  The  image  is  described  within 
some  window  of  rectangular  or  arbitrary  shape,  and  the  position  and 
dimensions  of  the  window  are  assigned  earlier  when  teaching  the  system 
or  are  determined  directly  during  the  manufacturing  process. 

The  developed  a  priori  glossary  of  component  parameters  includes  those 
calculated  on  the  set  of  component  elements,  contour  elements  and  local 
characteristics,  and  also  a  number  of  relative  values,  and  permits  one 
to  describe  objects  of  different  shape.  The  developed  algorithms  permit 
one  to  use  various  types  of  representation  of  contours  and  is  adapted  to 
the  shape  of  the  objects,  which  results  in  an  increase  of  the  speed  of 
the  system.  The  algorithms  for  analysis  of  the  connectedness  and 
description  of  components,  used  in  the  system,  guarantee  processing  of 
images  of  scenes  of  the  assembly  processes  in  real  time  (1-3  s). 

The  main  restrictions  imposed  on  the  composition  of  the  software  of 
industrial  robotic  vision  systems  are  the  low  main  memory  capacity  of 
the  computer  and  the  need  to  solve  problems  in  real  time  with  regard  to 
the  characteristics  of  the  manufacturing  process. 

To  automate  the  process  of  developing  the  robotic  vision  system 
software,  it  is  suggested  that  a  hierarchical  system  be  used,  the  task 
of  which  includes  a  flexible  interface  with  operator,  information 
support  of  the  decision  search  process,  evaluation  of  the  software  as  a 
whole  and  evaluation  of  individual  modules  according  to  the  used 
criteria,  storage  and  generalization  of  the  operator’s  action  and  a 
number  of  auxiliary  functions.  Each  program  module  is  evaluated  by  the 
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quantitative  and  qualitative  parameters  such  as  processing  time, 
capacity  of  the  memory  occupied,  degree  of  universality,  quality  of 
processing  and  so  on.  It  is  suggested  that  expert  analyses  be  used  with 
regard  to  the  possibility  of  total  formalization  of  parameters. 

Linguistic  variables  (LP),  designed  on  the  basis  of  odd  set  theory,  are 
considered  for  design  of  a  system  of  expert  evaluations.  A  context-free 
grammar  of  linguistic  variables  is  worked  out,  which  is  a  group  of  four 


G  =  < A\  TO' V),  U,  S), 


where  N  is  the  number  of  the  linguistic  variable,  T(N)  is  the  term  or 
term-set  of  the  values  of  the  linguistic  variable,  U  is  the  set  of 
values  of  the  basic  variable,  and  S  is  the  syntactical  rule  for 
formation  of  the  linguistic  variable. 

The  characteristic  feature  of  the  introduced  grammar  includes  the 
universality  of  the  image  of  both  even  and  interval  values  and  of 
qualitative  estimates. 

The  criterion  of  the  usefulness  of  the  module  within  the  robotic  vision 
system  software  is  an  unevenly  suspended  discriminate  function 


NEC=  Z  to, 


where  tfi  is  the  coefficient  of  preference  of  the  i-th  parameter  and  ei 
is  expert  analysis  of  the  i-th  parameter.  As  a  result,  the 
corresponding  values  of  the  preference  factors  were  found  for  different 
types  of  robotic  vision  systems.  Variants  of  the  configuration  of  the 
robotic  vision  system  software  were  selected  and  they  were  generated 
after  ranking  of  the  available  program  modules  by  the  value  of  the 
utility  factor. 

Force-moment  sensitization  systems  (SSO) .  The  experimental  robotic 
complex  includes  two  force-moment  sensitization  systems.  The  sensor  of 
the  first  system  is  located  in  the  base  of  the  assembly  table  and  is 
used  to  measure  the  components  of  the  main  force  and  moment  vector, 
occurring  upon  assembly  of  the  unit,  mounted  on  the  table.  Force-moment 
information  is  introduced  and  processed  in  the  robotic  complex  computer 
and  is  used  to  check  the  assembly  process  and  to  control  the  robots. 

the  second  system  is  made  in  the  form  of  a  separate  design,  located 
between  the  last  link  of  the  industrial  robot  and  its  grab.  The 
controller  built  into  the  housing  of  the  sensor  guarantees  input, 
processing  and  transmission  of  the  force-moment  information  in  the 
Sfera-36  controller  for  organization  of  external  force  feedback.  The 
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1816  VE48  single-chip  microcomputer  with  ROM  having  capacity  of  4  kbytes 
is  included  in  the  controller,  which  permits  implementation  of 
specialized  data  processing  programs.  The  indicated  programs  are 
designed  to  compensate  for  the  temperature-time  instability  of  the 
measuring  circuits,  to  take  into  account  the  static  components  of  forces 
and  force  moments,  caused  by  the  weight  of  the  grab  and  tool,  and  also 
to  present  a  system  of  forces  from  the  center  of  measurement  of  the 
sensor  to  a  point,  connected  to  the  end  executive  of  the  manipulator. 
Location  of  the  controller  inside  the  sensor  housing  permits  one  to  use 
the  concept  "intelligent  sensor,"  according  to  which  force-moment 
information  is  processed  autonomously.  This  frees  to  a  considerable 
degree  the  computer  resources  of  the  control  computer  of  the  robot. 

The  force-moment  sensors  of  both  systems  contain  strain  gage  resistor 
sensitive  elements  and  have  diagonal  stiffness  matrices,  which  permits  a 
level  of  cross-connections  between  channels  in  the  range  of  3-5  percent 
at  speed  of  the  force-moment  sensitization  system  of  not  worse  than 
100  /is. 

Mechanical  centering  devices  (MTsU) .  Along  with  force-moment 
sensitization  devices,  the  experimental  robotic  complex  includes  passive 
force-moment  adaptation  devices,  realized  in  the  form  of  mechanical 
centering  devices  with  remote  stiffness  center.  The  device  is  located 
in  the  wrists  of  the  industrial  robot  and  is  designed  to  compensate  for 
the  mutual  positioning  errors  of  conjugate  parts  of  the  "shaft-bushing" 
type  due  to  the  effect  of  forces  occurring  during  assembly.  The 
presence  of  a  remote  stiffness  center  reduces  the  probability  of 
jamming.  The  mechanical  centering  device  is  designed  in  the  form  of  a 
combination  of  linear  and  angular  pliancy  modules,  which  permits  one  to 
provide  small  values  of  linear  and  angular  stiffnesses  simultaneously 
(up  to  1,000  N/m  and  6  N* m/rad),  and,  accordingly,  to  provide  forces  of 
interaction  of  parts  at  a  significant  distance  of  the  stiffness  center 
of  the  device  (up  to  200  mm)  from  its  end  flange.  The  high  reliability 
of  the  device  is  achieved  by  using  elastomeric  structures  of  directed 
pliancy  as  the  elastic  elements,  which  are  subjected  to  shearing  and 
compression,  and  which  do  not  experience  tensile  loads. 

The  utilized  mechanical  centering  devices  contain  interchangeable 
modules,  having  different  stiffness  and  geometric  characteristics,  which 
permits  online  adjustment  of  the  robotic  complex  for  assembly  of  other 
types  of  products. 

COPYRIGHT:  "IZVESTIYA  VUZOV  SSSR  -  PRIBOROSTROYENIYE" ,  1990 
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[Text]  The  operating  characteristics  of  the  force 
sensor  of  an  industrial  robot  in  a  technological 
cleaning  operation  are  considered.  A  mathematical 
model  of  the  force  sensor  is  presented,  and  the 
operating  and  working  principles  of  individual 
functional  modules  of  the  force  sensitization  system 
are  described. 

Such  technological  operations  as  rough  grinding,  abrasive  and  milling 
cleaning  can  be  combined  into  one  group  of  operations  of  the  "trimming" 
type.  Regardless  of  the  type  of  material  to  be  machined,  the 
characteristic  features  of  this  group  of  operations  are  the  presence  of 
random  projection  and  use  of  a  rotating  milling  tool,  which  determines 
the  dynamics  of  interaction  with  the  part  to  be  machined.  An  abrasive 
disk  or  cutter  is  frequently  used  as  the  tool.  Practical  solutions  are 
implemented  slowly  with  the  obvious  importance  of  automating  these 
human-hazardous  technological  operations,  which  is  explained  by  the 
complexity  of  the  technical  task.  We  in  the  USSR  already  have  several 
industrial  robotic  trimming  robotic  complexes,  while  industrial  robots 
are  used  extensively  in  these  operations  abroad,  judging  by  the 
publications  and  recent  exhibitions.  Power  sensors  are  not  yet  used  in 
trimming  robotic  complexes  in  the  USSR.  For  example,  thermosetting 
plastic  parts  are  trimmed  at  the  Elektrosila  Association  without  a  power 
adaptation;  the  trimming  complexes  manufactured  by  NPO  [scientific 
production  association]  Granat  are  also  not  yet  equipped  with  force 
sensors.  The  use  of  force  feedback  permits  one  to  develop  a  trimming 
robot  control  system,  having  the  highest  machining  performance,  since 
the  extent  of  the  projection  is  indirectly  estimated  through  measurement 
of  the  cutting  forces,  which  permits  one  to  regulate  the  rate  of 
delivery  of  the  tool  along  the  contour  of  the  part  as  a  function  of  the 
real  value  of  the  projection.  The  use  of  adaptation  according  to  the 
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cutting  force  also  increases  the  safety  of  trimming  and  the  quality  of 
machining*  A  large  number  of  designs  of  various  force  sensors  for 
industrial  robots  are  described  in  the  Soviet  and  foreign  literature  [1, 
2],  but  there  are  no  descriptions  of  force  sensitization  systems  (SSO) 
for  trimming  robots. 

A  force  sensitization  system  has  been  developed  and  manufactured  at  MGTU 
[Moscow  State  Technical  University]  imeni  N.  E.  Bauman,  specially 
designed  for  trimming  operations.  The  sensitization  system  has  the 
following  specifications: 


Number  of  measured  components  (total  force 

vector  at  any  point  of  work  zone  of  sensor)  3 

Range  of  measured  forces  (for  each  component),  N  ±50 

Accuracy  by  channels,  pecent  3,  5,  5 

Maximum  frequency  of  input  signal,  kHz  50 


The  specifics  of  the  given  force  sensitization  system  includes  the  fact 
that  the  total  cutting  force  vector  in  the  region  of  contact  of  the  tool 
and  part  can  be  measured  at  any  point  of  the  circumference  of  the  tool. 
Since  the  force  vector  is  measured  at  a  specific  spatial  region  (work 
zone)  of  the  force  sensor  coordinate  system  rather  than  at  a  fixed 
point,  no  restrictions  are  imposed  on  the  diameter  of  the  tool,  and  the 
given  force  sensitization  system  can  be  used  not  only  with  small- 
diameter  end  mill  cutters,  but  also  with  disk  cutters,  and  with  abrasive 
and  diamond  disks.  During  trimming,  the  contact  zone  of  the  tool  and 
part  can  be  shifted  around  the  circumference  of  the  tool  and  can 
significantly  change  its  coordinates,  which  will  not  essentially  affect 
the  measurement  of  the  force  vector  components. 

A  block  diagram  of  the  force  sensitization  system  of  a  trimming  robot  is 
shown  in  the  figure.  According  to  this  layout,  several  force 
sensitization  systems  have  been  manufactured  at  MGTU  imeni  N.  E.  Bauman 
for  different  types  of  robots.  Let  us  consider  in  more  detail  a  version 
of  a  force  sensitization  system  for  the  RM-01  industrial  robot  with 
FUMA-NOKIA  electromechanical  actuating  device.  The  force  sensitization 
system  consists  of  three  functionally  complete  modules:  a  force  sensor, 
a  force  data  input  card  and  a  program  module,  which  additionally 
processes  the  force  signal. 

The  force  sensor  is  a  structurally  complete  product,  included  in  the  gap 
of  the  force  circuit  between  the  mounting  flange  of  the  robot  and  the 
working  tool  (a  pneumatic  turbine  with  cutter).  The  force  sensor 

-t 

measures  the  cutting  force  vector  F,  which  occurs  in  the  contact  region 
of  the  tool  and  part,  and  emits  in  analog  form  three  components  of  the 
cutting  force  in  the  coordinate  system  of  the  sensor.  KF-5P  foil  strain 
gages  with  base  of  2  mm  and  with  rated  resistance  of  200  ohms  are  used 
as  the  sensitive  elements.  The  principle  of  structural  separation  of 
force  factors,  similar  to  that  presented  in  [3],  is  used  in  the 
considered  sensor.  But  the  force  sensitization  sensor,  designed  on  the 
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basis  of  the  sensor  developed  at  MGTU  imeni  N.  E.  Bauman,  also  has  new 
properties,  manifested  in  measurement  of  forces  exactly  on  trimming 
operations . 


Block  Diagram  of  Robot  Force  Sensitization  System 


KEY : 
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Force  sensor 

9. 
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2. 

3. 
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The  primary  signal  vector  of  the  sensitive  elements  R,  the  dimension  of 
which  is  determined  by  the  selected  number  of  measuring  points,  is 
introduced  to  analyze  the  conversion  of  data  in  the  force  sensor*  Each 
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component  of  the  vector  R  is  an  electric  signal,  taken  from  the 
measuring  circuit  of  strain  gages,  glued  to  a  specific  region  of  the 
flexible  element*  The  strain  gages  are  ordinarily  connected  by  the 
measuring  bridge  or  semibridge  layout,  which  affects  only  the  level  of 
the  removed  signal  and  does  not  affect  the  functional  conversion  of  the 
signal;  therefore,  let  us  assume  that  a  bridge  or  semibridge  is  one 
generalized  sensitive  element,  which  measures  deformation  at  one  point 
of  the  surface  layer  of  the  flexible  element.  Let  us  also  introduce  the 
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relative  deformation  vector  E,  the  components  of  which  are  the  relative 
deformations  of  the  surface  layers  of  flexible  elements  at  the  measuring 
point.  With  regard  to  the  assumption  that  one  generalized  sensitive 
element  measures  deformation  at  one  point  of  the  surface  layer  of  the 
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flexible  element,  let  us  take  the  dimension  of  vector  E,  equal  to  that 
— +  — +  — * 

of  vector  R.  Vectors  E  and  R  can  then  be  linked  by  the  matrix  of  the 

-+  -4 

transfer  coefficients  of  sensitive  elements  B  by  the  relation  R  =  BE, 
where  matrix  B  is  a  diagonal  B  =  diag[bii,  b2  2*  b33,  b4  4,  b5  5,  b66l,  and 
its  elements  are  constants,  determined  by  the  coefficients  of  the 
strain-sensitive  strain  gages,  by  their  connection  circuits,  and  also  by 
the  power  supply  voltage  of  the  bridges.  The  sensor  generally  measures 

-+ 

the  total  vector  of  force  factors  G,  consisting  of  the  principal  vector 

-4  -4 

F  and  of  principle  moment  M,  reduced  to  the  measuring  center  of  the 
sensor : 


a 


If  the  matrix  of  sensitivity  of  the  flexible  element  of  sensor  A  is 
introduced,  the  vector  of  relative  deformation  E  will  be  determined  as 

-4  -4 

E  =  A*G. 

One  can  assume  for  most  trimming  operations  that  the  contact  region  of 
the  tool  and  part  is  small  and  can  be  considered  as  a  point.  The  main 
principle  of  the  force  sensor  is  that  its  measuring  center  is  mobile  and 
coincides  with  the  contact  point  of  the  tool  and  part.  The  main  moment 
of  external  forces  is  equal  to  zero  with  respect  to  this  point  and  the 

total  vector  of  the  external  force  factors  G,  measured  by  the  sensor, 

will  consist  only  of  the  principal  vector  F,  i.e.,  its  dimensionality 
will  be  equal  to  three.  The  dimensionality  of  the  primary  signal  vector 
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of  sensitive  elements  R  and  accordingly  of  the  relative  deformation 
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vector  E  in  the  sensor  is  equal  to  six,  i.e.,  to  the  number  of  selected 
measuring  points.  This  is  determined  by  the  need  to  measure,  besides 
the  three  components  of  the  force  vector,  three  additional  coordinates 
of  the  movable  contact  point  of  the  tool  and  part.  The  sensitivity 
matrix  of  the  flexible  element  of  the  sensor  A  is  thus  rectangular  with 
dimensionality  6x3.  the  primary  signals  of  the  sensitive  elements  are 

-4 

amplified  by  the  pre-amplifiers,  i.e.,  vector  R  is  transformed  to  signal 

-4  -4 

vector  S  =  CR,  where  C  is  the  diagonal  matrix  of  the  force  factors  of 
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pre-amplifiers,  and  C  =  diag[c| i ,  C22>  C33,  C44,  C55,  egg].  After  the 
six-dimensional  vector  S  has  been  processed  by  the  computer,  a  three- 

-»  -4 

dimensional  vector  of  the  measured  force  U  =  DS  is  found.  Here  D  is  a 
rectangular  3*6  matrix  of  the  computer.  Complete  transformation  of 

-t  -♦ 

the  data  by  the  force  sensor  can  be  expressed  by  the  relation  U  =  HF,  in 
which  H  =  DCBA. 

The  considered  sensor  has  a  sensitivity  matrix  of  the  flexible  element  A 


A 


~an  0  0 

an  0  0 

0  aSi  0 
0  ait  0 
0  0  a58 

_0  0  &63_ 


and  computer  matrix  D 


(1) 


D  = 


0 

0 


^12  0 
0  dK 
0  0 


0  (/  0  ’ 
du  0  0 

0  d 


(2) 


Expressions  (1)  and  (2)  show  an  ideal  type  of  matrices  A  and  D,  since 
realistically  all  the  elements  of  the  indicated  matrices  may  not  be  zero 
matrices,  which  results  in  a  mutual  effect  between  the  channels  and  an 
increase  of  the  total  measurement  error.  As  the  test  results  of  the 
sensor  showed,  the  measurement  error  with  regard  to  the  reciprocal 
influence  of  the  channels  does  not  exceed  5  percent  in  the  temperature 
range  of  15-30° C — this  is  sufficient  to  design  an  adaptive  force  control 
system  of  a  trimming  robot.  The  flexible  element  of  the  sensor  is  made 
of  40Cr  steel,  subjected  to  heat  treatment,  and  consists  of  three 
assembly  units,  pinned  during  assembly.  The  flexible  element  can  also 
be  made  a  unit  element,  but  this  complicates  the  manufacture  and 
subsequent  gluing  of  the  strain  gages.  The  overall  dimensions  of  the 
sensor  (60  x  65  x  80)  essentially  do  not  limit  the  service  angles  of  the 
actuating  device.  The  developed  design  is  technologically  effective, 
but  imposes  high  requirements  on  the  accuracy  of  the  production 
equipment  and  can  be  manufactured  at  any  plant.  The  electronic  part  of 
the  sensor  consists  of  an  additional  internal  stabilizer, 
pre-amplifiers,  and  a  computer,  realized  on  153UD5  operational 
amplifiers . 
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The  force  data  input  card  is  made  in  the  design  of  the  Sfera-36 
controller  and  is  a  specialized  high-speed  multichannel  analog-to- 
digital  converter.  The  input  card  is  connected  to  the  force  sensor  by  a 
cable,  through  which  power  of  ±15  V  is  fed  to  the  sensor  and  three 
analog  signals,  corresponding  to  the  components  of  the  measured  cutting 
-» 

force  vector  F,  are  received.  A  shown  in  the  figure,  the  basis  of  the 
force  data  input  card  is  an  analog-to-digital  converter,  realized  on  a 
K1107PV1  integrated  analog-to-digital  converter.  The  high  speed  of  the 
analog-to-digital  converter  and  of  the  K590KN6  multiplexer  guarantees 
switching  of  the  required  channel  and  subsequent  analog-to-digital 
conversion  within  time  less  than  10  /is,  i.e.,  less  than  the  maximum 
permissible  delay  of  the  "Response"  signal  of  the  passive  device  on  the 
channel  of  the  central  processor.  This  permits  one  to  do  away  with  an 
internal  data  buffer,  and  the  central  processor  immediately  reads  the 
code  of  the  converted  value  from  the  analog-to-digital  converter.  The 
input  card  in  the  address  space  of  the  central  processor  has  three 
addresses,  and  each  address  corresponds  to  its  own  channel  (the 
component  of  the  measured  force  vector).  After  decoding  the  address, 
the  multiplexer  switches  the  channel  corresponding  to  this  address  and 
starts  the  analog-to-digital  converter,  which  emits  the  received  code  to 
the  data  address  bus  through  the  channel  receivers  after  completion  of 
conversion. 

Final  processing  of  the  force  signal  occurs  in  the  program  module  for 
additional  processing  of  the  force  signal  (see  figure),  which  solves  the 
following  tasks:  it  filters  the  low-frequency  noise  component,  switches 
the  force  signal,  compensates  for  the  gravitational  component  of  the 
force  signal  and  determines  the  emergency  situation.  Commutation  of  the 
force  signal  guarantees  that  it  will  be  disconnected  for  the  time  of 
completion  of  transport  operations.  Compensation  of  the  gravitational 
component  of  the  force  signal  eliminates  the  effect  of  the  weight  of  the 
tool  as  a  function  of  its  orientation  in  space.  The  program  module  for 
additional  processing  of  the  force  signal  operates  under  control  of 
modified  ARPS  language  of  the  RM-01  robot. 

The  tests  of  the  developed  force  sensitization  completely  confirmed  the 
operating  principles,  embedded  in  all  its  functional  modules.  A  set  of 
technical  documentation,  which  permits  duplication  of  it  for  use  in 
trimming  complexes  based  on  the  RM-01  robot,  has  been  developed  for  the 
sensitization  system. 
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[Textl  The  problem  of  synthezing  a  system  for 
recognition  of  the  final  number  of  images  given  by 
their  brightness  functions  is  considered.  A  procedure 
for  determination  of  the  location  of  the  sensors  and 
of  designing  a  specialized  computer  for  recognition  of 
plane  images  is  described.  Examples  of  designing 
robotic  vision  systems,  based  on  the  given  recognition 
system,  are  outlined. 

One  must  solve  diverse  recognition  problems  when  developing  sensitized 
industrial  robots.  These  problems  are  ordinarily  solved  by  using 
integral  [1-5]  or  contour  [6,  7]  methods.  However,  integral  methods 
require  reading  of  the  entire  image,  while  contour  methods  require 
reading  the  contour  images,  which  leads  in  the  first  case  to  complex 
technical  realizations  and  in  the  second  case  to  low  noise  stability. 

The  proposed  method  is  applicable  for  recognition  of  a  finite  number  of 
images,  given  by  their  brightness  functions.  It  permits  one  to 
synthesize  m  <  n  sensors  and  a  specialized  computer,  which,  processing 
data  from  the  sensors,  performs  the  recognition,  by  the  brightness 
functions  of  n  images. 

Operating  principle  of  recognition  system.  Let  on  set  X  be  given  n 
mutually  different  functions  f i :X  -*  R  (i  =  l...n),  further  called 
brightness  functions.  One  of  the  brightness  functions  is  represented. 

It  is  required  to  know  which  function. 

Let  us  denote  by  x  the  point  of  set  X  and  let  us  determine  the  image 
F:X  -*  Rn  by  the  matrix  equation 
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F(X)  = 


/iW 

Jn  (*)  . 


Let  us  call  system  of  points  x 1 . . . xrn  of  set  X  the  separating  system,  if 
the  rows  of  matrix  [F(x> ) . . .F(xm) ]  are  mutually  different. 

Let  us  synthesize  a  recognition  system,  formed  of  sensors  and  a 
specialized  computer,  in  the  next  point  by  matrix  [F(x* ) . . .F(xm) ] ,  and 
let  us  describe  the  operating  principle  of  this  system. 

Let  us  assume  that  x1 . . .xm  is  the  separating  system  and  let  us  show  how 
it  is  proposed  to  solve  the  postulated  recognition  problem.  Let  f  be 
the  imposed  brightness  function.  Let  us  form  a  row  [f(x*)...f(xm)], 
corresponding  to  one  of  the  rows  of  the  matrix  [ F ( x ' ) . . . F ( xm ) ] ,  for 
example,  with  the  i— th  row.  In  this  case  f  -  fi  and,  accordingly,  the 
brightness  function  is  proved. 

Let  us  show  that  there  is  a  separating  system,  formed  of  m  <  n  points, 
in  set  X.  If  the  vectors  F(x> ) . . .F(xm)  form  the  basis  of  the  linear 
subspace,  generated  by  all  vectors  of  type  F(x)  at  points  x* . . .xra,  the 
system  of  points  xl . . .xra  will  be  the  separating  system.  Let  us  present 
the  proof  by  the  method  of  the  opposite.  Let  us  assume  that  rows  with 
numbers  i  and  j  coincide  in  the  matrix  A  =  [F(x* ) . . «F(xra) ] .  There  is 
the  vector  y  6  Rm  at  which  F(x)  =  Ay  for  any  x  6  X  by  specific  points 
x'...xra.  Therefore,  fi(x)  =  f j (x)  at  all  x  €  X,  and  this  can  not  be, 
since  the  brightness  functions  were  assumed  mutually  different.  This 
means  that  the  system  of  points  x1 . . .xm  is  the  separating  system. 

The  inequality  m  <  n  is  obvious. 

Technical  realization  of  recognition  system.  The  values  of 
f ( x 1 ) . . . f ( xm )  of  brightness  function  f  must  be  calculated  at  points 
x1 . . .xm  of  the  separating  function  to  realize  the  described  recognition 
method.  This  can  be  accomplished  by  using  sensors  that  operate  in  the 
following  mode: 


c  _  I  0  at 

r  1  1  at 

where  q  =  l...k,  /l  is  an  input  signal,  Xi<a2<  •  •  •  are  the 

quantification  levels,  and  are  output  signals.  The  necessary 

sensors  are  designed  by  matrix  columns  [F(x* ) . . .F(xra) ]  and  includes 

selection  of  quantification  levels  X(  . .  . and  of  their  number  kJ 

( j  =  1. .  .m) . 
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Let  I-1!  <  :A2  <  <  S-  ,+i  be  the  different  elements  of  column  F(xi). 

Let  us  assume  that  X/  =  {*/...  X^y=a^  as  the  quantification  levels  of 

the  j-th  sensor.  The  number  of  kJ  is  a  unit  less  than  the  number  of 
different  elements  in  column  F(xJ ) .  Let  us  denote  the  input  signal  by 

/U  and  let  us  denote  the  output  signals  by 

A  block  diagram  of  the  recognition  system  is  shown  in  Figure  1 ,  where 
D1 . . .D®  are  sensors  and  SV  is  a  specialized  computer. 

The  specialized  computer  performs  the  function  of  decoder:  if  signals 
l1  =  fi  (x1 )  J .  ./lm  =  fi(xra)  are  fed  to  the  inputs  of  sensors  D1 . . .D®,  the 
i-th  output  signal  of  the  computer  is  equal  to  one,  and  all  remaining 
signals  are  equal  to  zero . 

To  design  this  computer,  let  us  be  given  the  constant 


n  _  f  I  at  f’fx'XkJ,, 

L  irj~  \  0  at  f-ixi)>Kiq, 

and  here  i  =  1 . . .n,  j  =  1 . .  .m,  and  q  =  1 . ♦ .kj . 

The  computet  that  realizes  the  output  signals  ♦ . . rn  by  the  formula 

ti=  n  n  (b-.f&C'iti) , 

y=*t  (7-i 

where  i  =  l...n,  vdiile  ©  is  addition  modulo  2,  operates  in  the  required 
mode. 

It  should  be  noted  that  a  perturbed  function  f  +  g,  where  g  is  noise, 
rather  than  an  ideal  brightness  function  f  is  fed  to  the  inputs  of  the 
sensors  in  engineering  problems.  Therefore,  the  problem  of  protecting 
the  recognition  system  against  noise  occurs.  This  problem  can  be  solved 
by  selecting  the  values  of  the  quantification  levels  of  the  sensors. 

Let  us  be  given  the  number 

£  =  4\  m.in  min .  -  f4) 

and  let  us  use  as  the  quantification  levels  the  values  /ljq  =  fi jq  +  e, 
vdiere  j  =  l...m,  and  q  =  l...ki.  The  recognition  system  is  protected 
against  noise  that  does  not  exceed  the  number  e  in  absolute  value  at 
these  quantification  levels. 
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Moreover,  the  brightness  functions  are  not  equally  probably 

manifested  in  engineering  problems;  therefore,  the  problem  of  providing 
the  required  probability  of  correct  operation  of  the  recognition  system 
arises.  This  problem  is  solved  by  increasing  the  number  of  sensors  and 
of  modifying  the  specialized  computer.  Since  the  description  of  the 
general  solution  is  rather  cumbersome,  it  will  be  illustrated  further  on 
the  example  of  a  specific  problems. 

Examples  of  synthesis  of  recognition  system  [8].  The  described 
recognition  system  was  technically  realized  in  different  modifications 
for  recognition  of  large  parts  of  the  K701  tractor,  for  determination  of 
the  orientation  of  nonmagnetic  conducting  parts,  when  designing  the 
robotic  diesel  painting  complex,  to  determine  the  position  of  the  parts 
in  the  cold  stamping  section,  and  when  designing  the  recognition  systems 
and  a  nondestructive  quality  system.  Photosensors,  inductive  and 
strain-gage  sensors,  and  gamma-radiation  sensors  were  used  in  developing 
these  systems.  Specialized  computers  were  realized  in  all  systems  on 
microcircuits  of  series  K153  and  K155. 

The  proposed  method  will  be  illustrated  on  examples  of  design  of  flat 
part  recognition  systems. 

Example  1.  Flat  part  recognition  systems.  Let  us  determine  the 
brightness  function  f(x)  of  a  flat  image  G  C  R2 .  Let  D(x,  r)  be  a 
circle  with  center  at  point  x  and  with  radius  r,  and  let  fi( A)  be  the 
area  of  section  A.  The  brightness  function  f  of  image  G  is  given  on  set 
X  =  R2  and  is  determined  by  the  formula 


/(*)  = 


[j.  [o  n  d  (•*.  r)l 

n  r2 


Constant  r  is  the  radius  of  a  circular  section,  from  which  the  sensor 
perceives  information.  The  brightness  function  f  is  similar  to  the 
characteristic  function  of  set  G,  equal  to  one  at  G  and  equal  to  zero 
outside  G,  at  small  r. 

The  brightness  fucntion  f  is  measured  by  sensors  operating  in  the 
following  mode: 

f  0  at  i  , 

_  1 

°“[  1  at 


Let  us  consider  examples  of  synthesizing  a  system  for  recognition  of  the 
positions  of  a  fixed  part  (variant  1),  parts  with  one  degree  of  freedom 
(variant  2),  and  parts  with  two  degrees  of  freedom  (variant  3).  The 
part  has  the  shape  of  an  angle  in  all  cases. 
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Variant  1.  If  the  part  1  is  transported  by  belt  conveyer  2  (Figure  2, 
a) ,  It  begins  to  be  rotated  due  to  friction  forces  upon  approach  to  stop 
3.  It  is  clamped  by  jaws  4  after  completion  of  rotation  for  fixation  of 
the  part  and  occupies  one  of  three  positions,  shown  on  Figure  3,  a.  In 
this  figure,  axis  X  proceeds  along  the  stop,  while  axis  Y  proceeds  along 
the  left  jaw  of  the  grab.  The  superimposed  contours  of  the  positions  of 
the  part  are  shown  in  Figure  4,  a. 


KEY: 

1.  Specialized  computer 


Figure  2 

Let  f i ,  f 2  and  f3  be  brightness  functions  of  plane  images,  shown  in 
Figure  3,  a.  Let  us  assume  that  constant  r  is  sufficiently  small.  In 
this  case,  points  x1  and  x2 ,  shown  in  Figure  4,  a  by  the  asterisks,  form 
a  separating  system. 

Matrix  [F(x')»  F(x2)]  will  have  the  form 
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Accordingly,  the  specialized  computer  should  realize  the  following 
output  signals  to  recognize  the  positions  of  the  part 


T|  —  6  j 

T2=  (6i©1)62, 


t3=6i(62®1). 


If  P  is  the  probability  of  correct  operation  of  the  sensor,  the 
probability  of  correct  operation  of  the  synthesized  recognition  system 
will  be  P|  =  P2 . 


Let  us  place  additional  sensors  at  points  x3-x6,  shown  in  Figure  4,  a  by 
the  asterisks.  The  specialized  computer  that  realizes  the  following 
output  signals 

Ti=  ®  n  (6j9akj), 

jwm\ 


where  i  =  1,  2,  3,  while  11  x  6  matrices  c!  ,  c2  ,  and  c3 
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guarantees  the  probability  of  correct  operation  of  the  recognition 
system  /^  =  p6  +  6p5(l — pj  +4/r(l—  p)2.  At  p  >  0*54,  the  probability  P2  >  Pi  . 
Therefore,  the  second  recognition  system  operates  more  correctly. 


Figure  3 

Variant  2.  Let  part  1,  as  in  the  first  case,  be  transported  by  a  belt 
conveyer  2  (Figure  2,  b),  but  after  rotation  near  stop  3  is  not  squeezed 
by  the  jaws.  To  determine  the  position  of  the  part,  let  us  place 
sensors  5  on  push  rod  6,  having  the  possibility  of  forward  displacement 
parallel  to  the  stop  3.  The  mutual  arrangement  of  points  xl ,  x2,  and  x3 

of  the  sensors  is  shown  in  Figure  4,  b.  The  matrix  F(x2),  -F(*3)] 

has  the  form 


']  1  1" 
0  0  1 
1  0  1 


therefore,  the  output  signals  of  the  specialized  computer  are 


T1  =  616263, 

t2=  (6i©1)  (62©  1 ) 63, 
t3  =  6i  (62©  1)63. 
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We  note  that  the  sensors  placed  at  points  x1  and  x2  determine  the  number 
of  the  position  of  the  part,  while  the  sensor  located  at  point  x3 
determines  its  displacement  with  respect  to  the  stop. 

Variant  3.  Let  the  part  be  arbitrarily  located  on  a  plane.  Let  us  link 
coordinate  system  XY  to  the  plane  (Figure  3,  b) .  Let  us  measure  the 
coordiantes  (Xi ,  Y|)  of  the  topmost  and  (X2,  Y2)  bottommost  vertices  of 
the  part  (the  possible  ambiguity  is  easily  eliminated).  This 
measurement  can  be  made,  for  example,  by  using  an  integral  photodetector 
array. 


Figure  4 

Six  possible  positions  of  the  part,  when  its  different  vertices  are  the 
upper  or  lower,  are  shown  in  Figure  3,  c. 

Knowing  coordinates  (Xi ,  Y()  and  (X2,  Y2),  let  us  transform  the 
coordinate  system  so  that  the  equalities  V1  =  71  =  y2  =  0  are  fulfilled  in 

<v 

new  coordinate  system  XY.  The  contours  of  all  six  images  can  then  be 
combined  in  the  new  coordinate  system,  as  shown  in  Figure  4,  c. 

-v 

The  points  of  placement  of  sensors  x 1 ,  x2 ,  and  x3  are  shown  by  the 

circles  in  Figure  4,  c  in  coordinate  system  XY.  Let  us  apply  the 

inverse  transform  of  coordinates  to  points  x1 ,  x2,  and  x3  to  determine 
the  coordinates  of  cells  x1,  x2,  and  x3  in  the  photodetector  array,  from 


71 


which  the  information  must  be  read.  Regardless  of  the  transformation  of 
coordinates,  matrix  [F(x*)>  F(x2),  F(x3)]  will  have  the  form 


’I  1  0“ 

1  0  0 
0  1  0 
oil 
0  0  1 
.0  0  0. 

Accordingly,  the  specialized  computer  should  realize  the  following 
output  signals 


t,  =  M2(63©1),  t2  =  6i(62©1)(63©1), 

t3=  (6i©1)62  (  63©1),  X4=  (6i©1)6283, 

t5=  ( 6i ©  1 )  (62©1)63,  t6=  (6i©  1)  (62©1)  (63ffil). 

Example  2.  Recognition  system  for  three-dimensional  parts.  Let  us 
determine  the  brightness  function  f  of  a  three-dimensional  image  G  C  R3 . 
Let  X  =  R3  x  S2,  where  S2  is  a  unit  sphere.  Points  X  of  set  R3  x  S2  are 
pairs  (xi,  x'2),  where  xi  is  an  arbitrary  vector  and  X2  is  an  arbitrary 
unit  vector,  which  controls  the  direction.  Let  D(x,  r)  be  a  circle  with 
center  at  point  xt  and  with  radius  r,  located  in  plane  a,  perpendicular 
to  vector  x2  ,  and  let  ?r(x,  ,  G)  be  the  orthogonal  projection  of  the  body 
G  onto  plane  a.  The  brightness  function  f  of  the  image  G  is  determined 
by  the  formula 


/(*)  = 


OnOUi,  ri] 
71  r2 


The  sensors  that  measure  brightness  function  f  operate  in  the  previously 
described  mode. 


Figure  5 
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Figure  6 

Let  us  present  a  system  for  recognition  of  12  parts  of  the  K701  tractor, 
moved  by  an  overhead  push  conveyer.  Some  of  these  parts  on  the  conveyer 

are  shown  in  Figure  5.  The  points  *'=  (x,1.  *2'),  x2=  (xx2,  x22),  x3=  (x,3,  xi), 

r!=  (xi4,  -ly'j ,  xJ=  (X]  \  a:;),,  which  form  the  separating  system,  are  shown  in 

Figure  6  [ 8 ] . 
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MATHEMATICAL  MODEL  OF  POINTER  STROBE  TELEVISION  SYSTEM  UNDER  NOISY 
CONDITIONS 
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[Article  by  V.  M.  Artemyev  and  K.  V.  Glagolev ,  Moscow  State  Technical 
University  imeni  N.  E.  Bauman] 

[Text]  A  mathematical  model  of  a  strobe  pointer 
television  system  is  proposed.  The  effective  noise  on 
the  operation  of  the  system  in  the  steady  mode  is 
considered.  The  operation  of  a  two-circuit  system 
with  connection  between  channels  with  a  light¬ 
scattering  layer  in  the  working  space  was  analyzed 
within  the  model  by  the  statistical  linearization  v 
method . 

A  strobe  pointer  system  can  be  used  in  robotic  visual  sensitization 
problems  as  the  sensor  of  information  about  the  coordinates  of  the 
object  under  noise  conditions,  caused  by  the  presence  of  a  light¬ 
scattering  layer  in  the  working  space  of  the  system  [ 1 ] .  A  narrow 
visual  field  permits  one  to  reduce  the  effect  of  scattered  emission  and 
thus  to  enhance  the  contrast  of  the  tracked  object.  The  mathematical 
model  of  an  optical  angle  coordinator,  based  on  the  strobe  pointer 
principle,  is  proposed  in  the  article  for  operating  in  the  presence  of  a 
light-scattering  layer  in  the  working  space  of  the  coordinator.  The 
system,  shown  in  Figure  1,  operates  on  the  closed  cycle  and  has  two 
control  circuits — an  external  circuit,  closed  by  the  television  camera 
rotational  drives,  and  an  internal  circuit,  closed  by  the  strobe 
generator.  The  problem  of  the  internal  circuit  is  to  keep  the  strobe 
pointer  on  the  image  of  the  object. 

The  length  of  the  image  along  the  horizontal  and  vertical  is  denoted  by 
0X  and  0y  (Figure  2).  In  modern  systems,  the  dimensions  of  the  "window" 
are  automatically  adjusted  to  the  dimensions  of  the  image;  therefore, 
let  us  further  assume  that  the  horizontal  length  of  the  window  is  equal 
to  Ox,  while  the  vertical  length  is  0y .  The  errors  lx  and  ly  between 
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the  axes  of  the  "window"  and  the  axes  of  the  image  will  be  small,  if  the 
center  of  the  image  is  located  near  the  origin  of  the  television  screen. 

The  nature  of  the  signal  at  the  output  of  the  television  receiver  within 
one  line  is  shown  in  Figure  3.  The  signal  from  the  source  is  perceived 
on  some  background,  the  mean  level  of  which  is  denoted  by  U^.  The 

amplitude,  i.e.,  the  maximum  value  of  the  signal  with  respect  to  U^,  is 

denoted  by  A.  Each  signal  in  the  line  is  normalized  by  a  limiter,  at 
the  output  of  which  a  square-wave  signal  is  formed  with  amplitude  Uo , 
the  level  of  which  is  determined  in  the  background  level  estimation 
circuit  [2]. 


Figure  1.  Structure  of  Television  System: 

TVK— television  camera;  TVP— television  receiver 
with  image  scanning;  M — monitor;  OD — optical 
discriminator;  UF — amplifying  and  filtering  elements; 
SZ — delay  circuit;  GS— strobe  generator;  P— delay  time 
to  voltage  converter;  PPX  and  PPy — drives  for  rotation 
of  television  camera  through  coordinates  x  and  y 


KEY: 

1.  Television  camera 

2.  Television  receiver 

3.  Drives  for  rotation  of 
television  camera  through 
coordinate  y 

4.  Drive  for  rotation  of 
television  camera  through 
coordinate  x 


5.  Monitor 

6.  Delay  time  to  voltage 
converter 

7 .  Strobe  generator 

8.  Optical  discriminator 

9 .  Delay  circuit 

10.  Amplifying  and  filter 
elements 


Each  of  the  four  strobes  shown  in  Figure  2  "cuts  out"  part  of  the  image 
falling  into  this  strobe.  Four  voltages  Ugi,  Uga,  Ug 3 .  and  Ug4, 
proportional  to  the  mean  values  of  the  quantified  signals  impinging  on 
the  corresponding  strobes,  are  formed  at  the  end  of  the  frame. 
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Figure  2 .  Screen  of  Television  Monitor  With 
Coordinate  System  x,  y: 
a,  b — coordinate  system  of  strobe  pointer; 

1-4 — numbers  of  strobes 

For  example,  let  errors  lx  and  ly  lie  in  the  limits  0  <  |lx|  <  O.50x  and 
0  <  | ly |  <  0.56y.  It  is  this  case  of  the  image  that  is  in  Figure  2. 

The  averaged  voltages  in  the  strobes  are  then  equal  to 


Ug\—  (  UqITk)  0,50k  •  0,5%; 

Ug2=  (UoITk)  (0,50* — k)0,5%;  |  d) 

Ue 3=  (Uo!Tk)  (O,50k — /*)  (0,5% — f  (ly) ) ; 

U  g4==  (Uol'f  k)  O,o0k  (0,5%  f{ly))>  ^ 

where  Tg  and  Tq  are  the  time  of  the  frame  and  line  scanning, 

respectively,  ny  =  0y/To  is  the  number  of  lines  in  the  "window,"  and 
f(ly)  is  the  function  ly,  which  denotes  the  maximum  integer,  close  to 
the  ratio  ly/To .  The  error  voltages  are  generated  by  the  laws 


Ux  —  (Ugi-rUgi)  —  {Ug2+Ug^)t 

Uy  =  (Ugl  +  Ug  2) — (Ug2+Ug4). 


Generalizing  expressions  similar  to  (1)  and  using  (2),  one  can  write  the 
following  formulas  for  the  discrimination  characteristics  (DKh): 


It  is  obvious  that  the  continuous  error  function  lx  is  located  along  the 
horizontal  (line)  coordinate  of  the  discrimination  characteristic,  while 


77 


the  dependence  of  the  discrimination  characteristic  on  the  error  is  a 
discrete  step  function  along  the  vertical  (frame)  coordinate. 

The  conversion  factors  KP  (the  characteristic  slope)  of  the  channels  of 
the  discriminator  are  found  from  (3) 


(ny  -  |  /  (/,) 


Un 


T  k  ‘  c 


(6 


x~\h\ 


(4) 


It  follows  from  (4)  that  there  is  a  dependence  between  the  channels  of 
the  discriminator  through  their  conversion  factor,  which  is  manifested 
stronger,  the  greater  the  value  of  the  errors. 

Let  us  now7  consider  the  effect  of  noise  on  the  operation  of  the  system 
in  the  steady  mode.  The  presence  of  the  noise  component  of  the  signal 
at  the  input  fo  the  receiver  is  shown  in  Figure  3.  Let  the  mean  noise 
output  (dispersion)  by  time  with  respect  to  the  background  level  be 

equal  to  P^.  The  effect  of  noise  reduces  to  variation  of  the  moments  of 

response  of  the  threshold  device  t i  and  t2 •  As  a  result  of  the  internal 
noise  of  the  receiver,  the  lengths  of  the  pulsed  signals  in  the  strobes 
become  random  and  random  components  appear  in  voltages  Ug l -Ug  4  that  also 
lead  to  the  appearance  of  fluctuating  components  at  the  input  of  the 
discriminator.  Let  us  find  the  value  of  the  spectral  densities  of  these 
fluctuations  at  zero  errors  lx  =  ly  =0.  Consideration  of  the 
instability  of  the  response  of  the  threshold  circuit  is  described  in 
[3],  The  dispersion  of  the  instability  of  the  length  of  the  pulse  front 
is  approximately  equal  to  D  —  P^/ a 2 ,  where  a  is  the  characteristic  slope 

of  the  signal  pulse  in  the  region  of  threshold  level  Uo  [4] .  This  value 
can  be  taken  into  account  in  the  following  manner.  With  triangular 
shape  of  the  signal  with  amplitude  A,  a  =  2A/ •  Variation  of  the 
characteristic  slope  with  a  different  shape  of  the  pulse  is  denoted  b\ 
coefficient  ac  a  —  2aA/$x*  The  dispersion  of  the  instability  of  the 
pulse  length  will  then  be 


a -  4a*  An- 


Let  us  call  the  signal/noise  ratio  q  =  A2/Pffl.  The  dispersion  is  then 


D  —  B.x2/4a2q. 
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Since  there  is  0.5ny  lines  in  each  strobe  and  since  the  front  of  each 
pulse  in  the  line  fluctuates  statistically  regardless  of  other  lines, 
the  dispersion  of  each  of  the  voltages  Ugt~Ug4  will  be  equal  to 


0,5nyDUo 

T~k 


S*’l/  f  K 


Figure  3.  Signal  at  Output  of  Receiver 
Within  One  Line 


When  voltage  Ux  and  Uy  are  formed,  all  four  signals  Ugi  are  added 
algebraically;  therefore,  the  dispersion  of  the  amplitude  of  each  of  the 


voltages  will  be  fourfold  greater,  i.e. 


'la^qT  K 


Since  the  error 


voltages  Ux  and  Uy  are  filled  for  the  period  T  of  scanning  the  next 

frame,  the  fluctuation  component  forms  a  sequence  of  square-wave  pulses 
with  statistically  independent  amplitudes  (dispersions  D) . 


For  the  subsequent  narrowband  actuating  part  of  the  system,  this  process 
can  be  considered  white  noise  with  the  spectral  density 


A;=  D„,T..= 


‘^qT  y. 


If  the  value  of  the  conversion  factor  of  the  line  channel  is  written  in 
the  form 


k 


ex 


l/Uy)|  \ 

ny  j 


and  if  the  multiplifier  Uony/TR  in  the  conversion  factor  of  the 

actuating  part  of  the  system  is  taken  into  account,  the  parameters  of 
the  discrimination  characteristic  of  the  line  channel  will  be  equal  to 
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(5) 


In  similar  fashion,  for  the  frame  channel 


/J  "  T  Tl 

1  1  fi  .  A/  —  v  K  0 

1  —  «  »  y 


Or 


2a-q 


(6) 


l 

. 


Figure  4.  Block  Diagram  of  Television  Discriminator 

The  scattering  layer  in  the  working  space  of  the  system  leads  to  a 
decrease  of  the  optical  signal  at  the  input  to  the  system  and, 
therefore,  to  a  decrease  of  the  signal/noise  ratio  q.  If  this  decrease 
is  taken  into  account  by  introduction  of  coefficient  j)  =  q/qo,  where  qo 
is  the  signal/noise  ratio  in  the  absence  of  a  scattering  layer,  the 
spectral  densities  of  the  noise  components  in  the  line  and  frame 
channels  will ,  according  to  formulas  ( 5 )  and  ( 6 ) ,  be  inversely 
proportional  to  coefficient  /?: 


j\/  ——XJL 
*  f 


N. 


yo. 


(7) 


where  Nxo  and  Nyo  are  the  spectral  densities  of  noise  in  the  line  and 
frame  channels,  respectively,  are  determined  by  formulas  (5)  and  (6)  at 
q  =  qo. 

Let  us  carry  out  further  analysis,  assuming  that  a  =  const.  If  the 
useful  signal  is  determined  by  a  non-scattering  component  of  source 
emission,  the  characteristic  slope  of  the  signal  pulse  will  not  be 
dependent  on  scattering  and  coefficient  a  can  be  considered  constant. 

Of  course,  when  the  optical  coordinator  operates  on  the  scattering 
component,  coefficient  a  will  vary  as  scattering  increases  in  the 
located  space  due  to  variation  of  the  intensity  and  due  to  an  increase 
of  the  spreading  surface  of  an  image  of  a  light  spot.  However,  this 
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does  not  alter  the  methods  of  analysis  of  errors  in  the  system,  so  that 
one  can  assume  that  variation  of  coefficient  0  includes  not  only 
variation  of  the  signal/noise  ratio,  but  also  variation  of  the 
characteristic  slope  of  the  signal  pulse  due  to  scattering.  Comparison 
of  formulas  (4)  and  (7)  shows  that  an  increase  of  the  error  is 
equivalent  to  a  decrease  of  signal/noise  ratio  for  a  strobe  pointer 
television  system. 

As  one  can  see  from  expressions  (5)  and  (6),  there  is  a  relationship 
between  the  discriminator  channels  through  cofactors  (1  -  |f(ly)|/ny) 
and  (1  -  |lx|/#x).  As  a  result,  the  block  diagram  of  the  television 
discriminator  assumes  the  form  shown  in  Figure  4,  where  £x  and  £y  are 
white  noise  processes. 

A  two-circuit  interconnected  system  can  be  analyzed  in  two  steps.  In 
the  first  step,  the  properties  of  only  the  internal  circuits  are  studied 
with  regard  to  the  relationship  between  them,  and  in  the  second  step, 
the  system  as  a  whole  is  studied.  Let  us  consider  the  method  of 
calculating  the  fluctuation  errors  in  the  system  shown  in  Figure  4, 
under  the  following  conditions:  the  work  proceeds  on  linear  sections  of 
the  discrimination  characteristic,  the  number  of  television  scanning 
lines,  required  for  the  dimensions  of  the  image  along  the  vertical,  is 
greater  than  20,  which  permits  one  to  disregard  quantization  of  the 
signal  by  the  level  in  the  frame  channel  of  tracking  and  to  assume  that 
its  characteristics  are  continuous  functions.  Let  us  use  the 
statistical  linearization  method  to  analyze  this  system  in  the  steady 
mode  [ 5 ] . 

Since  only  the  fluctuation  effects  £x  (t )  and  ( y ( t )  are  taken  into 
account  in  the  form  of  white  noise,  tracking  errors  lx(t)  and  ly(t)  will 
be  random  processes  with  zero  mean  value  and  variance  Djx  and  Diy.  Let 
us  find  the  statistical  linearization  coefficients  for  a  system  with  a 
block  diagram  shown  in  Figure  4  by  the  equations 

#  i  .  /  2£)/v 

kox  ~  1  —  |/  n02  aMa'&0x(Dly)i 

t  ,  /”2ozr 

k oy  ~  u02  =  ^0y  (Ax). 

in  which  it  is  taken  into  account  that  the  mean  values  of  the  tracking 
errors  are  equal  to  zero  in  the  expressions  for  normal  probability 
densities . 

Using  the  Parseval  integral,  let  us  write  the  variance  of  errors  in  the 
steady  mode  [ 6 ] : 


Dlx — Nxlxt  Dly  Nyly, 


(9) 
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where  Ix  and  Iy  are  Parseval  integrals  for  transfer  functions  of  the 
errors  of  the  line  and  frame  channels,  respectively.  Relations  (9)  form 
a  system  of  two  equations,  connected  by  coefficients  kox  and  koy.  Since 
the  Parseval  integrals  Ix  and  Iy  can  be  found  in  explicit  form,  (9)  with 
regard  to  (8)  forms  a  system  of  nonlinear  algebraic  equations,  which  is 
solved  numerically  or  graphically. 

Further  analysis  concerns  study  of  a  television-optical  location  system, 
developed  by  V.  V.  Vasilyev  [1]  and  having  the  following 
characteristics:  transfer  functions  of  the  internal  circuits 


Ki(P)rkilP*  Ks(P)  =k2lp; 


transfer  functions  of  the  external  circuits,  containing  DC  electric 
motors 


a :«(*)- 


*3 

p(l+pT3)  ’ 


K<(P)  = 


*4 

PU+pTJ  ' 


As  the  experiments  showed,  the  mean  square  values  of  the  tracking  errors 
are  on  the  order  of  several  angular  minutes.  This  means  that  the  signal 
to  noise  ratio  under  conditions  of  a  clean  atmosphere  is  qo  >  5.  Based 
on  (8)  and  (9)  for  internal  circuits,  we  arrive  at  the  following  system 
of  algebraic  equations: 


(10) 


here  Dxo  and  Dyo  are  the  variants  of  the  errors  in  the  channels  of  the 
system  without  regard  to  cross-correlation  between  channels  and  in  the 
absence  of  a  scattering  layer. 

Graphical  solution  of  system  of  equations  (10)  is  presented  in  Figure  5. 
We  find  solutions  of  them  at  the  intersection  points  of  A  and  B  of 
curves  1  and  2.  According  to  the  condition  of  the  problem,  the  desired 
point  is  point  A,  since  the  values  of  Dix  and  Diy  increase  only  for  this 
point  as  Dxo  and  Dyo  increases.  At  some  values  of  Dxo  and  Dyo,  curves  1 
and  2  cease  to  intersect,  and  this  means  that  the  system  lost  stability 
due  to  cross-correlation  of  the  channels. 
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Figure  5.  Graphical  Solution  of  System  (10): 

1 — graph  of  function  Dlx  =  f(Dly);  2 — graph  of 
function  Dly  =  f(Dlx) 

It  is  more  convenient  for  practical  calculations  to  use  the  relative 
values  of  the  variance  of  errors  in  formulas  (10) 


(ID 


where 


Ax  —  Dxo/e*2,  Ay  —  Dyo/Qy2. 


Let  us  present  the  analytical  solution  of  equations  (11)  for  the  special 
case  when  the  parameters  of  the  channels  are  identical .  It  is  obvious 
that  here  Ax  =  Ay  =  Ao  and  Sx  =  Sy  -  6.  At  some  value  of  6  and  at  given 
parameter  Ao ,  we  calculate 


/= 


and  we  find  coefficient  as  a  function  of  the  variance  of  the  random 
component  of  the  error  in  the  channel: 


P=(/6)-1; 


(12) 
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6  is  determined  by  formula  (11).  The  dependence  of  the  variance  of  the 
error  in  the  channel  on  coefficient  ft ,  calculated  for  different  values 
of  parameters  Ao,  is  shown  in  Figure  6.  The  errors  in  the  system 
increase  strongly  at  the  same  values  of  ft  as  Ao  increases .  Analysis  of 
expression  (12)  shows  that  disruption  of  tracking  in  the  channels  of  the 
internal  circuit  begins  at  variance  of  error  in  the  channel  of  8  - 
-  2tt/9A2 .  Therefore,  curve  S (ft)  breaks  off  at  Ao  =  0.1  in  Figure  6, 
when  ft  -  0.43,  i.e.,  at  the  moment  of  disruption  of  tracking,  and  an 
increase  of  the  variance  of  the  error  begins  even  with  an  increase  of 
the  signal/noise  ratio  (the  dashed  curve  in  Figure  6). 


fi- 


Figure  6.  Dependence  of  Variance  of  Random  Component 
of  Error  in  Channel  for  One-Circuit  (Curves  1-3)  and 
Two-Circuit  (Curves  4,  5)  Systems  on  Coefficient  ft) 
the  Dashes  Correspond  to  Experimental  Results 

The  method  described  earlier  is  applicable  for  analysis  of  a  two-circuit 
system— only  the  structure  of  the  transfer  function  is  complicated.  The 
transfer  function  of  the  error  of  the  first  channel  has  the  form 


KM 


_ *1*3 _ _ 

*rjc*i*3  +  .P*<ur*i  +/’2  (1  +  *o,*.7' 3)  4-  p37 3 


Let  us  find  the  Parseval  integral  of  this  function  and  let  us  find  the 
next  expression  for  the  variance  of  the  error  of  the  line  channel: 


=  W, 


*i  (1  4  knxk\Tz) 


2*; 


oxi^ox^i^s  ^3^3  "b  1) 


The  expression  of  the  variance  of  the  error  for  the  frame  channel  is 
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Dlv  =  N. 


*4(1+  *oy  m 


y  2k0y(%ktTi~kiTi  +  ]) 


These  equations  can  be  solved  graphically.  Let  us  carry  out  specific 
calculations,  assuming  that  the  parameters  of  the  channels  are 
identical;  therefore,  Dix  =  Diy. 

If  there  is  no  light  scattering  (4=1)  and  if  there  is  no  relationship 
between  the  channels  (kox  =  1)>  the  variance  is 


Dxo  Nxq 


*a(l  +  *.7,3) 

2  (*i73 —  *3r3  +  1) 


Taking  (7)  and  (8)  into  account, 


_  Dix  _  1 
x  Dxo  p 


(I  +  *nr*i73)  {k^7\  —  *3 7* 3  -f-  1) 
k0x  (koxki  TS  ~  k3T3  +  1 )  ( 1  +  *3  T 3) 


(13) 


(14) 


where  hx=.Dxo  —  I®*2  >  the  same  as  for  a  one-circuit  system. 

One  can  find  a  clear  relationship  between  kox  and  from  equation  (14) 
at  some  parameter  of  Ax,  and  one  can  calculate  kox(^x)  for  the  selected 
value  of  Sx ,  while  one  can  calculate 


(1  +  *o.r*iT3)  (*iT3  —  ^3r3  -f  1 )  . 

k0X(k0ATz  -  kJz  +  !)  0  +  Vs)  < 15) 


— the  dependence  of  the  variance  of  the  error  on  the  coefficient  of 
statistical  linearization  in  the  absence  of  a  scattering  layer  at  some 
parameters  of  ki ,  k3,  and  T3  for  the  derived  value  of  kox •  According  to 
(13), 


P=f/6*, 


(16) 


here  Sx  is  the  earlier  selected  value  of  the  relative  variance  of  the 
error.  Relation  (16)  establishes  the  dependence  between  the  relative 
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variance  of  the  error  and  the  coefficient  of  variation  of  the 
signal/noise  ratio  and  permits  one  to  calculate  function  Sx{  3) • 

Let  us  perform  specific  calculations  at  k]  =  500  s'1,  k3  =  200  s'1,  T3  - 
=  0.2  s,  Ax  =  10'3  and  Ax  =  10'2.  The  functions  Sx(3)  for  a  two-circuit 
tracking  system  is  shown  in  Figure  6  by  the  dash-dot  lines.  The  nature 
of  function  6X(3)  is  the  same  as  that  for  a  one-circuit  system;  the 
variance  of  the  error  increases  rapidly  as  Ax  increases.  The  greater 
Ax,  the  earlier  disruption  of  tracking  begins  and  it  is  expressed  in  the 
form  of  an  increase  of  the  variance  of  the  error  as  coefficient  3 
decreases.  The  errors  at  the  same  values  of  3  are  greater  for  a 
two-circuit  system  than  in  a  one-circuit  system  at  Ax  =  A0.  Moreover, 
system  of  equations  (14)  and  (15)  generally  has  no  solution  that 
satisfies  the  physical  meaning  of  the  problem  at  Ax  =  10' 1 .  The  values 
of  formula  (16)  are  greater  than  unity  in  this  case. 

It  is  shown  in  Figure  6  that  the  values  of  are  small  at  3  -  1>  i»e., 
the  relationship  between  channels  in  the  system  is  weak  in  the  absence 
of  a  scattering  layer.  For  example,  if  Ax  =  10'3,  then  kox(^  =  D  = 

=  1.06.  Thus,  the  relationship  of  the  channels  upon  transition  from  a 
one-  to  a  two-circuit  system  is  intensified. 

The  mathematical  model  of  a  tracking  window  television  system  was 
checked  experimentally  on  a  laboratory  unit.  6-3N  and  8-3N  neutral 
light  filters,  initially  calibrated  on  a  spectrograph,  were  used  to 
attenuate  the  light  signal.  The  fluctuation  component  of  voltage  was 
taken  by  using  an  oscillograph  at  the  output  of  the  discriminator.  The 
resulting  experimental  values  (the  dashes  in  Figure  6)  are  easily 
described  by  a  model  function  at  Ax  =  1.4*10'2.  We  note  that  the 
probability  of  disruption  of  tracking,  observed  in  experiment,  was  0.5 
at  3  =  0.28  (8-3N  light  filter). 
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[Article  by  V.  I.  Solntsev,  Moscow  State  Technical  University  imeni  N. 

E.  Bauman] 

[Text]  Problems  of  pre-processing  of  multigradation 
images  of  typical  assembly  scenes  are  considered.  A 
context-free  grammar  of  a  linguistic  variable  is  used 
to  formalize  the  concept  "BRIGHTNESS."  It  is 
suggested  on  its  basis  that  the  brightness  scales  be 
converted  and  that  the  images  be  segmented. 

Information  about  the  course  of  the  manufacturing  process  must  be 
obtained  and  processed  under  conditions  of  a  f lexible  automated  plant . 
The  type  of  information  systems  is  selected  as  a  function  of  technology. 
Analysis  of  typical  robotized  assembly  operations  showed  that  a  robotic 
vision  system  (STZ)  has  the  greatest  information  content.  It  permits 
one  to  determine  the  position,  orientation  and  characteristic  dimensions 
of  the  assembly  unit  and  of  individual  parts,  and  also  to  transfer  this 
information  to  the  robot  control  system.  The  assembly  process  can  be 
checked  by  using  the  STZ. 

The  type  of  assembly  position  has  a  number  of  characteristic  features 
from  the  viewpoint  of  using  a  robotic  vision  system.  The  assembly  unit 
ordinarily  consists  of  uniformly  painted  parts.  Assembly  is  performed 
at  the  work  position  (assembly  table,  revolving  table  and  so  on),  having 
a  relief  surface  due  to  the  presence  of  clamping  devices,  fasteners  and 
other  hardware.  It  is  complicated  to  provide  special  lighting  of  the 
parts  under  production  conditions. 

Thus,  the  image  of  the  work  position,  obtained  by  using  optical  sensors 
(vidicons,  CCD-arrays,  and  photoarrays),  has  the  following 
characteristics : 


low  contrast  of  individual  parts  and  of  their  elements , 
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low  contrast  of  the  images  of  the  assembly  unit  and  of  the  work 
position; 

nonuniform  distribution  of  illumination  due  to  selection  of  the  local 
and  general  lighting  system; 

fluctuation  of  the  level  of  illumination  over  time; 

the  presence  of  contrast  interference  on  the  image  (contamination, 
highlights  and  so  on). 

Moreover,  the  assembly  unit  has  a  complex  structure,  which  leads  to  the 
appearance  of  structures  from  uniformly  lighted  regions  embedded  in  each 
other  on  the  image. 

Binary  robotic  vision  systems,  which  have  become  widespread  (the  number 
of  brightness  gradations  is  equal  to  two) ,  do  not  permit  processing  of 
scenes  of  this  complexity.  The  use  of  binarization  methods  (both  with 
fixed  and  with  adaptive  thresholds)  and  the  preliminary  processing  step 
of  the  images  requires  uniformity  and  stability  of  lighting  and  high 
contrast  between  the  " object"  and  "background."  The  resulting 
silhouette  images  contain  only  an  insignificant  part  of  the  information 
about  the  scene.  Thus,  despite  the  low  cost  and  high  speed  of  binary 
robotic  vision  systems,  the  sphere  of  their  application  for  automation 
of  the  assembly  process  is  limited. 

This  problem  can  be  solved  by  using  multigradation  STZ  (the  number  of 
brightness  gradations  is  from  8-16  to  256  and  above).  A  monochromatic 
half-tone  image  is  subjected  to  various  types  of  transformations 
(correction,  improvement,  filtration,  segmentation  and  so  on)  at  the 
preliminary  processing  step  to  intensify  contrast,  to  compensate  for 
noise  and  to  distinguish  its  characteristic  elements. 

A  large  number  of  methods  and  algorithms  for  processing  multigradation 
images  has  now  been  developed  [1-3].  However,  robotized  assembly 
imposes  a  significant  restriction  on  the  speed  of  the  robotic  vision 
system.  All  the  processing  steps  from  input  of  the  image  of  the  scene 
to  recognition  of  it  should  occur  within  1-3  s  [1].  Therefore,  many 
algorithms  used  with  STZ  to  check  assembly  operations  require 
development  of  special  high-speed  hardware,  which  increases  the  cost  of 
the  system. 

Making  a  histogram  of  the  image  frame  by  brightness  is  most  frequently 
used  in  industrial  robotic  vision  systems  for  adaptation  to  variable 
lighting  conditions.  The  threshold  or  interval  of  binarization  is 
determined  by  the  type  of  histogram,  and  the  image  is  further  processed 
the  same  as  in  binary  STZ  with  adaptive  threshold.  This  method  is 
applicable  for  high-contrast  images.  The  brightness  distribution 
histogram  for  images  of  standard  assembly  units  ordinarily  has  no 
characteristic  "dark"  and  "light"  peaks,  which  does  not  permit  one  to 
determine  with  sufficient  accuracy  the  binarization  threshold. 
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Gradient  methods  yield  good  results  in  the  case  of  an  absence  of  noise 
on  the  image  [2].  But  the  low  contrast  results  in  disruption  of  the 
contour,  approximation  and  connection  of  segments  of  which  require  large 
computer  expenditures.  Moreover,  information  is  lost  after  calculation 
of  any  gradient  function  (Sobel  or  Roberts  operator  or  something  else), 
i.e.,  it  becomes  more  difficult  to  determine  to  which  of  the  parts  it 
corresponds  when  determining  the  only  the  contour  of  parts. 

Studies  of  man’s  visual  perception  system  [4]  indicate  the  psychological 
nature  of  perception  of  the  illumination  and  color  distribution 
function.  The  concept  "BRIGHTNESS”  is  constructed  on  the  basis  of 
subjective  impressions  of  man  and  is  dependent  on  the  total  level  of 
illumination  of  a  scene,  attention  to  its  individual  parts,  and  on  a 
number  of  other  factors.  Determination  of  regions  on  an  image,  uniform 
in  illumination,  is  one  of  the  preliminary  processing  problems  of 
multigradation  images  for  robotic  vision  systems.  Segmentation  of 
images  can  be  regarded  as  a  problem  of  making  a  decision  on  the 
affiliation  of  the  element  of  the  image  of  some  region  to  the  level  of 
illumination,  typical  for  it. 

It  is  suggested  that  linguistic  variables,  constructed  on  the  basis  of 
odd  set  theory,  be  used  to  formalize  the  concept  "BRIGHTNESS"  [5].  From 
the  mathematical  viewpoint,  the  degree  of  the  affiliation  of  the  image 
Hi  €  [0,  1]  to  the  following  odd  set  corresponds  to  the  i-th  element  of 
an  image  with  brightness  code  bi 


Hi — fAjin(^i)»  (1) 


where  /^(b)  is  a  function  of  the  affiliation  of  the  odd  set.  The  key 

aspect  is  determination  of  the  type  of  affiliation  functions.  Basic 
linguistic  variables  "LIGHT,"  "GRAY,”  and  "DARK"  were  selected  for  the 
concept  "BRIGHTNESS."  Based  on  analysis  of  expert  estimates  of 
brightness  for  separate  areas  of  images  of  a  number  of  assembly  units 
and  parts,  it  was  determined  that  the  affiliation  functions  have  a 
nature  close  to  quadratic  (Figure  1).  The  concept  "BRIGHTNESS"  can  be 
approximated  by  a  context-free  grammar,  which  is  the  set  of  four 


G=<N,T(N),U,S), 

where  N  is  the  name  of  the  linguistic  variable,  T(N)  is  the  term-set  of 
its  values,  U  is  the  set  of  values  of  the  basic  variable,  and  S  is  the 
syntactic  rule  of  formation  of  the  linguistic  variable.  The  structure 
of  the  grammar  (2)  is  presented  in  Figure  2. 
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Figure  1.  Functions  of  Affiliation  of  Basic 
Linguistic  Variables: 

-  - - "LIGHT"; - "DARK", - "GRAY" 

An  alphabet  of  linguistic  variables  is  formed  by  using  the  proposed 
grammar.  The  number  of  values  of  the  variables  is  dependent  on  the 
postulated  problem  and  should  not  exceed  the  Miller  number  7  ±  2 ,  which 
expresses  the  maximum  number  of  gradations  of  the  parameter, 
simultaneously  stored  in  the  human  memory  [4],  Expression  (1)  can  be 
considered  as  one  of  the  methods  of  improving  images  by  conversion  of 
the  brightness  scale. 

The  result  of  conversion  of  the  i-th  element  of  the  input  image  for  the 
selected  linguistic  variable  is  determined  by  the  formula 


V=5|ljin(&i)i  (3) 

where  B  is  the  number  of  brightness  gradations. 

The  "brightness  cutoff"  of  the  initial  image  is  determined  by  using 

transformation  (3).  Adaptation  to  the  dynamic  range  of  brightness  is 

made  upon  formation  of  the  affiliation  function  zt _ (b)  for  one  or 

wm 

several  image  frames.  If  the  linguistic  variable  "DARK"  is  used,  the 
result  of  the  transformation  will  be  a  more  contrasting  image,  and  if 
the  "LIGHT"  variable  is  used,  the  result  will  be  an  inverted  contrast 
image. 

The  image  is  segmented  according  to  an  alphabet  of  n  linguistic 
variables.  The  k-th  variable  that  corresponds  to  the  maximum  to  it,  for 
which  the  following  condition  is  fulfilled 


pft  =  max|ij(&i),^  at  /=  1,  (4) 

is  then  selected  for  each  i-th  element  of  the  image.  This  linguistic 
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variable  is  characterized  by  a  brightness  code  Bk  and  is  calculated  by 
the  centroid  method  as  the  abscissa  of  the  center  of  gravity  of  the 
figure,  limited  by  the  affiliation  function  [5]. 


Figure  2.  Structure  of  Grammar  of  Concept 
"BRIGHTNESS" 


KEY : 

1 .  Brightness 

2 .  Light 

3 .  Gray 

4 .  Dark 

5.  Not 


6. 

Vary 

7. 

Lighter 

8. 

Darker 

9. 

More 

10. 

Less 

The  proposed  algorithm  for  segmentation  of  multigradation  images  using 
odd  linguistic  variables  satisfies  the  principles  of  the  completeness, 
connectedness,  homogeneity  and  other  requirements,  imposed  on  the 
segmentation  process  [1].  It  can  be  used  in  the  case  of  low-contrast 
images  with  contrast  noise,  under  conditions  of  fluctuation  and 
irregular  distribution  of  the  level  of  illumination,  along  with  other 
preliminary  processing  algorithms.  Segmentation  of  images  measuring 
128  x  128  elements  with  64  gradations  of  brightness  yielded  good 
results . 


The  advantages  of  the  given  algorithm,  compared  to  other  known 
algorithms,  include  the  possibility  of  adaptation  to  variation  of  the 
conditions  of  illumination,  greater  information  content  of  the  resulting 
image,  and  high  speed  and  simplicity  of  realization.  The  disadvantages 
may  include  some  degree  of  subjectivism  in  selecting  the  type  of 
affiliation  functions  and  of  the  alphabet  of  linguistic  variables  as  a 
function  of  the  image.  Studies  of  the  possibility  of  automating  this 
selection  are  now  under  way. 
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[Text]  Use  of  the  principles  of  associative  data 
processing  in  the  area  of  computer  vision  is 
considered.  A  method  of  realizing  associative 
concepts  for  pattern  recognition  is  proposed,  which 
includes  selection  of  the  method  of  encoding  the  image 
contours  of  objects,  and  formation  of  the  basic 
dictionaries  of  the  most  typical  elements  of  a 
silhouette  image  of  an  object. 

Recognizing  the  images  of  objects  is  a  process  of  analyzing  them  on  the 
basis  of  selected  features  and  working  out  decision  rules,  which  permits 
one  to  relate  each  of  the  imaged  objects  to  one  of  previously  formulated 
classes  according  to  a  specific  criterion.  The  theoretical  methods  of 
image  recognition  are  based  on  two  mathematical  approaches .  the 
discriminant  and  structure-linguistic  [1]. 

Discriminant  methods  of  recognition  comprise  the  basis  of  the  software 
of  serially  manufactured  robotic  vision  systems  (STZ),  which  are  used 
successfully  to  analyze  working  scenes,  containing  isolated  objects  that 
have  a  limited  number  of  stable  equilibrium  positions.  Complication  of 
the  scenes  and  objects  causes  specific  difficulties  in  the  use  of 
discriminant  methods,  since  not  only  the  numerical  characteristics  of 
the  objects,  but  also  information  in  terms  of  parts  and  relations  that 
characterize  their  mutual  location  must  be  taken  into  account  in  these 
scenes.  Structure— linguistic  methods  [2,  3],  based  on  the  structural 
relationships  between  "primitive"  and  their  ordered  combinations, 
correspond  to  the  greatest  extent  to  this  situation. 

However,  despite  a  number  of  successful  applications  of  structure- 
linguistic  methods  [1],  practical  realization  of  their  potential 
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capabilities  is  now  very  low.  This  is  explained  by  the  following 
factor:  specially  developed  methods  of  transforming  physical  objects  to 
mathematical  structures,  convenient  for  using  classification  algorithms, 
do  not  yield  a  clear  description  of  different  objects  [4] .  In  other 
words,  different  methods  of  structural  description  are  required  for 
different  classes  of  objects. 

Thus,  a  central  problem  in  structure  recognition  is  development  of  a 
universal  method  of  describing  the  shape  of  objects  (within  the  limits 
of  two-dimensional  projections  of  three-dimensional  objects)  using  the 
principles  of  associative  representation  of  information,  realized  by 
associative  memories.  The  use  of  dynamic  associative  memories  (DAZU), 
which  store  the  dynamics  of  input  events,  is  of  greatest  interest.  The 
distinguishing  feature  of  DAZU  is  the  capability  to  structure  the  input 
data .  Let  us  briefly  present  the  main  aspects  [ 5 ,  6 ] . 

Presentation  of  information  in  associative  memory.  The  DAZU  can  be 
described  by  a  network  of  parallel-connected  neuron-like  elements  (NE) 
with  time  summation  of  input  signals.  A  neuron-like  element  of  this 
type  (Figure  1)  is  a  device  containing  a  multibit  shift  register, 
threshold  converter  and  memory  element,  and  all  the  cells  of  the 
register,  having  exciting  or  inhibiting  input,  are  connected. 


KEY: 

1.  Threshold  element  2.  Memory  element 

The  first  part  of  the  NE  operates  as  an  address  part.  A  sequence  of 
pulses,  which  can  be  interpreted  as  the  sequence  A  =  (ai :ai  €  (0,  1)} 
(the  0  means  the  absence  and  1  means  the  presence  of  pulses  during  some 
time  cycle  i),  is  fed  to  the  input  of  the  neuron-like  element.  Symbol 
ai ,  arriving  at  moment  i ,  is  fed  simultaneously  to  all  cells  of  the 
shift  register  and  is  added  to  the  contents  of  each  cell  with  its  own 
weight  (+1  or  -1).  The  contents  of  the  cells  of  the  register  are  then 
shifted  by  one  position  to  the  right  (according  to  Figure  1).  The 
maximum  sum  stored  in  the  last  element  of  the  register  may  reach  a 
value,  equal  to  the  sum  of  the  ones  in  the  address,  and  this  occurs  only 
if  the  input  sequence  repeats  the  distribution  of  signs  at  the  inputs  of 
the  register . 

The  dynamic  associative  memory  contains  2n  neuron-like  elements,  each  of 
which  simulates  one  of  the  vertices  of  an  n-dimensional  single  hyp)ertube 
in  the  signal  space. 
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The  binary  sequence  A,  the  input  sequence  for  the  dynamic  associative 
memory,  is  transformed  to  the  n-dimensional  address  space  in  the  form  of 

A 

a  sequence  of  points — trajectories  A  =  F(A),  from  which  the  initial 

A 

sequence  A  =  F" 1  (A)  can  be  restored.  Transform  F  has  the  property  of 
associative  access  to  the  information  by  association  with  respect  to 
content  (to  n  symbols  of  the  input  sequence).  The  trajectory  is  stored 
by  introducing  the  memory  function  H  at  the  vertices  of  the  hypercube 
(the  next  symbol  of  the  sequence  is  stored). 

Allotment  of  function  H  by  threshold  properties  permits  one  to  realize 
statistical  processing  of  the  information  stored  in  the  dynamic 
associative  memory,  including  formation  of  a  dictionary  of  statistically 
reliable  events  {)?}  at  the  input  of  the  DAZU  in  the  form  of  a  set  of 
sequences  {A} : {^}  =  F({A}). 

Depending  on  the  threshold  of  function  H,  the  words  {/(}  will  either  be  a 
combination  of  trajectories  of  realizations  of  events  (all  the 
information  is  stored) ,  or  will  be  a  diffuse  or  point  intersection  of 
them  (a  more  or  less  general  part  of  the  information  is  allocated) .  The 
coinciding  sections  of  the  sequences — events — are  transformed  to  the 
same  chain,  and  different  sections  are  transformed  to  different  chains, 
thus  forming  the  graph  word. 

The  information  contained  in  dictionary  {[))  is  eliminated  from  the  input 
sequence  by  introduction  of  transform  F  * 1 ;  only  the  connections  of  the 

words  of  the  dictionary  are  retained  in  this  case,  and  a  syntactical 
sequence  C  =  Fc',[F(a),  {/?}  ]  is  formed.  The  mechanism  of  formation  of 

the  syntactical  sequence  permits  one  to  use  the  dynamic  associative 
memory  for  structured  data  processing  within  the  framework  of  the 
hierarchical  structure  from  the  DAZU,  since  it  becomes  possible  to 
combine  several  data  flows . 

Playback  of  the  information  using  transform  F" 1  permits  recognition  of 
the  input  sequence  by  comparing  it  to  the  reproduced  sequence. 

Coding  of  input  information  sequence.  The  selected  method  of  structured 
data  processing  is  realized  only  by  the  presence  of  an  active  method  of 
coding  it,  i.e.,  the  input  information  should  be  described  in  the  form 
of  a  sequence  of  codes.  Thus,  any  information  about  the  image,  for 
example,  contours  that  limit  the  zones  of  equal  elimination  or  zones  of 
identical  texture,  can  be  represented  by  this  method.  The  external 
contour,  a  description  of  which  is  also  selected  for  program  realization 
of  the  proposed  structured  processing  method,  carries  the  basic 
information  about  the  shape  of  the  object. 

Among  the  methods  of  contour  coding,  used  in  computer  vision, 
differential  Freeman  codes  are  most  suitable  for  solution  of  the  given 
problem,  since  the  information  about  the  contour  in  them  is  already 
represented  in  the  form  of  a  closed  series  of  symbols.  However,  this  is 
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inadequate,  since  coding  for  a  dynamic  associative  memory  should  be 
binary  (the  permissible  symbols  are  only  zero  and  one)  and  position- 
independent  . 

Fulfilling  the  first  requirement  permits  considerable  simplification  of 
the  methods  of  data  processing  in  a  dynamic  associative  memory,  and 
fulfillment  of  the  second  requirement  permits  a  significant  increase  of 
its  operating  reliability. 

The  binary  method  proposed  by  the  authors  describes  a  straight  line  in 
the  form  of  a  sequence  of  alternating  zeros  and  ones,  a  turn  to  the  left 
in  the  form  of  a  combination  of  two  zeros,  and  a  turn  to  the  right  in 
the  form  of  a  combination  of  two  ones ,  according  to  Freeman  coding . 

Formation  of  level-by-level  dictionaries.  The  following  sequence  of 
writing  information  to  the  hierarchical  structure  of  the  dynamic 
associative  memory  has  been  suggested  in  the  paper.  Code  sequences  that 
describe  a  set  of  straight  line  segments  are  used  as  the  elements  of  the 
lowest-level  (first-level)  dictionary.  The  arc  of  circles  of  different 
radius  are  included  in  the  dictionary  of  the  second  level.  The 
dictionary  of  the  third  level  contains  information  about  different 
angles.  This  representation  is  explained  by  the  fact  that  the 
statistically  more  reliable  elements  of  objects  (straight-line 
fragments)  are  primarily  allocated,  then  the  less  reliable  (arcs  of 
circles)  and,  finally,  more  variative  (angles)  are  then  allocated.  The 
proposed  representation  can  also  be  interpreted  in  terms  of 
two-dimensional  orthogonal  transforms,  where  the  process  of  detailing 
the  description  proceeds  from  the  lower  to  the  higher  harmonics. 

To  study  the  processes  occurring  in  the  dynamic  associative  memory  and 
the  methods  of  structured  processing  of  visual  data,  the  authors 
developed  a  research  complex  that  includes: 

a  machine  graphics  program  package  for  synthesis  of  the  images; 

an  image  preprocessing  program  package  that  converts  silhouette  images 
of  objects  to  contour  descriptions  (in  the  form  of  Freeman  catenary 
codes) ; 

a  program  for  recoding  a  catenary  code  to  a  binary  code; 
a  program  model  of  the  DAZU; 

a  program  for  visualization  of  the  trajectories  written  in  the  dynamic 
associative  memory. 

The  complexity  of  transforming  the  trajectories  in  an  n-dimensional 
signal  space  to  a  two-dimensional  visualization  plane  was  the  reason  for 
representing  the  graph  of  trajectory  transitions  in  the  DAZU  in  the  form 
of  a  tree,  in  which  branching  of  the  trajectory  in  the  signal  space  is 
depicted  in  the  form  of  a  fork,  while  the  fusion  of  two  trajectories 
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into  one  is  depicted  in  form  of  the  end  of  the  branch  of  the  tree. 
Continuation  of  the  trajectory  coincides  with  one  of  the  points  of  the 
tree.  The  branching  nature  of  the  tree  qualitatively  reflects  the 
degree  of  complexity  of  the  graph.  Thus,  if  the  graph  has  no 
intersections  (the  degenerate  case),  its  representation  will  be  a 
straight  line. 
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Figure  2 

An  example  of  teaching  the  DAZU  is  presented  in  Figure  2,  where  2,  a  is 
a  synthesized  silhouette  image  of  the  object,  and  also  its  contour 
representation  in  binary  code  (A  is  the  initial  point),  and  2,  b  is  the 
internal  representation  in  the  form  of  a  binary  tree • 
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Institute  of  Automated  Control  Systems  and  Radio  Electronics] 

[Text]  A  robotic  vision  system  has  been  proposed  for 
checking  the  shape  of  semiconductor  crystals,  realized 
on  the  universal  Elektronika-60  microcomputer.  A 
functional  diagram  of  the  visual  system  is  described. 

Algorithms  for  operation  of  the  system  and  its 
specifications  are  presented. 

Much  attention  is  now  being  devoted  to  development  of  robotic  vision 
systems  (STZ),  designed  to  automate  product  checking  operations  of  the 
electronics  industry  (semiconducting  crystals  and  plates).  One  of  the 
methods  of  realizing  STZ  is  to  develop  these  systems  on  the  basis  of 
universal  microcomputers  [1-3],  the  stock  of  which  contains  a  large 
number  of  various  types  of  machines. 

The  robotic  vision  systems  now  developed  on  the  basis  of  universal 
microcomputers  are  either  cumbersome  and  expensive  (made  in  the  CAMAC 
design)  [1],  are  inconvenient  to  service  (there  is  no  visual  check  of 
the  images  to  be  analyzed)  [3],  or  are  incapable  of  remembering  the 
necessary  number  of  reference  images ,  and  can  not  be  connected  to  a 
computer  network  [2] . 

The  robotic  vision  system  described  in  the  article,  designed  for 
automatic  quality  control  of  manufacture  of  semiconducting  crystals 
(check  of  the  conformity  of  the  shape  (design)  to  a  standard  and  of  the 
presence  of  mechanical  damage  ( scratches ) ) ,  does  not  have  the  indicated 
deficiencies.  An  Elektronika-60  microcomputer  is  used  as  the  control 
computer  in  the  robotic  vision  system.  An  Electronika  MS0585  or 
Elektronika  MS1212  microcomputer  can  also  be  used  in  the  STZ. 
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A  functional  diagram  of  the  robotic  vision  system  is  presented  in  the 
Figure.  The  STZ  consists  of  the  following  devices:  an  optical 
microscope  (M) ,  television  camera  (TK) ,  buffer  storage  (BZU), 
programmable  ROM  (PPZU) ,  electronic  disk  (ED) ,  long-range  communication 
module  (BDS),  an  interface  (I),  a  controlled  power  source  (UIP),  a  key 
module  (BK),  lighting  module  (BO),  color  monitor  control  module  (MUTsM) , 
a  color  monitor  (TsM) ,  black-white  monitor  (ChBM)  and  a  computer.  The 
optical  connections  are  shown  on  the  diagram  by  dashed  lines,  while 
electric  connections  are  shown  by  solid  lines. 
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The  total  operating  cycle  of  the  STZ  contains  two  steps:  teaching 
(writing  the  image  of  the  reference  object  to  the  memory)  and  operation 
(obtaining  the  necessary  data  on  the  object  to  be  checked).  The 
teaching  step  can  be  fulfilled  either  by  writing  the  information  to  the 
memory  of  the  STZ  (to  the  electronic  disk)  by  direct  reading  of  the 
image  of  the  reference  object,  or  by  recording  the  reference  information 
(to  the  electronic  disk  or  PROM)  by  program.  The  first  method  of 
teaching  is  used  in  the  STZ,  which  begins  with  establishement  of  the 
necessary  lighting,  provided  by  the  interface  (a  1-2  device),  controlled 
power  source,  key  module,  and  lighting  module  (assembled  on  light  diodes 
of  type  AL107A)  in  the  working  zone  of  the  microscope  (of  type  PME-1  or 
MBS-9).  The  light  module  contains  different  numbers  of  light  sources, 
the  light  intensity,  operating  frequency  and  frequency  of  switching-on 
of  which  are  established  from  the  computer.  The  light  sources  are 
arranged  about  the  control  object  ( crystal ) .  An  image  of  the  crystal  in 
the  working  zone  of  the  microscope  is  perceived  by  the  television 
camera,  is  read  (written)  to  the  buffer  storage,  which  stores  one  frame 
of  the  image  measuring  256  x  256  and  16  brightness  gradations  (64K 
four-bit  words ) .  The  reference  image  is  then  copied  to  the  electronic 
disk  (with  capacity  of  2,000  kbits),  which  stores  eight  frames  of  the 
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images.  The  image  is  processed  (filtered  or  segmented)  when  the 
reference  information  is  written  to  the  electronic  disk  according  to  the 
algorithms  described  below.  This  permits  one  to  present  a  higher 
quality  image  and  to  enhance  considerably  the  accuracy  and  operating 
speed  of  the  STZ.  The  process  of  teaching  the  STZ  terminates  at  this 
point.  The  STZ  program,  developed  in  Assembler  language,  is  written  to 
the  programmable  ROM.  Communications  are  established  through  the  common 
bus  of  the  robotic  vision  system  and  central  computer  by  using  the 
long-range  communication  module  (of  the  BDS-60  type),  which  permits  one 
to  provide  problems  of  simulation  of  the  operating  algorithms  of  the  STZ 
or  program  support  of  the  STZ  if  necessary.  A  local  data  transfer  net 
is  used  for  this.  The  original  devices  of  the  STZ  are  the  buffer 
storage  and  electronic  disk. 

The  buffer  storage  operates  in  the  following  modes:  writing  the  image 
frame  from  the  television  camera  to  the  memory  module  (BP)  with  a 
simultaneous  check  of  the  image  introduced  to  the  black-white  monitor, 
reading  the  information  from  the  memory  module  to  the  computer,  and 
writing  the  information  from  the  computer  to  the  memory  module  (to  check 
the  processed  image)  to  the  computer.  The  buffer  storage  is  embedded  in 
the  design  of  the  microcomputer  and  occupies  one  complete  card.  The 
memory  module  of  the  buffer  storage  is  realized  on  16  565RUZ 
microcircuits . 

The  electronic  disk  operates  in  the  following  modes:  writing 
information  from  the  microcomputer,  regeneration,  and  reading  the 
information  to  the  microcomputer.  In  the  write  mode,  the  electronic 
disk  generates  all  the  necessary  signals  for  storage  of  the  line  and 
column  addresses  in  the  internal  registers  of  the  memory  microcircuits, 
and  also  for  writing  of  data  to  the  corresponding  group  of 
microcircuits .  After  the  write  cycle  is  completed,  the  electronic  disk 
switches  to  the  regeneration  mode,  which  is  achieved  by  accessing  each 
of  256  lines  of  all  frames  at  a  speed  of  2 *10" 6  s.  The  read  mode  is 
similar  to  the  write  mode  and  is  controlled  from  the  computer.  The 
electronic  disk  is  structurally  made  in  the  microcomputer  standard  (one 
half-height  card) .  The  main  memory  of  the  electronic  disk  is  assembled 
on  32  565RU5  microcircuits. 

The  software  of  the  robotic  vision  system  permits  the  operations  of 
checking  semiconducting  crystals  or  plates:  estimating  the  correctness 
of  shape  ( arrangement  of  the  structural  components  of  the  crystal ) , 
measurement  of  the  geometric  characteristics  (area,  linear  dimensions, 
for  example,  length  and  width,  and  perimeter  of  the  contour),  and 
determination  of  the  mechanical  damage  (the  presence  and  dimensions  of 
chips  and  scratches ) . 

One  operating  cycle  of  the  STZ  consists  of  the  following  steps. 

determination  of  the  boundaries  (the  distance  of  the  edges)  of  the  image 
frame; 
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preliminary  processing  of  the  image; 
removal  of  the  necessary  information. 

The  boundaries  of  the  image  are  determined,  achieved  by  eliminating  the 
distortions  on  the  edges  of  the  image  frame,  in  the  first  step  of 
functioning  of  the  robotic  vision  system.  In  the  second  step  of 
operation  of  the  STZ,  the  contours  of  the  image  are  filtered  and 
determined.  The  filtering  operation  is  used  to  remove  noise, 
interference  and  distortions  of  the  images,  which  occur  due  to 
imperfection  of  the  hardware  and  due  to  the  process  of  formation  of  the 
image.  This  operation  is  achieved  by  one  of  three  methods:  by  the 
threshold  method,  by  low-frequency  filtration  or  by  median  filtration 
[4].  The  threshold  method  is  based  on  measuring  the  brightness  of  all 
the  elements  of  the  image,  on  comparing  this  brightness  to  the  mean 
brightness  of  adjacent  elements  of  the  image,  and  on  replacing  the 
measured  brightness  by  the  mean  brightness.  Low-frequency  filtration  is 
achieved  by  reducing  the  initial  image  file  with  the  correcting  file. 
Such  methods  of  filtration  as  a  window  in  three  or  five  elements,  a 
cross  in  five  or  nine  elements,  and  the  weighted  mean  (along  the 
vertical  or  horizontal  sequentially)  in  three  or  five  elements  are  used 
with  median  filtration. 

The  operation  of  discrimination  of  contours  is  based  on  the  use  of  one 
of  the  following  methods  [4,  5]:  of  a  differential  Roberts  or  Sobel 
operator,  of  a  logarithmic  Walsh  detector,  of  formation  of  the  contour 
by  analysis  of  four  or  eight  adjacent  points  of  the  image. 

Selection  of  the  filtration  method  and  of  discrimination  of  contours  is 
determined  by  the  geometric  and  light  characteristics  of  the  crystal, 
and  by  the  required  speed  and  operating  accuracy  of  the  robotic  vision 
system. 

All  the  necessary  measurements  of  the  visible  upper  plane  of  the  crystal 
are  made  on  the  third  operating  step  of  the  robotic  vision  system.  The 
algorithms  of  the  STZ,  which  realizes  the  method  of  comparison  to  a 
standard,  consists  of  the  following  operations:  "review"  of  the  crystal 
plane,  measurement  of  the  length,  width,  perimeter,  area  and  singles, 
determination  of  the  correctness  of  the  ratio  of  the  design  (for 
example,  displacement  of  the  contact  surface  from  the  center  of  the 
surface  of  the  crystal ) ,  and  estimation  of  the  presence  and  value  of  the 
area  of  chips  and  scratches  of  the  crystal.  A  reference  file, 
formulated  during  teaching  of  the  robotic  vision  system,  is  used  to 
analyze  the  geometric  characteristics.  Estimation  of  crystal  defects  is 
based  on  the  correlation  method  of  comparison  of  images  [5]. 

The  following  characteristics  were  obtained  when  testing  the  robotic 
vision  system:  resolution  of  0.2  percent  of  the  size  of  the  image  frame 
(for  example,  the  resolution  is  no  worse  than  100  ftm  at  frame  size  of 
25  x  25  mm,  the  accuracy  of  estimating  the  linear  dimensions  (total)  is 
random  (mean  square  deviation)  and  systematic,  which  comprise  a 


measurement  error  of  0.1  percent  of  the  size  of  the  image  frame  at 
checked  size  of  the  crystal  (0.3  x  0.3)10  “  ro,  and  the  time  of  checking 
one  crystal  of  0.1-10  s  as  a  function  of  the  required  measurements  and 
complexity  of  the  shape  of  the  crystal .  The  described  variant  of  the 
robotic  vision  system  is  being  operated  successfully  at  TIASUR  [not 
further  identified]  and  at  other  USSR  enterprises.  The  cost  of  the  STZ 
(including  the  computer)  is  approximately  13,000  rubles. 

Thus,  the  described  robotic  vision  system  is  less  expensive  and  is 
distinguished  from  known  systems  [1-3]  by  better  characteristics  due  to 
the  capability  of  storing  a  sufficient  number  of  reference  images,  due 
to  connection'  of  the  STZ  to  a  computer  network,  and  due  to  convenience 
of  servicing  because  of  selecting  the  image  on  black-white  or  color 

monitors . 
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[Text]  The  working  and  operating  principle  of  a  fiber 
sorption  muscle  and  the  structure  of  a  one-step  drive 
system  are  considered.  The  law  of  control  of  the 
drive  is  synthesized.  Promising  areas  of  using  the 
module  are  indicated. 

There  are  now  two  parallel  creative  processes.  One  is  related  to  the 
evolutionary  development  of  nature,  unobtrusive  in  daily  life,  and  is 
the  original  fiber  universality  of  the  muscular-receptor  development  of 
organisms.  The  other  is  more  dynamic  and  is  expressed  in  the  results  of 
the  technical  activity  of  man.  The  observed  attenuation  of  the  uniting 
mutually  enriching  bionic  approach  in  the  engineering  of  machines, 
mechanisms  and,  specifically,  of  robots  and  manipulators,  is  explained 
by  the  absence  of  satisfactory  analogues  of  muscles,  receptors  and  so 
on,  which  are  possibly  the  key  linking  elements  for  the  compatibility  of 
the  natural  and  technical  approaches.  The  particularly  theoretical 
papers  on  development  of  bionic  approaches  acquired  the  form  of  specific 
technical  developments  with  synthesis  of  the  world’s  first  fiber 
sorption  electronic  muscle  in  the  Soviet  Union  in  1988. 

The  main  elements  of  a  one-step  controlled  module  are  kinematic 
anthropomorphic  links  and  an  antagonistic  pair  of  fiber  sorption 
muscles.  A  layout  of  the  working  principle  of  the  muscle  is  presented 
in  Figure  1 .  It  contains  a  cord  of  an  anisotropic  sheathing  1  in  the 
form  of  a  hollow  cylindrical  cord,  inside  which  there  is  a  closed 
elastomeric  space  2.  Adsorbent-adsorbent  fillers  3  and  4  are  located  in 
this  space.  The  adsorbent  3  is  at  the  same  time  a  thermocouple. 

Electric  contacts  5  and  6  are  hermetically  led  to  the  outside  from  the 
opposite  ends  of  the  sheathing  1.  The  surface  of  the  adsorbent  3  is 
heated  due  to  the  effect  of  electric  current  delivered  to  contacts  5  and 
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6.  As  a  result,  the  adsorbate  4  leaves  the  surface  of  the  adsorbent  3, 
switching  to  a  vaporous  state.  The  pressure  in  the  elastomeric  space  2 
is  increased.  The  filaments  in  the  internodes  of  the  anisotropic 
sheathing  1  are  rotated  and  the  muscle  is  contracted.  If  the  current  is 
switched  off,  the  adsorbent  3  is  cooled  by  convection  or  preventively, 
and  the  adsorbate  4  is  adsorbed  due  to  reversibility.  The  pressure 
drops  and  the  muscle  does  not  prevent  its  own  lengthening.  Thus,  the 
force  developed  by  the  muscle  is  determined  in  the  final  analysis  by  the 
pressure  of  the  adsorbate. 

Two  main  elements  can  be  distinguished  in  the  muscle  unit,  a 
micromechanism  of  a  sheathing  force  converter  and  a  sorption 
microengine.  The  sheathing  transforms  the  internal  energy  of  the  muscle 

to  force  contraction. 

Studies  of  the  parameters  of  the  cord  1  revealed  a  specific  dependence 
of  the  maximum  value  of  the  relative  longitudinal  deformation  on  the 
nature  of  representation  of  braiding  of  the  sheathing.  Specifically, 
the  following  calculating  function  is  valid  for  braided  sheathing. 


C0S»inax  \7  {  ‘<0,21, 

COS  «rnm  l\  c  '  /  , 


(1) 


where  e  is  the  derformation  of  the  braided  cord  of  the  muscle,  a  is  the 

angle  in  the  internodes  of  the  braiding  of  rhombic  elementary  cells ,  c 
is  the  width  of  the  initial  filament,  and  hR  is  the  length  of  the  side 

of  the  rhombus  of  the  elementary  cells  of  the  cord. 


For  bound  sheathing,  the  length  of  the  edge  seemingly  increases  by 
approximately  the  width  of  the  filament  which  is  oriented  from  above, 
while  the  value  of  the  internodes  becomes  twofold  greater  than  that  in 
braided  sheathing.  On  this  basis,  a  similar  function  is 


£b  max 


COS  gmav 
COS  &min 


(2) 


The  maximum  angles  in  the  internodes  of  the  rhombic  cell  of  the  cord 
layer  of  the  muscle  sheathing  were  established 


54°>a>24V.- 

The  principles  that  link  the  pressure  inside  the  muscle  in  the  static 
mode  to  the  deformations  and  forces  developed  by  it  in  the  absence  of  an 
external  load  are  expressed  by  the  formula 
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(4) 


q  =  K  =  (x  -  1  j  ( ' 1  +  A.  ^(r2  -  A.)  +  ( l  ^x)  ( 1  +  A-)"'  X 


where  is  the  thickness  of  the  elastomeric  volume  of  the  muscle,  x  is 

the  relative  variation  of  the  length  of  the  muscle  (specific  length),  E 
is  the  modulus  of  elasticity  of  the  material  of  the  elastomeric  volume 
of  the  muscle,  and  r2  is  the  external  radius  of  the  elastomeric  volume 
of  the  muscle. 

The  force  developed  by  the  muscles  is  related  only  to  compressive 
strain.  Therefore,  the  external  load  on  it  has  physical  meaning 
exclusively  in  the  direction  of  the  tension  of  its  sheathing.  The 
dependence  of  the  steady  internal  pressure  qp  on  the  external  load  P  on 

the  muscle  is  derived.  The  values  of  q  and  qp,  according  to  the 

principle  of  superposition,  determined  the  total  internal  pressure 


qz  =  q  +  qp, 


(5) 


where 


qp  —  P  sin2  drain  [  ( 1 — X2  COS2  drain)  Jir !2] 


(6) 


Here  rj  is  the  internal  radius  of  the  elastomeric  volume  of  the  muscle. 


Figure  1.  Working  Principle  of  Electronic  Muscle: 

1 — cord  of  anisotropic  sheathing;  2 — closed  elastomeric 
volume;  3 — adsorbent;  4 — adsorbate;  5,  6 — electric  contacts 

The  mechanisms  presented  by  formulas  (1)— (6),  which  in  combination 
permit  numerical -step  variation  of  the  main  parameters,  specifically 
permitted  calculation  of  the  sheathing  of  the  first  electronic  muscle. 

Reversible  sorption  microengines  in  electronic  muscles  essentially 
perform  the  role  of  direct  converter  of  electric  energy  to  mechanical 
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energy.  Adsorption-desorption  equilibrium  is  in  the  given  case  the 
initial  physical  phenomenon.  It  is  expressed  in  regulated  concentration 
of  materials  on  the  surface  or  of  solid  micropores  in  the  volume. 

Matter  of  the  solid  phase,  on  which  the  material  to  be  absorbed  is 
concentrated,  is  called  the  absorbent.  The  material  to  be  absorbed, 
present  in  the  gaseous  phase  outside  the  adsorbent,  is  called  the 
adsorbtive,  while  its  part  directly  concentrated  on  this  surface  or  in 
the  pores  is  called  the  adsorbate.  Construction  of  the  sorption 
microengine  is  closely  related  to  the  concept  on  micropores  as  on  the 
regions  of  space  in  a  solid,  comparable  in  dimensions  to  the  adsorbed 
molecules .  The  field  of  the  adsorption  space  includes  three  types  of 
interactions  in  the  given  case.  The  main  interaction  is  dispersion, 
manifested  between  nonpolar  molecules.  The  electrons  move  in  a  complex 
manner  about  nuclei  in  atoms  and  molecules.  The  dipole  moments  of 
nonpolar  molecules  over  average  time  are  equal  to  zero,  but  the 
instantaneous  value  of  the  dipole  moment  may  be  distinct  from  zero .  The 
instantaneous  dipole  creates  an  electric  field,  which  polarizes  adjacent 
molecules.  As  a  result,  there  is  interaction  of  the  instantaneous 
dipoles.  The  energy  of  interaction  of  nonpolar  molecules  is  the  mean 
result  of  interaction  of  these  instantaneous  dipoles.  In  physics, 
forces  of  this  class  of  interaction  are  called  London  forces.  A  common 
nature  distinguishes  them— they  can  occur  between  any  atoms  and 
molecules,  by  which  their  universality  is  also  explained.  An 
orientation  effect,  caused  by  dipole-dipole  interaction  due  to  the 
effect  of  Keesom  forces,  occurs  in  the  presence  of  constant  dipoles  in 
the  molecules  of  the  adsorbate.  Intermediate  between  the  dispersion  and 
dipole-dipole  type  of  interaction  is  the  induction  effect,  which  is 
manifested  between  the  polar  and  nonpolar  molecules  under  the  effect  of 
Debye  forces.  London,  Keesom  and  Debye  forces  are  mutual  attractive 
forces.  Determination  of  the  energy  parameters  of  electronic  muscles  of 
fiber-sorption  type  is  semi— empirical  in  nature.  An  adsorption  isotherm 
is  initially  constructed  for  the  selected  type  of  adsorbate-adsorbate 
filler  of  the  electronic  muscle  in  the  form 


B=f{q)  at  r=const, 


one  of  the  experimental  variants  of  which  is  presented  in  Figure  2. 


The  work  of  adsorption  A,  equal  to  the  differential  variation  of  the 
Gibbs  adsorption  energy,  taken  with  opposite  sign 


A=—AG=RT\n(qa(T'), 


(7) 


where  R  is  the  universal  gas  constant,  T  is  the  absolute  temperature, 
and  qs  is  the  tabular  value  of  the  saturated  vapor  pressure  above  the 
liquid  adsorbate,  is  then  determined.  The  corresponding  qo  and  To,  and 
also  the  differential  molar  work,  which  expresses  the  characteristic 
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adsorption  energy  Bo,  are  found  in  the  absence  of  adsorbent  at 
temperature  T.  The  maximum  adsorption  for  temperature  To  is  determined 
by  the  equation 

Bq=Boi  exp[ — y(T — Tj)], 

here  7  is  the  thermal  coefficient  of  maximum  adsorption. 


Figure  2.  Experimental  Isosteres  for  Different 
Fractions  of  Filling  of  Volume  of  Adsorbate-Adsorbent 
System  of  Electronic  Muscle: 

V — volume  of  pores  of  adsorbent;  Vo — total  volume 


Figure  3.  One-Step 
Antagonistic  Pair 

KEY: 
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Drive  System,  Containing 
of  Electronic  Muscles 


4.  Position  sensor  of  output 
element  of  drive 

5.  Computer  of  difference  of 
mean  and  instantaneous  power 
of  heating  of  muscle 
adsorbents 
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Substituting  the  calculated  values  of  Bo  and  Ao  into  the  thermal 
equation  of  Dubinin-Astakhov  adsorption 


the  adsorption  isotherms  are  calculated  for  a  number  of  values  T,  lying 
in  the  temperature  range  of  practical  interest.  The  graphs  of  the 
adsorption  isosteres  are  then  constructed,  i.e.,  functions  of  the 
following  type  are  constructed 

q=h(T)  at  Bo = const. 

Adsorption  increases  as  the  pressure  of  the  adsorbtive  increases .  The 
differential  work  on  a  background  of  an  increase  of  absorption  with  an 
increase  of  the  pressure  of  the  adsorbtive  decreases.  The  thermodynamic 
state  of  the  electronic  muscle  is  expressed  by  adsorption-desorption 
equilibrium  between  the  adsorbate  and  adsorbtive,  which,  as  in  the 
classical  variant  of  a  heat  engine,  is  determined  by  three  values: 
pressure,  temperature  and  the  volume  of  the  working  gas.  The  difference 
is  only  that  the  volume  in  the  given  case  has  adsorption  properties  7. 
Differential  variation  of  entropy  can  be  expressed  by  finding  the 
differential  of  Ag  by  temperature: 

«  —  TT-  =  -  * "> -  (8) 

—  RiMw-')-RT(-l£A-)+RT^p-. 

The  direct  derivative  is  correctly  written  in  the  last  term,  since  qs  is 
independent  of  adsorption  B.  Let  us  find  the  differential  variation  of 
enthalpy,  having  substituted  into  the  equation 


ah=ag+tas 


the  corresponding  expressions  from  equations  (7)  and  (8): 

=  —  RT\n{qsq'1)  -  RT\n(qsq  gr  — 
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The  first  term  of  the  right  side  of  equation  (9)  expresses  the 
differential  heat  of  adsorption  of  the  adsorbent,  while  the  second  term 
corresponds  to  the  heat  of  condensation  of  the  adsorbtive.  Thus, 
construction  of  isosteres  is  used  to  calculate  the  differential  entropy 
and  differential  enthalpy  of  the  sorption  microengine  of  the  electronic 
muscle . 

If  there  is  no  external  interactive  force,  the  length  of  the  muscle  is 
dependent  on  its  internal  thermodynamic  state.  Without  regard  to 
thermal  losses,  it  is  characterized  by  the  mean  power  of  the  heating  of 
the  adsorbent 


Nc=±j.Nd^,  (10) 

where  N  =  I2r  is  the  instantaneous  heating  power,  I  is  the  electric 
current  in  the  adsorbent,  r  is  the  resistance  to  current,  and  T\  is 
time. 

If  the  given  contraction  of  the  muscle  was  established,  one  can  in  the 
first  approximation  permit  a  link  between  length  la  and  power: 


la  =  KyNc+lo, 

lo  is  the  length  of  the  muscle  at  P  =  0  and  K  is  a  coefficient 

C  3 

determined  experimentally.  The  steady  force  P^,  developed  by  the 
muscle,  is  represented  by  the  formula 


Pm =Kz  (/.-/)  =  Kz  (K*NC + *o-l) , 


in  which  1  is  the  length  of  the  muscle  deformed  from  outside  and  is 
the  longitudinal  stiffness  coefficient.  With  regard  to  (10),  we  find 

The  structure  of  a  one-step  drive  system,  containing  an  antagonistic 
pair  of  muscles,  is  shown  in  Figure  3:  M|  and  M2  are  the  muscles,  U  is 
an  amplifier,  UU  is  the  controller,  DP  is  the  position  sensor  of  the 
output  element  of  the  drive,  and  VM  is  a  computer  of  the  difference  of 
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the  mean  and  instantaneous  power  of  heating  of  the  adsorbents  of  the 
muscles . 


For  an  antagonistic  pair  of  muscles  in  the  steady  state,  the  forces 
developed  by  them  are  equal  to 


here  Ii  and  I 2  are  currents  of  the  first  and  second  muscles  and  1\  and 
I2  are  their  lengths.  Let  us  assume  that  lo  -  li  =  1^,  lo  -  I2  =  1m 

1  £  0  are  the  position  of  the  output  element  of  the  drive.  The  values 

1  =  0  and  li  =  I2  =  lo  correspond  to  the  off  state  of  the  pair  of 

muscles . 

The  force  of  drive  P  is  equal  to  the  difference  of  forces  Pi  and  P2 , 
acting  oppositely.  In  the  steady  state 

P  =  KZ  I [*.-£  I  (/?  -  /!)  dz,  -  24  . 

In  the  transient  process  from  the  initial  (unheated)  state  of  the 
muscle,  P  is  calculated  by  the  equation 

P (z)  =  KZ |tf.  j  (/? -  II)  dzt  -  2/„ j  ,  (ll) 

where  t  -  T\  -  To  is  the  time,  as  a  function  of  which  force  P  developed 
by  the  drive  with  time  delay  varies,  and  r o  is  the  time  delay  in  the 
appearance  of  force,  corresponding  to  the  difference  of  the  mean  heating 
power  of  the  adsorbents  of  the  muscles  and  to  the  position  of  the  output 
element  of  the  drive  1  at  moment  of  time  r i . 

In  (11),  ~dx=dx\\  /i=/i (ti),  h=h(xi),  4=4( t).  The  forces  of  currents  I, 

and  I2  are  related  to  the  value  of  signal  Us  t  at  the  input  of  the 
amplifier 


I\  —KyUstku  II—  KyU8i%2i 
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where  Ky  is  the  transfer  coefficient  of  the  direct  channel  of  the  drive, 
/li  is  a  coefficient  equal  to  unity  at  Ust  >  0  and  equal  to  zero  at  Us  t  < 

<0;  ^2  =  1  _  ^  i  • 

Let  us  substantiate  the  law  of  control  Us  t »  realized  in  the  controller 
(UU).  The  signal  Us  (Figure  3)  corresponds  to  a  given  position  of  the 
output  element  of  the  drive  1B ,  i.e.,  to  the  positions  of  the  mechanical 
link  Lg  =  Knlg  at  =  const. 

The  different  processes  of  the  tasks  Us(r)  are  found,  having  minimized 
by  function  ln ( r )  the  following  functional,  known  in  optimal  control 

theory 


/-  (12) 

Then  varying  the  value  of  coefficient  K&.  Here  1B  =  lB(r)  and  lB<l>(r) 
are  given  functions,  and  L(1,(t)  =dl,/dx,  /n(1)(T) T3‘dtjd,z.  jf  t*-*-  oo,  l, — q( 
ls(l)—ln0) -*■  0,  the  necessary  condition  of  a  minimum  of  (12)  will  be  an 
Euler  equation  [2] 


(/.(»— /„('))-  0.  (13) 

Let  us  find  the  relationship  of  the  speed  of  the  output  element  of  the 
drive  ln(  1  *  to  power  Nn  and  N^.  To  do  this,  let  us  first  assume  that 

there  is  no  longitudinal  deformation  of  the  muscle  by  external  force  P^. 

This  means  that  the  force  of  the  drive  P  is  equal  to  zero.  Then 

(/?-/?)  *. 

Having  differentiated  ln( r )  by  variable  r,  we  find  the  desired  function 


(14) 


in  which  ^/12__/aa^ r  is  the  instantaneous  power  and  H nc ™  (7js — 

-  I22 )dr  is  the  mean  power.  For  the  position  mode  of  development  of  the 


113 


task  Us  ( r) ,  let  us  assume  that  derivative  ls ( 1 ) ( t)  is  equal  to  zero. 
After  substitution  of  (14)  into  (13),  we  find 


*  (4  -  /.)'+ (AT«  -  Nn)  =0. 


(15) 


The  relation  1  -  1  is  insignificant  at  small  values  of  r.  If  r  -*  oo, 

s  n 

the  difference  Nnc  -  Nn  can  be  disregarded. 

The  law  of  control  of  the  drive  is  derived  by  combining  (15)  with  the 
equation  of  motion  of  the  output  element  of  the  drive 


mlnw  =  P—Pw' 


where  m  is  the  reduced  mass  of  the  structural  elements,  displaced  by  the 
drive,  Pw  is  the  force  of  resistance  to  the  drive,  and  P  is  the  force  of 

the  drive  in  the  transient  process.  However,  if  the  indicated  equations 
are  combined,  we  arrive  at  a  control  law  Us  t ,  to  realize  which  is  it 
possible  only  by  having  determined  Pw  and  time  derivatives  of  ln.  This 

is  very  difficult  in  practical  realization. 

The  approach  to  control  of  the  drive  using  only  an  Euler  equation,  valid 

at  P„  =  0,  appears  more  rational.  In  this  case,  the  law  of  control  is 

w 

formulated  so  that  the  drive  system  approaches  the  state  described  by 
equation  (15)  at  each  moment  of  time  r.  Let  us  assume  that  this 
equation  is  not  exactly  observed.  Then  either  a  positive  or  negative 
value  is  found  in  its  right  side.  In  both  cases,  the  right  side 
approaches  zero,  if  the  control  signal  is  sufficiently  large  in  absolute 
value  |  Us  1 1 ,  and  its  polarity  corresponds  to  the  direction  of 
development . 

Let  us  introduce  into  consideration  the  control  law 


Uat  =  j  F{o)dv, 


where  q=Tn  (/,_/„)  +  J (Nnc—Na) ,  1  • e  ■  >  of  the  left  side  of  equation 

(15);  variable  r  is  replaced  by  =  const,  r n  is  the  time  segment,  on 

which  the  transient  process  of  operation  of  the  drive  essentially  ends 
with  given  accuracy,  F(^)  =  U>  Oat  $  >  0,  F(*r)  =  U<  Oat  <  0, 

F( (r)  =  0  at  <r  =  0,  and  |u|  =  const. 
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Replacing  r  by  constant  means  that  increased  sensitivity  to  the 

difference  of  lengths  1  -  1  is  given  to  the  drive  system. 

Implementation  of  the  integral  control  law  Ust(r)  creates  smoothness  in 
development  of  the  task  by  the  drive  system  and  contributes  to 
elimination  of  autooscillations  in  the  steady  process.  For  direct 
adjustment  of  the  drive,  it  is  recommended  that  the  proposed  control  law 
be  used  in  the  form 


T  .  ! 

Ust  —  j  U  dx  at  cr=  U — In + K (Nac — N„)  >  0; 

Ust=  /  ( — U)dx  at  a<0;  . 

o 

T 

Ust=  f  0  dx  at  c=0, 

£/=const>0. 


(16) 


Adjustment  of  the  drive  reduces  to  finding  coefficient  K. 

Let  us  assume  that  Us  t »  after  integration,  reaches  saturation  U  >  0  at 

H 

<r  >  0  and  1  -  1_  <  0  and  that  it  satisfies  the  condition  N_  <  N„^.  K 

s  n  n  nc 

then  increases  due  to  the  appearance  of  a  large  difference  of  the 

intensities  of  currents  I\  and  I2.  As  a  result,  1^,  which  approaches 

ls,  increases  with  a  delay  of  r.  The  difference  |N^c  -  |  and 

|lg  -  ln|  decreases  in  absolute  value  and  0  approaches  zero.  Initially 

ff  -  0  at  some  small  value  of  |l'g  -  1^1  and  then  <r  <  0.  The  process  of 

a  decrease  Us  t  begins ,  and  a  decrease  of  power  corresponds  to  this . 

Thus,  power  will  increase  on  the  final  section  of  the  transient 

process  when  1^  -»  1  ,  and  will  decrease  due  to  a  decrease  of  Us t .  As 

a  result,  N  -  "*  ^ ’  allc^  a  0*  At  the  end  of  the  process,  in  the 

stable  system,  =  N^,  Us t  -  const,  and  lg  =  1^.  The  transient 

process  should  proceed  in  similar  fashion  if  r  =  0  and  tr  <  0  and  1  - 

-  !n  <  °. 

One  can  judge  the  stability  of  the  drive  by  equation  ( 16 ) .  Let  us 
assume  that  the  mean  power  N  does  not  vary  significantly  during  the 

nc 

transient  process,  a  delay  occurs  after  moment  of  time  r  =  to. 
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The  equation  of  dynamics 


2KzU + mlaw = KyrK,Kz  V\t-Pw 


will  then  have  the  solution 


<„(')  =  [-£*--  U« +  /« (°)]  ■ '  “s  ]/- 


2  K, 


m 


(17) 


w 


U-—2t 


j-j-  follows  from  (17)  that  the  system  is  stable f  but  autooscillations  are 

possible.  The  delay  To  almost  disappears  for  adequately  heated  muscles. 

Each  muscle  approximates  an  aperiodic  link.  The  controller  changes  Ust> 

maintaining  the  equality  1_  =  1  .  The  amplitude  of  the  autooscillations 

fl.  s 

of  heated  muscles  reduces  to  a  minimum  or  disappears  by  selecting 
coefficient  K  in  formula  (16). 


Figure  4.  Experimental  Curves  of  Transient  Processes 

1^( t )  of  Drive : 

1 — process  of  development  of  task  at  =  0.6  N  and  at 

initial  weakened  state  of  muscles;  2  similar  process 

at  P  =  1  N;  3 — process  of  development  of  task  by  heated 
W 

muscles  at  =  1  N 

The  experimental  curves  of  the  transient  processes  ln( r )  of  the  drive 
are  presented  in  Figure  4. 
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Drives  based  on  fiber-sorption  muscles  with  range  of  developed  forces 
from  0  to  40  N  have  been  manufactured  since  1990.  The  efficiency  of  the 
first  models  of  their  drive  is  approximately  1  percent.  The  time 
constant  of  the  drive  for  unheated  muscles  in  the  initial  state, 
considered  as  an  aperiodic  link  (without  delay) ,  is  at  the  level  of 
30  s.  Work  is  now  under  way  to  increase  the  speed  of  the  muscles  and  of 
their  efficiency. 
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[Text!  The  construction  and  design  principles  of 
multioutput  actuating  modules  with  several  independent 
linear  and  angular  motions  are  considered.  A 
simulation  model  of  a  two-output  actuating  module, 
designed  on  the  basis  of  Petri  nets,  is  presented. 

When  developing  industrial  robots  and  manipulators,  it  is  necessary  to 
ensure  simplicity  of  design  and  minimum  values  of  mass-size  indicators 
[1,  2].  One  of  the  methods  of  solving  this  problem  is  the  use  of 
multioutput  actuating  modules  (MIM)  in  the  actuating  subsystem  of  the 
robot,  which  permit  one  to  produce  a  set  of  independent  linear  and 
angular  motions  at  the  output.  MIM  can  be  located  directly  in  the 
joints  of  the  structural  layout  of  the  robot ^/or  on  a  fixed  base. 

The  design  principle  of  the  multioutput/  actuating  module  is  based  on 
supplying  the  moving  parts  of  the  motor  with  couplings,  controlled 
separately  from  the  controller  [3],  A  generalized  block  diagram  of  the 
MIM  is  presented  in  Figure  1.  The  MIM  includes  a  controller, 
reciprocating  motor,  2n  (n  >  2)  working  couplings,  n  fixing  couplings, 
mechanically  connected  to  the  reciprocating  motor,  linking  elements, 
made  in  the  form  of  nonmagnetic  rods  and  flexible  cables,  and 
transforming  members,  made  in  the  form  of  pulleys  connected  to  the 
moving  links  of  the  robots. 
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Figure  1.  Generalized  Block  Diagram  of  Multioutput 
Actuating  Module: 

UU — controller;  DVPD — reciprocating  motor;  EKS — fastening 
and  joining  elements;  M — couplings;  SPR — degrees  of 
freedom  of  robot;  PO — converting  member;  DP — position  sensor 


KEY: 
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The  controller  controls  the  reciprocating  motor,  2n  working  and  n  fixing 
couplings,  and  is  connected  to  the  robot’s  control  system.  The 
reciprocating  motor  can  be  based  on  electromagnets  or  pneumatic  and 
hydraulic  cylinders  or  other  similar  devices,  and  is  used  to  ensure 
reciprocating  motion  of  the  working  couplings.  The  working  and  fixing 
couplings  are  controlled  separately  from  the  controller  and  transmit 
reciprocating  motions  through  the  linking  elements  either  directly  or 
through  transforming  members  to  the  movable  links  of  the  robot.  The 
designs,  shapes,  dimensions  and  operating  principles  of  the  working  and 
fixing  couplings  are  dependent  on  the  shapes,  dimensions  and  materials 
of  the  linking  elements,  where  the  latter  are  used  to  transmit  the 
motions  of  the  couplings  to  the  degrees  of  freedom  of  the  robot  and  are 
different  shafts,  rods,  and  flexible  elements  (cables). 

The  linking  elements  and  together  with  them  the  corresponding 
electromagnetic  couplings  can  be  located  on  the  outside,  inside,  or  a 
combination  as  a  function  of  the  kinematic  layout  of  the  MIM  and  the 
mutual  spatial  arrangement  of  the  degrees  of  freedom  of  the  robot. 

The  transforming  member  is  used  to  convert  the  linear  motions  of  the 
elements  and  is  kinematically  connected  to  the  moving  links  of  the  robot 
manipulator. 

The  multioutput  actuating  module  operates  in  the  following  manner. 
Producing  several  angular  and  linear  motions  is  based  on  conversion  of 
reciprocating  motions  of  two  moving  parts,  moving  in  the  opposite 
direction  in  n  (n  >  2)  independent  linear  and  angular  motions  of  the 
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rods  and  cables  by  alternately  switching  the  groups  of  2n  working  and  n 
fixing  couplings  on  and  off,  controlled  separately  for  each  independent 

motion.  The  functions  of  converting  X  to  xj  (i  =  1,  n)  are  established 
by  the  laws  of  variation  of  the  signals  of  the  controller  i ,  2 >  • • • * 

Un  over  time. 

As  a  result,  the  reciprocating  motions  of  the  reciprocating  motor  are 
transformed  to  a  combination  of  linear  (xt,  x2  ,  . . . ,  xi),  and  angular 
Ul  +  ,»  »1  +  2 .  •  •  • ,  <»n)  motions,  which  are  transmitted  by  transfer 
mechanisms  to  the  moving  links  of  the  robot  manipulator  or  of  the 
robotic  complex  (RTK) . 

The  controller  of  the  MIM  is  based  on  a  microprocessor  complex  of  series 
K580  and  generates  control  signals  for  the  reciprocating  mo  or  an 
windings  of  the  electromagnetic  couplings  according  to  a  given 
algorithm.  The  control  algorithm  of  the  MIM  is  based  on  the  phase 
engagement  method,  which  permits  uniform  variation  of  the  dynamic  force. 
The  controller  of  the  multioutput  actuating  module  is  linked  to  the 

control  system  of  the  robot  and  generates  control  signals  U  =  OJi ,  U2 , 
...,  Un)  according  to  deviation  of  the  given  values  of  vector  q3  from 

the  actual  values  of  q^. 

Petri  nets,  which  are  distinguished  by  clearness  of  the  model  in 
combination  with  rather  broad  functional  capabilities,  have  been  used 
extensively  of  late  among  the  set  of  simulation  models,  oriented  toward 
solution  of  various  problems  of  simulation,  analysis  and  synthesis  o 
processes  of  different  physical  nature. 

The  Petri  net  (SP)  is  formally  determined  as  a  set  of  five  of  type  [4] 

'  :  5  =  <P,  T,  I,  0,  po>,  ^  (1) 


where  P,  T,  I,  and  0  are  finite  sets  of  transitions  (events),  positions, 
input  positions,  and  output  positions,  respectively ,  /to  is  the  initial 
marking?  and  fi3: P  -»  N,  N  =  (0,  1,  2,  ...}  is  a  set  of  natural  numbers. 

Problems  of  design  of  simulation  models  of  multioutput  actuating  modules 
of  robots  and  robotic  systems  are  considered  in  the  paper  on  the  basis 
of  the  apparatus  of  expanded  Petri  nets,  designed  for  description  and 
analysis  of  the  functioning  of  the  latter. 
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Figure  2.  Multioutput  Actuating  Module  With 
Two  Degrees  of  Freedom 


KEY: 

1 .  Multioutput  actuating  module 

A  multioutput  actuating  module  of  the  electromagnetic  type,  designed  for 
horizontal  and  vertical  motions  in  the  structure  of  the  robot  or  robotic 
system,  is  simulated.  The  MIM  includes  an  electromagnetic  module  with 
movable  elements,  a  coupling  module,  specifically  electromagnetic 
couplings,  which  are  divided  into  working  and  fixing,  depending  on  the 
functional  designation  in  the  MIM,  a  controller  and  output  rods. 

Let  us  consider  the  design  step  of  the  interpreted  Petri  nets  on  the 
example  of  a  multioutput  actuating  module  with  two  degrees  of  freedom 
(Figure  2).  One  degree  of  freedom  is  realized  by  using  a  rod  I,  which 
performs  horizontal  displacement  to  the  left  (to  the  right),  and  the 
second  degree  of  freedom  is  realized  by  using  a  flexible  rod  II  (or 
ca^->le) »  which  performs  vertical  motions  downward  (upward) . 

The  positions  of  the  rods  are  fixed  along  the  vertical  and  horizontal  by 
using  position  sensors  in  coordinates  Yl,  Y2  and  XI,  X2,  respectively. 
One  of  the  variants  of  the  model  that  describes  the  functioning  of  the 
multioutput  actuating  module  is  the  Petri  net,  shown  in  Figure  3. 

Positions  1  and  2  describe  the  starting  conditions  of  the  electromagnets 
of  the  MIM,  3,  4,  and  5,  6  reflect  the  states  of  these  electromagnets, 

,  8  and  9,  10  reflect  the  states  of  the  couplings  which  are  rigidly 
connected  to  the  moving  elements,  11  and  12  are  introduced  to  simulate 
the  processes  of  control  of  the  electromagnetic  couplings,  13  and  14 
reflect  the  state  of  the  output  rod,  15-17  are  used  to  simulate  the 
processes  of  control  of  the  motor  with  regard  ot  the  information  of  the 
position  sensor,  and  18  simulates  stopping  of  the  rod.  Transition  t| 
reflects  the  process  of  variation  of  the  state  of  the  electromagnet, 
transition  t2  reflects  the  condition  of  connecting  the  coupling,  t3 
means  transition  of  the  coupling  from  one  state  to  another,  and  t4 
reflects  transition  of  the  rod  from  one  position  to  another. 
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Figure  3.  Simulation  Model  of  Multioutput 
Actuating  Module 


The  presence  of  markers  at  positions  13  and  14  denote  the  response  of 
the  sensor  at  position  XI,  and  the  presence  of  markers  at  positions  15 
and  16  denote  a  response  of  the  sensor  at  position  X2.  The  marker  at 
position  17  denotes  response  of  the  sensor  at  position  of  either  X  or 
X2,  while  a  marker  at  positions  4  and  5  denotes  motion  to  the  left  or'  to 
the  right,  respectively.  The  delay  time  of  the  marker,  equal  to  the 
time  of  displacements  between  extreme  positions,  corresponds  to  these 
positions.  Let  the  markers  be  located  at  positions  15-17.  The 
electromagnets  of  the  multioutput  actuating  module  can  be  started  only 
after  a  marker  appears  at  position  1.  Upon  response  of  transition  ti ,  a 
marker  appears  at  position  4,  which  means  motion  of  the  movable  element 
to  the  left.  Markers  arrive  at  positions  13,  14,  and  17  after  response 
of  transitions  t2 ,  t3 ,  and  t4 .  The  resulting  marking  is  the  initial 
marking  for  motion  of  the  rod  to  the  right. 


The  work  of  the  second  output  of  the  multioutput  actuating  module  can 
also  be  simulated  by  a  Petri  net,  presented  in  Figure  3.  In  this  case, 

positions  I,  3,  4,  7,  8,  13,  and  14  simulate  the  motion  of  the  output 

rod  downward,  while  positions  2,  5,  6,  9,  TO,  15,  and  16  simulate  motion 
of  the  output  rod  downward. 

The  resulting  Petri  net  of  the  multioutput  actuating  module  is  formally 
described  in  the  following  manner: 


P=  {Pu  Ps,.:.,  Pib},  T=  {*j,  t2, . . .,  4};  •  • 
1(h)  =  {Pi,  Pb,  Pis,  Pi*}.'.  0(h)  =  {P4};  : 
I(t2)  =  {Pi},  0(t2)  =  {P1},\I(t3)  =  {Pu  Pn}; 

0 (t»)  =  {Pa},  I(U)  =  {Pa},  0 (U)  =  {P,3;  Pu,  Pn) \ 
I (*»)  —  {P2,  Ps,  Pis}  >  0  (4)  —  {Pa} ,  1(h)  ~  {Pa}  ‘, 
/(4)  =  {P9,  P12},  0(t7)  =  {P1(j); 

I  (ts)  =  {P.o} ,  0(4)  =  {P .5,  Pie,  P17} : 

-  /(4)  =  {Pio},  0(4)  =  {Pia}. 
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New  designs  of  electromagnetic  and  pneumatic  multioutput  actuating 
modules,  used  as  the  actuating  members  of  robots  and  robotic  complexes, 
have  been  proposed  on  the  basis  of  the  generalized  block  diagram  of  MIM. 
The  proposed  designs  are  distinguished  from  existing  robot  drives  by 
simplicity  of  design  and  minimal  mass-size  indicators. 

The  structural  layouts  of  industrial  robots,  operating  in  different 
coordinate  systems,  have  been  worked  out  on  the  basis  of  multioutput 
actuating  modules.  The  proposed  designs  of  MIM  can  be  used  in 
construction  of  the  actuating  members  of  adaptive  robots,  robotic 
complexes,  and  of  systems  of  different  designation. 
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